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ABSTRACT 


Several Important aspects of the Regulator with a Prescribed Degree of 
Stability (RPDS) methodology and Its applications are considered. The 
solution of the time-varying RPDS problem as well as the characterization 
of RPDS closed- loop elgenstructure properties are obtained. Based on the 
asymptotic behavior of RPDS root-loci, a novel one-step algorithm for 
designing Regulators with Prescribed Damping Ratio (RPDR) Is developed. 

The robustness properties of RPDS are characterized In terms of the pro- 
perties of the return difference and the Inverse return difference matrices 
for the RPDS state feedback loop. This class of regulators Is found to 
possess excellent multi- loop margins with respect to stability and degree 
of stability properties. The ability of RPDS design to tolerate changing 
operating conditions and unmodelled dynamics are also Illustrated with a 
multi-terminal DC/ AC power system example. The output feedback 
realization of RPDS requires the use of Linear-Quadratic-Gauss Ian (LQG) 
methodology. Several Interesting Issues that arise In the application of 
robustness recovery procedures to the design of RPDS based LQG compensators 
are also examined. 
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CHAPTER I 
INTRODUCTION 


1.1 Thesis Motivation 

A linear time-* Invariant (LTI) system of the form 

i(t) =* A xCt) (1.1) 

is said to have a degree of stability a H all the poles of A are located 
to the left of the line a * - oi (where a is real and positive) in the left- 
half complex plane (see Figure 1.1). The Regulator with a Prescribed 
Degree of Stability (RPDS) problem is one of determining the weighting 
matrices 2 . and R in the cost functional of a Linear Quadratic (LQ) regxilator 
problem. 

[x^(t) S x(t) +^^(t) R u(t)Jdt (1.2) 

. '^o 

subject to the dynamic constraint 

x(t) = Ax(t) + B u(t) (1.3) 

so that the resulting state feedback design has a prescribed degree of 
stability. 

The RPDS problem for LTI systems was first formulated and solved 
by Anderson and Moore [An IJ ; their work is one of the many attempts to 
combine the poleplacement techniques and LQ methodology in some useful 
fashion. Other major contributions in this area include the asymptotic 
LQ poleplacement technique of Harvey and Stein [Ha 1] and the sequential 
poleplacement method of Solheim [So Ij . Both of these methods address 
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FIG. 1.1 Region (shaded) of Allowable Closed-Loop Poles 
for a RPDS Design with Stability Factor Equals 
to a . (a > 0) 


-3- 


the problem of finding the appropriate ^ and R matrices that correspond 
to a prescribed set of eigenvalues. 

The RPDS problem Is useful from the application point of view. In 
many engineering problems of Interest, the exact position of the closed- 
loop poles may be of secondary Importance. As a result, the designers 
are Interested only In locating all the closed-loop poles In some regions of 
the left-half complex plane. 

One such example which motivates the research In this thesis Is 
the design of a state feedback law for a reduced order model of a m-termlnal 
DC/AC power system([Gr 1], [Ch 1]). This is the preliminary stage of 
design for a decentralized output feedback control scheme. A standard 
power system model of this type, with n aggregated areas, has (2n-2) 
modes of Iterarea oscillations^. Such modes are usually poorly damped 
In the absence of compensation. The two remaining modes are real. The 
one located at the origin Is the clock error mode. The one located 
slightly to the left of the origin is the average frequency mode. I''or 
physical reasons. It Is desirable for both of these modes to remain 
unchanged under state feedback. This could be trivially accomplished 

2 

In the LQ design by making such modes unobservable In the cost functional. 

The primary design objective Is to find a state feedback control law that 
results in sufficient damping for the (2n“2) oscillatory modes. One 
natural way to specify the damping criterion Is simply to require all the 

^ See Figure 5.1 of Section 5.1 for an example of such pole configuration 
2 

A detailed discussion of such a poleplacement technique Is given In 

Section 5.2 > 

i 

t 
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closed-loop poles to lie to the left of some line <7 • in the complex 
plane. This class of problem can be effectively handled using non- 
asymptotic LQ poleplacement techniques such as RFDS and Solhelm's method 
[So 1]. With respect to the computational requirements, RFDS Is the 
simpler of the two methods. In every design Iteration using RFDS method, 
one needs to solve a Rlccatl equation and an eigenvalue problem. In the 
case of Solhelm's method, one may have to solve as many as n Rlccatl 
equations and n eigenvalue problems for each Iteration. Although the 
above mentioned design objective can also be handled by asymptotic 
techniques such as Harvey/Stein's method [Ha 1] with relative ease; the 
use of which In the multltermlnal DC /AC power system Is however Inappro- 
priate. From an engineering point of view. It Is undesirable to speed up 
some of the oscillatory nodes as required by the asymptotic poleplacement 
method In order to allow the poles of the remaining modes to approach their 
specified locations. 

Despite the potential usefulness of RFDS as evidenced In the above 
discussion, research on this class of regulators has been largely over- 
looked In the literature. Apart from a few brief remarks found In the 
original work of Anderson and Moore [An 1] , little Is known to date 
about the relation between the prescribed degree of stability and the 
feedback properties of a RFDS system. This Is partly a result of the 
lack of appropriate tools for effective analysis of multivariable feed- 
back systems. It Is only recently. In the context of studying the robust- 
ness properties of controllers derived using Llnear-Quadratlc-Guasslan 
(LQG) technique that an appropriate formulation has emerged ([Le 1], [Do 4]). 
One major objective of this thesis Is therefore to analyze the feedback 
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properties of RPDS with special emphasis on robustness properties. 

Anott;er Important problem of Interest here Is the possible 
application of RPDS to the design of LQ regulators with specifications 
other than prescribed degree of stability. One such example Is the 
design of regulators with prescribed damping ratio (RPDR) . This class 
of LQ regulators has all Its poles lying Inside a convex cone centered 
at the origin and s}nnmetrlc about the negative real axis (see Fig. 1.2). 
Such designs are of potential application In the power system example 
mentioned above. 

It Is hoped that the new Insights on properties of RPDS obtained 
In this thesis will help the control system designers to better appreciate 
the benefits that they may expect from this class of regulators. 

1.2 Thesis Organization 

This thesis Is organized as follows: 

In Chapter II we first formulate and solve the generalized RPDS 
problem for linear time-varying systems. A precise notion of degree of 
stability that applies to all finite dimensional linear system Is 
Introduced for this purpose . We then specialize these results to the time- 
invariant case which was the original form of RPDS problem addressed by 
Anderson and Moore. Extension of RPDS techniques to the dual problem of 
designing a Kalman Bucy filter with a prescribed degree of stability 
(FPDS) Is also considered, and Interpretations of the FPDS design procedure 
In the context of estimation are given. 

In Chapter III, the elgenstructure properties of time- Invariant RPDS 
are examined. We first derive some sensitivity formulas for the closed- 




FIG. 1.2 Region (shaded) of Allowable Closed^Loop Poles 
for a Regulator with a Prescribed Damping Ratio 
CO80 (0 < 9 < j) 
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loop poles with respect to the stability factor cx • This Is followed 
by a study of the behavior of RFDS root-loci as the relative weight 
between the state and control weightings varies. It turns out that this 
problem can be formulated In a fashion similar to the conventional LQ 
root-locus problem [Kw 1], The control weighting matrix Is chosen to 
be of the form PR. The relative weight between the state and control 
weightings can be adjusted by varying the scalar p . We are particularly 
Interested In the asymptotic behavior of the root- loci as p(p > 0) becomes 
very large or close to 0. Based on the asymptotic properties of the 
RPDS root-loci as p approaches Infinity, a novel algorithm for solving 
the regulator with prescribed damping ratio problem Is derived. 

The robustness properties of RFDS are discussed In Chapter IV. Adopting 
the framework developed In [Do 2] and [Le 1], we use the minimiim singular 
value of the return difference and the Inverse return difference matrices 
as basic robustness measures for multiple-input multiple-output (MIHO} 
feedback systems. The dependence of the robustness properties on the 
stability factor a is analyzed by using the frequency domain Inequalities 
for the RFDS return difference and Inverse return difference matrices. 

The dual robustness problem for FFDS Is next formulated, and the FFDS 
robustness properties are Interpreted In the context of estimation. 

In Chapter V, we consider an example where the RFDS technique Is 
used to design a state feedback control law for a reduced order model 
of a multi-terminal DC/AC power system. Besides validating theorems and 
conjectures stated In previous chapters, the results obtained here also 
reveal several desirable features of RFDS designs not predicted by the 
theorems developed In this thesis. This provides considerable justification 
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£or the use of RPDS methodology In power syeteos control. 

Since full stato. feedback can never he realized In real life, we 
are thus motivated to study the feedback properties of RPDS based LQG 
compensators In which the feedback gain is designed by using the RPDS 
technique and where a KBF is used to provide the state estimates. This 
problem is considered in Chapter VI. LQG compensators are known to 
possess no guaranteed uncertainty tolerance in general [Do 4] . As a 
result » we only focus our discussion in this thesis on the class of RPDS 
based LQG compensators that are designed with the robustness recovery 
procedures ([Do 1] and [Rw 1]). Such procedures allow one to approximate 
the state feedback transfer matrix (and consequently their robustness 
properties) with the LQG loop-*transfer matrix In a systematic manner. 
Several difficulties encountered In attempts to combine the RPDS technique 
and robustness recovery methods In designing LQG compensators are 
Illustrated with numerical examples. 

Chapter 7 consists of conclusions and suggested directions for 
future research. 

1.3 Contributions of This Thesis 

The major contributions of this research are: 

(1) Formulation and solution of the tlme-'varylng RPDS problem. 

(2) Development of a simple procedure for studying the asymptotic 
behavior of the RPDS root-loci. 

(3) Development of a novel algorithm for solving the regulator with 
prescribed damping ratio problem. 

(4) Clarifying the robustness properties of RPDS 



(5) Uncovering some pocenclal problems chat one may encounter 
In combining the RPDS technique with robustness recovery 
methods In designing LQG compensators. 
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1.4 Notation 


SISO 

single- Input slngls-output 


MDfO 

multiple- Input multiple-output 


LQ 

llner-quadratlc 


LQG 

linear-quadratic Guasalan 


KBF 

Kalman Bucy filter 


KPDS 

regulators with e prescribed degree of stability 

FPDS 

Kalman Bucy filter with a prescribed degree of stability 

LTI 

linear time- Invariant 


T(») 

loop transfer’ matrix 


T^U) 

loop transfer matrix for RFDS with a 

degree of stability a 


loop transfer matrix for FFDS with a 

degree of stability a 

L(») 

multiplicative perturbation of X(>) ot 

\<A) 

an eigenvalue of A 



complex conjugate transpose of A ^ 



maximum singular value of A ^ ^nax^ 

1 

mlnufflum singular value of A ^ ^isln^ 

(i* A) 


(/ A) 

4 

defined as 


A > B 

A - ^ is a positive definite matrix 


A > B 

A - £ is a positive seml-definlte matrix 

1 

Identity matrix 
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CHAPTER II 

THE REGULATOR WITH A PRESCRIBED DECREE OF 
STABILITY PROBLEM AHD ITS DUAL 

2«1 Introduction 

In this chapter, w« formulate and solve the RPDS problem for finite 
dimensional linear system. By 'linear system' we mean a pair of 
equations of the form 

x(t) • ^(t) x(t) + u(t) (2.1) 

^(t) - C(t) x(t) (2.2) 

We make an assumption at this point that Is to hold throughout this 
chapter. The elements of the mtrlx A(t) are continuous and bounded 
functions of time defined on -« < t < • . The elements of the mtrlces B(t) 
and £(t) are piecewise continuous and bounded functions defined on -«kc<m . 

The formulation and solution of the time-varying RPDS problem Is 
given In section 2.3. In order to study this problem, a precise definition 
for 'degree of stability' that applies to both time-varying and time- 
invariant system Is Introduced. The time- Invariant version of the RPDS 
problem studied by Anderson and Moore [An 1] becomes a special case of 
our formulation. 

The problem of designing Kalman Bucy Filter with a prescribed degree 
of stability (FPDS) Is considered In Section 2A . In view of the dual 
nature between the Kalman Bucy Filter problem and the LQ regulator problem, 
the design techniques for RPDS can be readily applied to FPDS. Several 
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Inctrprecaclons of the fornulatlon and aoluclon of che 7PDS problam are 
given. 


2.2 Degree of Stability for Linear Sviteaa 

The concept of 'Degree of Stability* for linear tlae- invariant (LTI) 
systems is directly '.'elated to the system pole locations. 

Definition 2.1 An autonomous system x(t) ■ A x(t) is said to possess a 
degree of stability a, for some real positive constant a, if all the 
eigenvalues of A are located to the left of the line o ■ -a in the left* 
half complex plane. 

Intuitively speaking, a LTI system with a degree of stability oi will 
attenuate any given initial state perturbation at a rate faster than or 
equal to e~^^. To motivate the extension of such concepts to the time- 
varying case, let us first recall the following definition of exponential 
stability for finite dimensional linear system. 

Definition 2.2 ([Br 1 1), Section 29) 

The system x(t) ■ A(t) x(t) is said to be exponentially stable if 
there exists some positive constants Y and X , such chat for all c 
and t^ in the half plane t > t^. we have 



ll<t,t^)! ^Y« 

(2.3) 

where 

• 



Mt,t^) - A(t) £(t,t^) 

(2.4) 

and 


d(t ,t ) ■ I 

“ O 0 ~ 

(2.5) 
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It is well known that Che exponential stability of LII systems 

can be easily established by checking the pole locations of the system 

matrix. The following theorem clarifies the relationship between the poles of 

a LTI system x(c) ■ A x(c) and the Induced norm^ of Che respective 

At 

transition matrix e— . 


Theorem 2.1 Consider the autonomous system 
x(t) » A x(t) 


( 2 . 6 ) 


Let a be some given positive constant. Then 


max (Re(\.(A)) < -ol 


if and only if there exists a positive constant Y such that 


A(t-t ) 
o I 


< Y ® 


-a(t-t^) 


for all t > t 


(2.7) 


( 2 . 8 ) 


Proof ; Sufficiency ; Suppose that (2.8) is true for some positive y and 

Re(X^(A^ 


a given ci. Observe that e 


is the spectral radius 


A(t-to) 

s 

always less than or equal to its spectral norm, we get 


of e ° . Using the fact that the spectral radius of a matrix is 




A(t-t ) 

O ! 




< Ye 

Since Che inequality (2.9) holds if and only if 


(2.9) 


All the vector norms considered here are the Euclidean norm. The 
corresponding induced norm is Che 2-macrlx norm CSee(De IJ). 
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max Rea.(A))<-a , (2.10) 

i 

this completes the proof of safflclency. 

Necessity ; Suppose that max(Re^.(A)) < -a . We can find 

i ^ 

a similarity transformation that reduces A to Its Jordan form 
Aj * P. A jP ^ ([B 1] Section 12). For system in Jordan Form, the transition 
matrix is block diagonal with each of the nonzero element being one of 
the three forms 


k C7t ^ ..h ^t . .k Xt 

t e cos cut, t e sin cut or t e 


( 2 . 11 ) 


where a- Jcu and X are eigenvalues of A. It follows from our assumption 

on the eigenvalues of A that each element of e^^ approaches 0 at a rate 
-^t 

faster than e as t approaches infinity. 

/S 

This in turn implies the existence of some positive constant Y 
such that 


), 

T o 


Y-e 


-a(t-to) 


( 2 . 12 ) 


Since 


A - P'^ Aj P , 
A(t-t ) 

lie- ‘’ll 


it follows that 



° P 


1 Lt (t-t ) 

i iip^iiii.^ °iiiipii 

1 -®(t-t ) 

< Y llP'^l! Hill e 


(2.13) 


This completes the proof of necessity. 



-15- 


Corollary 2«1 The autonomous system (2.6) has a degree of stability oi 

If and only If there exists a positive constant Y such that for all 
t > t 




Y e 




(2.14) 


Proof ; This Is a trivial consequence of Definition 2.1 and Theorem 2.1. 


It seems natural from the result of Corollary 2.1 that we should 
define the stability of a linear system ° A(t) x(t) In terms of the 
norm of Its transition matrix. 

Definition 2.3 

A linear system x(t) ■ A(t) x(t) Is said to have a degree of stability 
a If there exists some positive constants Y such that for all t ^ t^, 

-a(t-t ) 

lli(t.t^)il 1 Y e ° (2.15) 

The following theorem characterizes the degree of stability using 
Lyapunov's Direct Method. 

Theorem 2.2 Let V(x(t),t) be a Lyapunov function of the form 
T 

V(x(t),t) ■ X (t) 5. i(t) for x(t) » A(t) x(t) on the whole state space. 

If for some constant e > 0, 

Q(t) > e I for all t > t (2.16) 

— — — O ) 

then the system described by x(t) » A(t) x(t) has a degree of stability a 
If for all t ^ t^. 


V(x(t),t) < - 2a V(x(t),t) 


(2.17) 
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Proof ; It follows from (2.17) that V(x(t),t) is bounded above by 
e Since V(x(t),t) llx(t)ll^ this means that 

/. i\ -a(t-t ) 

I lx(t) 1 i <(^ V(x(t^) je (2.18) 

Dividing both sides of (2.18) by I 1 1 taking the supremum 

on the right, we obtain 


liuTI ' 




(2.19) 


Now, by taking the supremum on the left side of (2.19), we 


obtain 


||£(t,to)|| £^^Amax (£(t^)))^ e 


( 2 . 20 ) 


which establishes the degree of stability. 


Remark In situations such as designing feedback laws with optimal control 
methods where a Lyapunov function is available. Theorem 2.2 provides 
a convenient mean of evaluating the degree of stability for a feedback 
design. Tnls point will be made clear in Theorem 2.4. 

2.3 Formulation and Solution of the Continuous Time RPDS Problem 

In this section, we formulate the RPDS problem for finite dimensional 
linear systems. Although the RPDS problem studied by Anderson and Moore 
is a special case of that considered in this thesis, their solution of 
time- invariant problem extends readily to time-varying case. 
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2.2.1 The RPDS Problem Statement 

Consider the following LQ regulator problem 


min J 
u(t) 


11m I X 

t I 
1 


^(t) £(t) xCt) +u^(t) R(t) u(t) dt 


( 2 . 21 ) 


subject to the dynamic constraint 


x(t) » A(t) x(t) + Mt) ^(t) (2.22) 

x(t ) * X (2.23) 

— o*^ -o 

where (A(t), ^(t)) Is uniformly completely controllable. 

The RPDS problem Is to find the appropriate weighting matrices £(t) 
and R(t) such that the steady state control law that minimizes the cost 
functional (2.21) has a degree of stability a. 


2.2.2 Solution of the RPDS Problem 

The solution of the RPDS problem Is given In the following 
theorem . 

Theorem 2.3 Consider the deterministic LQ regulator problem In 
Section 2.2.1. The weighting matrices are chosen as 

a(t) - £ e^“^ (2.24) 

R(t) - R (2.25) 

where *£(t) and R(t) e^^*" are piecewise continuous and bounded on 

[t^,«] with 

a(c) i i and R(t) 1 ^2 I 


for all t (2.26) 
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1 

2 

where 6^ and 3^ are some positive constants. Furthermore, (^(t) , A(t)) 
is uniformly completely reconstructable. Then the solution ^(t) of 
the Rlccatl equation 

“Vt) - (A(t) + a K^(t) + K^(t) (A(t) + a I) + £(t) 

- K^(t) B(t) R"^(t) B^(t) £^(t) (2.27) 

with terminal condition ^ converges to JK^(t) as tj^ approaches “ 

for any ^ ^ £ and the steady state optimal control law given by 

u(t) - -r^(t) B^(t) K^j^(t) x(t) (2.28) 

has a degree of stability ex. 


Proof ; The assumptions on A(t), ^(t) , £(t) and R(t) relating to 
continuity, boundedness, uniform complete controllability and uniform 
complete reconstructability guarantee the existence of an exponentially 
stable steady state optimal control law that minimizes the given cost 
functional. To establish the degree of stability of the control law, 
we Introduce the following transformations 


i(t) 

£ 

a(t-to) 

e 

i(t) 


A 

a(t-t ) 


u(t) 

m 

O 

e 

u(t) 


Then, x(t) and u(t) are related dynamically by 

dt 


> . <^(t-t ) 

x(t) • jr (.^ x(t)) 


“ (A(t) + a 1) x(t) + B^(t) u(t) 
The initial condition of x(t) is given by 


(2.29) 

(2.30) 

(2.31) 

(2.32) 
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x(t ) - X 
— o — o 

The cost functional J In (2.21) can now be written as 


(2.33) 




[x(t)^ £(t) x(t) + u^(t) R(t) u(t)]dt (2.34) 


Observe the strong connection between the transformed LQ regulator 

problem described by (2.32) and (2.34) and the original problem described 

by (2.21) and (2.22). Suppose that ^(t) Is the optimal control at time 

t for the original problem, then ^(t) ■ e^ u (t) Is the value of the 

optimal control at time t for the transformed problem. The resulting 

value of the state x*(t) Is given by x*(t) « x (t), provided that 
Qtt . 

ft Q it 

X (to)** e X Thus, a feedback control law obtained for the 

transformed problem readily yields a feedback control law for the original 
problem. Moreover, the resulting minimum value of the cost functional 
Is the same for each problem. 

Our next step Is to study the transformed LQ regulator problem 


min J 
u(t) 


11m 




’^(t) 2(t) x(t) + u'^(t) R(t) u(t)] 


dt 
(2.35) 


subject to dynamic constraint 


x(t) » (A(t) + a i) x(t) + ^(t) ^(t) ; i(tg) • Xq (2.36) 

We need to check all the technical conditions that ensure the 
existence of an exponentially stable steady state optimal control law 
for the above optimization problem. It is trivial to see that A(t) -4*a^ 
Is bounded and continuous. Moreover, the piecewise continuity and 
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boundedness of ^(t) , £(t) and R(t) follows from the assumption of the 

original LQ regulator problem. To guarantee the existence of an 

exponentially stable steady state control^law» we also require the 

2 

uniform complete reconstructability of (£ (t) , A(t) + a I^} . The uniform 
complete controllability of (A(t) a B^(t)) follows from that of 

(A(t), ^(t)) by observing the equivalence of the following four state- 

, 1 
ments. 


(1) (A(t), B^(t)) is uniformly completely controllable 

(11) There exists some positive numbers such 

that 


(a) Y, 


i -k 


^(t„ +a,T)B(T) b^(t) (ji^Ct^ +C.T) dr i 


(2.37) 




I J ■‘'<7*''^) B(t) B^(T) £^(t^+ 0,T)JdT 


£ 1 i *11 tc 


(2.38) 


(111) There exist some positive real numbers C7, Y > Y-, > B , B. 

o 1 o 1 

such that 

to-Hj 


(a) Y. 


■ 4 [. 


a(t +a-T) 


l(t +<7,T) B(T) r(T) 4>^(t -W.T) 


a(t^+a-T) 


]dT Yj^ I for all t<. 


(2.39) 


Readers are referred to Chapter I of [Kw 1] for a detailed discussion 
on uniform complete controllability and uniform complete observability 
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(b) 3 I < ♦(t ,t +0) 

o — — — o o 


L 0 

B(t) b\t) O.T) e ° JdT 


f^o^r a(t + 


a-T) 


+ o) e"°^ < 6. I 
— 0 0 — 1 •- 


for all t. 


(2.40) 


(Iv) (A(t) + a ^,B(t}) is uniformly completely controllable. 

The equivalence between (111) and(lv) follows from the observation 
a(t-t ) 

that ^(t,t^) e 
described by (2.32). 


Is the state transition matrix of the system 


By duality, the uniform complete recons tructablllty of (j^‘‘(t), A(t)) 
also Implies the uniform complete recons tructablllty of (£^(t),A(t) +al). 

Thus, It follows from Theorem 3.6 of [Kw IJ that ^^(t) , the 
solution of the Rlccatl Equation (2.27) with terminal condition 
K^Ctj^) ■ 1^, converges to jCQj(t) as tj^ * for any K ^ 0. Moreover, 

K^(t) Is also a solution of (2.27). The resulting steady state optimal 
control law is given by 

u(t) — R"^(t) B^(t) K^(t) x(t) (2.41) 

It follows from the uniform complete observability of 

1 

2 

(2. (t) , A(t) + al) that the closed- loop system 

x(t) - (A(t) + - B(t) R‘^(t) ^"“^(t) £a^t))x(t) (2.42) 

is exponentially stable. Applying the feedback law (2.41) to the original 
system (2.22), we get 


x(t) - (A(t) - B(t) V^(t) B^(t) K^(t))x(t) 


(2.43) 
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Let 6 (t,t ) and 6 (C,t ) be the stete transition natrlces of 
-ca o *c o 

systems (2.42) and (2.43) respectively. Since system (2.42) Is exponen< 
tlally stable, there exists some positive constants X and Y such that 
for all t > t , we have 




(2.44) 


Since ) e 

.-C o 


a(t-t J 


<^^(t,t ), this means that 
_cu o 




.-a(t-to) 


(2.45) 


'X(t-t ) -a(t-t-) 


^Ya ” e 

which establishes the prescribed degree of stability of system (2.43) 
and this completes the proof. 


Remark The proof of the above theorem Is parallel to that given by 
Anderson and Moore [An 1] for the time- Invariant version of the theorem. 

It Is also possible to prove Theorem 2.3 without resorting to a trans- 
formed LQ regulator problem. Consider the Riccatl equation of the 
minimization problem defined In the statement of the theorem 
• 

K(t) - -ii(t) A(t) - A^t) K(t) + e K(t) 3(t) R(t) ^ B(t)^ K(t) 

- S(t) e^“*^ (2.46) 

Let us postulate a solution of the form 

K(t) - K (t) e^^ (2.47) 

where (t) Is some positive definite matrix function. Substituting 
~a 

(2.47) into (2.46), we obtain 
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-^(t) • Kjit) (A(t) +ap + (A(t) +ai)^ Ka(t) +g(t) 

- K^(t) B(t) R“^(t) B^(t) ^(t) (2.48) 

which is the saoe Rlccatl equation aa (2.27) obtained from the trans- 
formed LQ problem. If we further assume that the boundary conditions 

- ~ 2at 

at t^ for equations (2.48) and (2.27) match, then K(t) ■ ^ 

(t ,t. ) where K (t) Is the positive definition solution of (2.27). As 

0 1 “tl 

t^ approaches we have l^(t) approaches ^(c) for any 

- - 2at 

l^(tj^) * i 2.* H«nce K(t) approaches J^(t) e for 

K (tj^) • ^ It follows that the feedback control law Is 

given by 

u(t) - -(R(t) e^“^)“^ B^(t) K (t) x(t) 

” ~Cl “ 

- -i'^(t) B^(t) K (t) x(t) (2.49) 

-u 

which agrees with that given In Theorem 2.3 and this establishes our 
claim. 

An alternate way to verify the degree of stability for the closed- 
loop system 

i(t) - (A(t) - B(t) R"^(t) B^(t) l^(t))x(t) (2.50) 

Is to use the result of Theorem 2.2. The following theorem Identifies 
the Lyapunov function required for such purpose. 

T 

Theorem 2.4 V(x(t),t) • x (t) l^(t) *(t) is s Lyapunov function 

for the closed-loop system (2.50) with the property that 
V(x(t),t) 1 -2a V(x(t),t). 
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Proof: To show that V(jc(c),t) is a Lyapunov function, wa have to 

establish the following: 

(i) V(x(t),t) is once differentiable with respect to time 

(ii) There exists some positive constants ir^ , ir^ and 8 such that 


and 



(2.51) 


3V(i,t) 

3t 



for all t > t (2.52) 

0 


(iii) V(x(t),t) < 0 along the trajectory of 

£(t) - (A(t) - B(t) R‘^(t) B^(t) ^(t)) xCt) 

for all t > t (2.53) 

— o 

For any given x e r“, (x, t) • x^ io(t) x. Since ^(t) is the 

steady state solution of the Riccati equation (2.27), it follows that 
V(£,t) is once differentiable with respect to t. By Theorem 3.4 of 
[Kw 1], we have 


V(x,t) - 11m X 
t,-« “ 


where 


1 


t) l£(T) +^(T) B(X) R'^^j ^^(t) 1^(t)1 


L 


(2.54) 


^ ^(x,t) - (A(t) + C^I - B(t) R"^(t) B^(t) ^(T))^(T,t) (2.55) 

Since x(t) - (A(t) + - B(t) R’^(t) B^(t) Koj(t)) x(t) is exponen- 

tially stable by problem definition, and 

[2(t) B’^(t) B^(t) ^(t) ] is a bounded positive definite 
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macrlJt function by the propnrtie* of the tine-varying Riccati equationa, 
it then foUowa from Theorem 6 in section 3.1 of I&r IJ that there exists 
constants and such that for all t and x c R°, 

0 < 1*1^ 1 V(at,t) < lil^ C2.56) 

3V 

To show the boundedness of ^ (x,t) » take the norm on both sides 

of the Riccati equation (2.27) and obtain 

lli^t)|| < 2(a+l|A(t)l|) |(K/t)|| + l!(i(t)l| 

+ I 11^1 li(t) 1 1 ^onixi (R(t) ) for all t > t^ (2.57) 

Since all the aatrices on the right hand side of the above inequality 
are bounded from above and since ^(t) is also bounded from below as 
assumed in Theorem 2.3, it follows that there exists a constant 6 such 
that 

H<i*t)|< llK^(t)l| e 1x1^ (2.58) 

To cootpute the derivative V(x(t).t) along the trajectory of 
i(t) - (A(t) + a I, - ^(t) R~^(t) B^(t) ^(t)) x(t), observe that 
V(x(t),t) - J(t) [^(t) A(t) - l^(t) B(t) R‘^(t) B^(t) I^(t) 

+ A^(t) K (t) - K (t) B(t) R'^(t) B^(t) K (t) 

— -rx -o — ” "" -a 

+ K^t)]x(t) (2.59) 

- x“^(t) 1-2 l^(t) - 2(t) - J^Ct) B(t). i'^tt) B'^(t) l^(t)] x(t) (2.60) 

< - 2olx^(t). K^(t) x(t)] (2.61) 


20 V(x(t),t) < 0 


(2.62) 
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In going from (2.59) to (2.60), w« us«d the Rlccaci cquaclon (2.27). 

The Inequality (2.61) la obtained by obaervlng that the natrlx function 
(a(t) B(t) R"^(t) B’^(t) K^(;)) Is positive seal-definite for 

all t ^ t^. Inequality (2.62) follows from the positive definiteness 
of V(x,t) and the assusqttlon that a la positive. Hence V(x,t) Is a 
Lyapunov function for x(t) “ (A(t) - B^(t) R ^(t) B^(t) ^(t))x(t) 
with Che property V(x(c),t) ^ -2a V(x(c),c) and this completes the proof. 

Remark By Theorem 2.2, the condition V(x(t) ,t) < - 2 aV(x(c) ,c) Implies 
that the system x(t)«(A(.t) - B(t) R*^(t) B^(t) 5j^(t))x(t) has a degree of 
stability a. This agrees with the result given In the proof of 
Theorem 2.3. 

Remark The result of this theorem is not new. The LTI version of 
Theorem 2.4 was used by Anderson and Moore [An 1] to establish the 
asymptotic stability of LQ regulators. Safonov has employed the same 
type of Lyapunov function to Induce conic sectors Chat were used to 
analyze the stability robustness of LQ regulators. 

If Che matrices A(t), B(t), ^(c) and R(c) in Tl-x^ irem 2.3 are time- 
invariant, then the resulting state feedback control law Is also time- 
invariant. Observe that the cost functional is still tlms'^arylng because 
It includes the factor e^^. Tlxat the feedback law should be constant 
In this case Is not at all obvious. The proof of this fact simply 
follows from that of Theorem 2.3. Applying the transformation (2. 29) 
and (2.30) to the cost functional 
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lim 

t,-** 


/;■ 


2o(t-t^) 




^(t) £x(t) + (/(t) Ru(t)j 


dt 


(2.63) 


and Che dynamic constraint 


x(t) “ A x(t) + B u(t) 


we get 
J » 

and 


^ t n 

3 ,1 j ^ ^x(t)^ £ x(t) +u"^(t) Ru(t)jdt 


x(t) ■ ( A + al) x(t) + ^ u(t) 


(2.64) 


(2.65) 


( 2 . 66 ) 


respectively. 


Since the transformed minimization problem Is time- Invariant, the 
resulting feedback law for this problem Is of the form 

(2.67) 


u(t) - -r"^ ^ x(t) 


where ^ Is the unique positive definite solution to Che algebraic 
Rlccatl equation 

^ (A + ap + (A + al)^ ^ B 5,"^ ^ + a “ £ (2.68) 

The controllability of (A + a^, p , which Is Implied by that of 

(A,p , ensures the existence of . The stability of the control law 

~l/2 

(2.67) follows from the observability of (^ » A + oi) which Is Implied 

~l/2 

by that of (£ ,A) . Reversing the transformation (2.29) and (2.30) 

we obtain the optimal control law for the original problem 

-a(t-t ) ^ ) 

u(t) - -e Sa ® 


“ -R"^ K^^x(t) 


(2.69) 
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whlch Is time- Invariant. That the closed- loop system 
• -IT 

x(t) ■ ( A - B R B x(t) has a degree of stability a follows from 
* -1 T ~ 

the stability of x(t) ■ (A + al - B R B K ) x(t). In terms of the 
closed- loop eigenvalues, we have 

Re(X^(A + al - B r"^ B^ I^)) < 0 
which in turn loq>lies that 

Re (X^(A - B r"^ B^ ^)) < - a . (2.70) 

and this establishes the degree of stability of the feedback lai: (2.67). 

We summarize the results of the above discussion in the following 
theorem 


Theorem 2.5 Consider the LQ regulator problem 


min 

u(t) 

^^o 


11m 




(t) £(t) x(t) + u^(t) 


R(t) u(t)J 


dt (2.71) 


subject to the dynamic constraint 


x(t) ■ A x(t) + B u(t) 


(2.72) 


where (A, ^) is controllable. 

The weighting matrices are chosen as 

£(t) - £ and R(t) » R e^^ (2.73) 

where £ > 0 and R > 0. 

~l/2 

Furthermore, supposed that (£ , A) is observable. Then the 

optimal control law for the given LQ regulator problem is given by 


(2.74) 


uCt) - -R-^ ^ Kc^(t) 

where ^ Is the unique positive definite solution of the algebraic 
Rlccati equation 

K (A + al) + (A + al)’^ K - 0 (2.75) 

Moreover, the closed-loop system 

i(t) • (A - B r"^ ^) x(t) (2.76) 

has a degree of stability at least a. 

In view of the tlnif^ invariant control law that results from the 
minimization problem considered in Theorem 2.3, one might suspect the 
possibility of constructing a cost functional with time-invariant matrices 
£ and R such that the control law obtained from minimizing the cost 
functional is the same as that derived in the preceding theorem with 
2 e^“^ and R e . Pairs of weighting matrices £, R with the above 
mentioned properties do indeed exist. The following corollary to 
Theorem 2.5 shows how such matrices may be chosen. 

Corollary 2.5 If S_(t) and R(t) in Theorem 2.5 are chosen to be 

2(t) - + 2a (2.77) 

and 

R(t) -R (2.78) 

respectively, where ^ is the unique positive definite solution of the 
algebraic Riccatl equation (2.73) > then the resulting optimal feedback 


law is given by 
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• -IT 

and the closed- loop system x(t) * ( A - ^ r ^ x(t) has a 
degree of stability at least a . 

Proof ; Note that (A, B^) Is controllable by assumption and that 
■" 1/2 

(C5, + 2a ,A) Is observable as a result of ^ > 0. 

It follows that the algebraic Rlccatl equation 

K A + A^ K + (£ + 2a - K B K « 0 (2.80) 

associated with the given minimization problem has a unique positive 

definite solution. Comparing (2.68) with (2.80), It Is clear that 

which Is the solution of (2.68) Is als? a solution of (2.80) . This completes 
the proof. 

Remark While constant matrices ^ and R can be chosen such that the 
associated regulator problem leads to a closed-loop system with degree 
of stability a. It does not seem possible for us to obtain such 
matrices without first solving a LQ regulator problem with weighting 
matrices of the form 2 . R • 

Remark It Is clear from the preceding corolla;^ that RPDS for LTI 
systems Is Just a particular class of LQ regulator with a very special 
choice of weighting matrices. This Is an Important observation for 
It Implies that all the feedback properties of LQ regulators are shared 
by RPDS. 

2.4 Kalman-Buov Filter with a Prescribed Degree of Stability (FPDS) 

In this section, we turn to the problem of designing a Kalman-Bucy 
filter with a prescribed degree of stability. By stability of a filter, 
we refer to the stability of the estimation error dynamics. Mathematl- 
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cally, this Is a problem dual to the RFDS problem considered In the 
last section. To make clear the connection between the regulator problem 
and the filter problem, we consider a linear system of the following 
type 


x(t) ■ A(t) x(t) + 

C (t) 

(2.81) 

2 (t) ■ C(t) x(t) + 

i(t) 

(2.82) 


where £.(t) and^(t) are uncorrelated zero mean white noises with 
spectral Intensity H(t) and 0(t) respectively. We wish to obtain a 
linear estimate of x(t) given 2 . (t) , ^ T < t, such that the mean square 

estimation error Is minimized. Under the assumption that , 

(1) A(t) Is continuous and bounded, C(t) , H(t) and 9^(t) are 

piecewise continuous and bounded for all t (2.83) 

(11) £(t) ^ and 0(t) > $2 i ^ where and 

are some positive real constants (2.84) 

(ill) ( C(t),A(t)) is uniformly completely recons true table and 

(A(t), H^^^(t)) is uniformly completely controllable (2.85) 

It Is well known (see Chapter 4 of [Kw 1]) that the optimal linear state 
estimate x(t) Is peclfled by 

£(t) - A(t) x(t) +Z(t) C.(t) 0(t)*-^ (2(t) -C(t) x(t)) 

( 2 . 86 ) 

where^ (t) Is the steady state solution of the Rlccatl equation. 
lit) - A(t) Z(t) +i(t) A^(t) +5(t) -i(t) cj(t) H"^(t) C.(t) i(t) 


(2.87) 
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as approaches - ® for any initial condition ^ > 0. Moreover, the 
state estimate error dynamics specified by 

e(t) - (A(t) - H(t) C. (t))e(t) (2.88) 

where ^(t} «• x(t) - x(t) is the state estimation error and where 
H(t) ■ r(t) cj (t) 0(t)”^ » the filter gain, is asymptotically stable. 

Given an Initial estimation error, the rate of decay of e(t) is 
determined by that of (|»p(t,t^) where 

^(t.t^) - (A(t) - H(t) C(t)) ijCt.t^) (2.89) 

The decay rate of| |(|>_(t,t )|1 is in turn dependent on the ncise Intensity 

•r o 

matrices ^(t) and ^(t). 

In some applications, one may require a state estimate error dynamics 

that is faster than the one specified by the noise characteristics 

given by the physical systems. This can be accomplished by making 

appropriate adjustment of ^(t) and ^(t) in the filter design equation 

(2.87) so that for a given positive constant ci, there exists some positive 

-a(t-t ) 

Y such that 1 li-(t,t^)l 1 < Y e ° for all t > t . We call a 

filter with such property a KBF with degree of stability a . In view 
of the dual relationship between the KBF and the LQ regulator problem, 
the technique for designing RFDS discussed in the previous section is 
readily applicable to the design of FPDS. This is Illustrated by the 
following theorem. 

Since the duality between the LQ regulator problem and Che KBF 


problem is well known, we will omit all Che dual proofs in this section 
and simply formulate Che dual problems and state the corresponding 
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results. Special emphasis will, however, be given to result Inter- 
pretetlons that are unique to estimation problems. 

Theorem 2.6 Consider the linear dynamical system described by the 

equation pair (2.81) and (2.82). Suppose that ^(c) and G(t) are zero 

"2ott 

mean, uncorrelated white noises of spectral Intensity ^(t) e and 

O^(t) e respectively, and that ct Is some positive constant. The 

matrices £(t) and 0(t) are £tssumed to be bounded and piecewise continuous 

with H_(t) ^ I and 0(t) ^ ^2 I for all t where 3^^ and $2 are positive 

constants. Furthermore, suppose that (£(t) , A(t)) Is uniformly completely 

1/2 

recons tructable and that (A(t), £ (t)) is uniformly completely control- 

able. Then the KBF gain ^(t) obtained from solving the linear least 
square estimate of x(t) given y.(T), - « < T< t,is 

H (t) - Z ft) C^(t) 0"^(t) (2.90) 

where Z (t) Is the steady state solution of the Riccatl equation 

a 

L(t) - L(t) (A(t) + al)^ + (A(t) + al) r (t) + 5(t) 

- L(t) C^(t) 0"^(t) C.(t) z„(t) (2.91) 

as t approaches -»for any E^(t ) > 0. 
o o — 

Moreover, the state estimate error dynamics has a degree of 
stability of at least a . 

As we specialize the results of Theorem 2.6 to the time-invariant 
case by setting A(t), B(t), £(t) and 0(t) equal to some constant matrices 
A, that satisfy the required controllability, observability and 

positive definiteness conditions, the KBF gain obtained above Is also 
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time- Invariant. A precise statement o£ the results Is stated In the 
following theorem which Is the exact dual of Theorem 2.5. 


Theorem 2. 7 Consider the linear dyncunlcal system 


x(t) - A x(t) + 5^(t) 

(2.92) 

2 .(t) - £ x(t) + Q(t) 

(2.93) 


Suppose that 6(t) and ;(t) are zero mean white noise with spectral 
Intensity matrices Ge and He respectively and a being some 
positive constant. The matrices G and H are assume to satisfy 


G > 0 and (2.94) 

and Is controllable. Furthermore, suppose that a) Is 

observable. Then the KBF gain ^ obtained from solving the linear 
least square estimate of x(t) given -» < t < t la given by 


Sa 


r T --1 


(2.95) 


where ^ is the unique positive definite solution of the algebraic 


Rlccatl equation 


(A + al) E + Z (A + al)*^ + H ^ G'^ C Z «0 (2.96) 

~ -€t —a — — — €t — •— — €t 

Moreover, the estimation error dynamics described by 

i(t) - (A - H c) e(t) (2.97) 

has a degree of stability a . 


Remark Observe that both noises specified In Theorem 2.7 have spectral 
Intensity matrices that decay exponentially with time. The resulting 
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FPDSwlll therefore give more emphasis Co the recent data In vle» of Its 
greater accuracy. This has the desirable effect of enhancing the filter's 
tolerance of modelling errors* In the classical Weighted Least Square 
approach to estimation, the same effect of giving more emphasis to recent 
data Is achieved by penalizing the recent estimation errors more than 
the old ones [Sc 1]. 

We now state the dual of Corollary 2.3. 

Corollary 2.7 If ^(t) and ^(t) In Theorem 2.7 are chosen to be 

5(t) » H + 2^ ^ and £(t) - 0 (2.98) 

respectively, where Is the unique positive definite solution of the 
algebraic Rlccatl equation (256) then the resulting KBF gain Is given 
by 

(2.99) 

and the. error dynamics of the filter 

e(t) - (A C) e(t) (2.100) 

has a degree of stability at least a. 

Remark It Is clear from Corollary 2.7 that one needs to Introduce 
a very special choice of process noise In order to design a KBF with a 
prescribed degree of stability. One common approach to speed up the 
filter error dynamics is to use an observation noise v/lth spectral 
Intensity matrix of the form p0(p > 0). The value of p Is then 
decreased toward zero until a satisfactory design Is obtained. Intuiti- 
vely, decreasing the value of p will Improve the speed of the state 
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reconstruction since less attention Is now paid to the filtering of 

observation noise. However, this procedure does not always lead to 

Improvement of error dynamics. Consider the case where £Ct) Is of the 

form c (t) ■ F^(t) with F_ having the same dimension as £ and £(t) 

being a zero mean white noise. Let the spectral Intensity of £(t) be £ 

where £ Is positive definite. The spectral Intensity of £(t) Is then 

T 

given by £ • ^ Z • Then, as p approaches 0,the P filter poles will 

A 

approach the p values where 

lfRe(v^):i.O (2.101) 

A 

'^l " 

-Re(v.) + j Im(v.) If Re(v.) > 0, 

^ 1-1,7„....P 

and v^'s are the values of s (complex) that give rise to rank deficiency 
of the matrix 



C 0 


It Is possible that some of the may be located very close to 
JO)-axls. The filter poles driven toward such v^'s as p approaches 0 
will result In poor error decay for certain Initial values of the filter 
error. 

There are however situations where a filter designed with the 

latter procedure Is preferable to one designed using the method suggested 

In Theorem 2.3. If the observation noise Intensity P 6 is small compared 

- 1/2 

with the process noise intensity £ , and If (£, A, £ ) constitutes a 

minimum phase system. It can be shown that 


J 
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11m E { (x(t) - x(t)^ (£(t) - x(,t)) > » 0 (2.102) 

p+0 ' ' 

la chls limiting casct the structure of the system to be observed is 
so exploited that the estimate error £(t) cannot be driven Into the 
subspace spanned by the eigenvectors of those filter poles that are 
located In the neighborhood of the j(o-axls. Employing a FFDS with an 
arbitrarily large degree of stability In such (.Ituationa may In fact 
give rise to a large error covariance. 

2.5 Concluding Remarks 

We have studied In this chapter the extension of the RPDS method 
to the linear time-varying systesis. Given an appropriate definition of 
'degree of stability' for such systems, the exponential weighting 
technique developed by Anderson and Moore for solving the time-invariant 
RPDS problem Is readily applicable to the time-varying case. 

We also obtain a characterization of the degree of stability In 
terms of Lyapunov function. This Is a natural generalization of the 
known results In Laypunov stability theory and Is useful for establishing 
the degree of stability of feedback laws which are designed using RFDS 
technique . 

The formulation and solution of the time-varying FPDS problem 
follows dually from that of RPDS. For the time-invariant case, the 
FPDS solution admits Interpretation of Interest to design of fast 
response filters. It Is not possible In general to speed up the filter 
dynamics by simply scaling up the process noise. The solutions to the 

FPDS problem specify a class of spectral Intensity of the process noise 
that Is useful for such purpose. 
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CHAPTER III 

EIGEWSTRUCTURE PROPERTIES OF THE TIME- INVARIANT 
RECllLATORS WITH A PRESCRIBED DEGREE OF STABILITY 


3.1 Introduction 

This chapter explores several important elgenstructure properties 
of time- invariant RPDS and their potential applications. 

Some equations are obtained for the derivatives of the RPDS closed- 
loop eigenvalues with respect to the stability factor c. These equations 
are useful for the purpose of recomputing the closed- loop eigenvalues 
given a small change of . 

Equations that describe the asymptotic behavior of RPDS root-loci 
are also derived. Speciflcallyt we consider a cost functional of the 
form 


lim 


I 


^2atj^T^^) 2 x(t) + vJ(t) pR u(t)]dt; 


p > 0 


( 3 . 1 ) 


and examine the various branches of the loci traced by the closed-loop 
eigenvalues as p varies. 

To close this chapter, we employ the asymptotic property of RPDS 
poles as p approaches infinity to obtain a novel solution of the 
Regulator with Prescribed Damping Ratio (RPDR) problem. 


3.2 Properties of the Solution of the RPDS Algebraic Riccati Equation 


In this sectioo, we state two lemmas that will be used extensively 
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in the remainder of this chapter and in the next one. Both lemmaa addreaa 
the behavior of the matrix ^ , which ia the unique poaitive definite 
solution of the RPDS algebraic Riccati equation (2.68) in Theorem 2.5, 
with respect to the increment of the stability factor a. 

Lemma 3.1 Let K and K be the unique positive definite solutions 

-til -Di2 

of the algebraic Riccati equatlona 

1^^(A + a^I) + (A + + 2 - 1 r"^ - 0 (3.2) 


I^(A + a2i) + (A + U2l)^ + a - i -C<2 * - 

respectively. 

The following conditions 

(i) 0 and R > 0 

(il) (A, h) controllable 
1/2 

(iii) (2 • A) observable 


(3.3) 


(3.4) 

(3.5) 

(3.6) 


that are sufficient to guarantee the existence of Kg and Kg are assumed 

“ 1 2 

to be satisfied. Then one has ^ ^ if ou > ou« 

12 ^ ^ 


Proof: Let AK > K - K . We need to show that A K is positive definite 

- -c«i -o<2 

if a, > a.. Substitute K ■ K 4>AK into (3.2) to obtain 

1 2 — Olj^ —0^ ~ 

(K +AK) (A + a,I) +(A clI)^ (1^ +AK) + £ 

- cl 2 - - -“2 

- (K +aK) B R*^ B^(K + AJO - 0 
”®2 


(3.7) 
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Rttarranglng Che ceciu In (3.7), we obcein 


“ 1 ST^ + (“l - a2^*iy + (A - B if ^ ^ (ttj^ - ttjj) I) AK 


- MU B r“^ B^ AK + 2(a, - a-)K - 0 
- - - - 1 2-02 - 


(3.8) 


Observe thee (3.8) is an algebraic RiccaCl equation in AK. To find 
out if there exists a unique positive definite solution AK for (3.8), we 

need to check Che controllability of (A + (Oj^ - oi2)i " i ^ ^ 

1/2 -1 T 

and the observability of ((2(02 - ^2^— a ^ . A ” 1 ^ 1. " 0i2)P 

The controllability of (A - B R ^ B^ K + (a. - o«)i, B) is implied 

Ot j X ^ 

by the controllability of (A + (o^ - 02) I.* 8) , since controllability is not 

affected under state feedback. The controllability of the latter in cum 

follows from chat of (A, B) as shown in the proof of Theorm 2.3. Since 

Che controllability of (A, is assumed in the original problem state- 

-1 T 

ment, this proves the controllability of (A + (o^ “ ^2)^ * 8, A A* *!.)• 


The 


observability of ([2(a, - a,)K ISa ^“1 ” “‘2^— ^ 

1 Z -^2 2 


,-l ,T 


follows from the posltiva definiteness of K 

-02 

In order to prove the positive definiteness of AK> we need to 
demonstrate the stability of the matrix (see Theorem 3.7 (Kw 1]) 

A i A - B R"^ B^ K + (a, - a,)I - 8 r"^ B^ AK 

~ — -02 1 2 — “ 

By definition of AA> rewrite A as 


(3.9) 


^ - A " 8 r’^ B^ K + (a, - a,)I 

— z — ■ 


(3.10) 
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-1 T 

Since A-BR BK 4* a. I is stable by definition and > 0 > the 
stability of A follows iomedlately and this completes the proof. 


Remark ; Lemma 3.1 establishes the fact that ^ is an increasing 

function of a . We shall make extensive use of this result in our study 
of the REDS robustness properties in Chapter 4. 


Lemma 3.2 Let ^ be the unique positive definite solution of the 
algebraic Riccati equation 

K (A + al) + (A + al)”^ K - K B r"^ B^ K + ^ - 0 (3.11) 

~ct ~ — — — —a 

The following conditions 


(i) 2 > 0 and R > 0 

(ii) (A, controllable 

(iii) A) observable 


(3.12) 

(3.13) 

(3.14) 


that are sufficient to guarantee the existence of IC are assumed to be 

3K 

satisfied. Then the derivative matrix positive definite for all 

3a 

values of o. . 


Proof; Differentiating the left side of (3.11) with respect to a, we 
obtain 


(A + al) + ^ + al)^ + 2 K 



■1 -r 3K 
br-^b ) ^ 


+ £ 


0 . 


(3.15) 


After rearranging terms, we obtain the following Lyapunov equation 
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3 K 

gP (A+ai- B + (A+al - B -^ + (£ + 2 I^) • 0 

(3 .16) 

-1 T 

where ^ * R B ^ is Che state feedback gain. 

— ct — — —a ^ 

3K 

To show the existence of a positive definite solution to (3.16), 

we need to check the stability of (A +al. - B G^ ) (l.e. that all eigenvalues 
of (A + ai - B^ G^) have negative real part) and the positive definiteness 
of (2+20. 


The stability of (A + a jC - ^ G^) is a direct consequence of conditions 
(1), (11) and (ill) (see Theorem 3.7 [Kw IJ). The positive definiteness 


simply follows from that of Since we have thus far made 

3K 

no reference to the sign of a, the matrix derivative is therefore 


of (2 + 2 


positive definite for all values of a and this completes the proof. 


3.3 Eigenvalue Sensitivity with Respect to the Stability Factor a. 

In this section, we present two methods for finding , the 

derivative of the RPDS poles with respect to the stability factor. 


Le"nm 3.3 For any value of a and for any distinct eigenvalue X^ of 


,we have 


9X 

IT 


H „ „-l -T 
” 2 .^ 5 . — 5 . 

ij 


3K 

3a 


(3.17) 


where 3^^ and ^ are the right and left eigenvectors of A - 2 B 1. 
which are defined in the usual way by 
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(A - B R"^ K - X.I) X, « 0 
-a i- -i - 


(3.18) 


z." (A - « »’■ K 

— — — ~a 

respectively. 


x^i) -Ji 


(3.19) 


Proof ; First, we differentiate (3.18) with respect to a to obtain 

9x 

(A - B r’^ K - X. I)x, + (A - B r"^ B^ K - X. I) “T— - 0 

(3.20) 

Multiplying both sides of the above equation by the left eigenvector 

U 

2 ^ cancels the second term on the left hand side and one gets 

(A-BR"^b'^I^ - Xj^I) - 0 (3.21) 

By rearranging the terms, we obtain 




t (A-BR-^f 



(3.22) 

(3.23) 


The existence of for an arbitrary real constant a follows from 

3a 

Lemma 3.2, and this completes the proof. 


3X. 

Remark ; Lemma 3.3 is adapted from the standard result for finding 


where X . is a closed-loop eigenvalue of the system matrix (A - ~ ^ K C) 
^ 3K - ? 


(see [Th 1] Chapter 3). The derivative 


3a 


can be computed from a 
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Lyapunov equation of the form (3.16). An alternate way to compute 
Is via the Hamiltonian system matrix we define now. 



Definition 3.1 Consider a LQ regulator problem with cost functional 




(t)^ 2 xCt) + (t) R u(t)]dt 


(3.24) 


and dynamic constraint 


k(t) « A x(t) + ^ u(t) 


(3.25) 


Suppose that all the conditions required to guarantee the existence 
of a stable state feedback law minimizing J are satisfied, then the 
Hamiltonian system associated with the LQ regulator problem is given by 


where 


£(t) » ^ _z(t) 







(3.26) 

(3.27) 


z(t) 


x(t) 

c(t) 


(3.28) 


The matrix ^ Is of Interest because Its elgenstructure describes 

the solution of the LQ regulator problem defined by (3.24) and (3.25) 

(see Chapter 3 of [Kw Ij). The eigenvalues of ^ are symmetric about the 

j<i)-axls. Moreover, the eigenvalues In the left-half complex plane are 

exactly those closed-loop eigenvalues of the LQ regulator. If 
T T T 

[Xj^(t) . S^(t) ] Is a right elvenvector corresponding to a left-half 
plane eigenvalue of then x. (t) Is also a closed- loop right eigenvector 
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o£ the LQ regulator and 

i^(t)-Kx^(t) (3.29) 

where K Is the positive definite solution of the algebraic Rlccatl 
equation 

KA + A^K-KBr"^B^K + ^ - 0 (3.30) 

3X. 

In order to conpute r— for the RPDS problem solved In Theorem 2.5* we 

9ct 

consider the following Hamiltonian system 


£(t) - Z z(t) 


where 


Z 

"a 


-1 T 
- B R 


1 A + a£ 

|-a - (A + al) 

This system Is obtained by considering the LQ regulator problem with 
cost functional 


(3.31) 

(3.32) 


J “ lim [x 

‘rJ. ■ 


^(t) 2. x(t) +vj(t) R u(t)]dt 


(3.33) 


and dynamic constraint 

x(t) ■ (A + al) x(t) + B u(t) 


(3.34) 


It follows from the proof of Theorem 2.3, and the property of Hamiltonian 

Systems, that the eigenvalues of the matrix (^ - al) positioned to 

the left of o m -a are the closed loop eigenvalues of the RPDS obtained 

3X^ 3X^ 

in Theorem 2.3. We can therefore obtain by computing for 

3X. 

Xf ■ X^. To determine , simply use the following lemma which follows 

easily from the result of Lemma 3.3. The proof of which Is omitted. 



values of a > 0 and any distinct eigenvalue 

of Z_ Qj, we have 



»“i +v“; 


where 


and 



[Wi“ Z^+al] - 0 


(3.35) 

(3.36) 

(3.37) 


Remark If the right eigenvectors [x? ?. and the left eigenvectors 

-a T. 

[v ^ V ^1^ of Z - al available as is the case when solving the algebraic 

^ "H -a ”* 

Riccati equation using diagonalization method (see Chapter 3 of [i^w 1]), the 

3Xi . 

computational effort required to obtain using (3.35) is negligible. 

This is a definite advantage over the method introduced in Leamma 2.3 
which involves solving a nth order Lyapunov equation. 


3.4 Asymptotic Behavior of RPDS Root-Loci 

In this section, we examine the behavior of RPDS root-loci as the 
weight on the control vector goes to zero(or infinity). More precisely, 
the following problem is considered. 

Given a LQ regulator problem with cost functional 


J « lim I e 

t,-« J n 


[x(t)^5,x(t) +u^(t) pRu(t)]dt; 


P > 0 


(3.38) 


and dynamic constraint 
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Our objective is to study the root-loci traced by the closed-loop 
regulator poles as p varies. We first consider the case of a single- 
input KPDS In section 3.4.1. Then we extend these results to the multiple- 
input case In section 3.4.2. 


3.4.1 The Single- Input Case 

Assume a state weighting matrix of the form 


1 


T 


» c c 


(3.40) 


T 

where {c , A) la observable. If a ■ 0, this reduces to a conventional 
optimal root-locus problem. Either the root-square locus method of 
Chang [Chg 1] or the root-locus method of Kwakemaak and Slvan [Kw 1] may 
be used for this purpose. 

In order to adapt these techniques to cases where a > 0, we will 
used the transformed LQ problem employed In the proof of Theorem 2.3. 

The cost functional and the dynamic constraints are given by 


and 


11m 



[x(t) £x(t) + u(t) R u(t)]dt 


(3.41) 


x(t) - (A + al) x(t) + B u(t) (3.42) 

respectively, where x(t) ■ x(t) and u(t) ■ u(t) . This transformed 
problem Is In a form where conventional optimal root-loci techniques 
are applicable. Moreover, recall from our discussion of Theorem 2.4 
that the closed-loop regulator poles of the original problem can be 
easily obtained by subtracting a from those of the transformed problem. 

We can therefore plot the optimal root-loci of RPDS as a functional of 


p In two steps. 
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Step 1 ; Plot the optimal root-locus of the txansfomed LQ regulator 
problem using conventional techniques. 

Step 2; Shift the entire root loci to the left by -a . 

The following theorem summarizes the asymptotic properties of SFDS 
root-loci In the single-input case. 


Theorem 3.1 Consider the LQ regulator problem with cost functional (3.41) 

and dynamic constraint (3.42) where u(t) Is a scalar. Suppose that 

T T 

(A, y is controllable and that £ has the form £ ■ £ £ where (£ , A) 

T —1 

is observable. Let the transfer function £ (s^ - A) ^ be 


c'^ (si - A)“^ b - k u (s - V J ; k ^ 0 (3.43) 

1-1 


7 (s - u ) 
1-1 ^ 


where v'^s and irj^s are the zeroes and poles of the transfer function 
T -1 

£ (s ^ - A) b respectively. 

Then the following properties hold 


(a) If p approaches 0, p of the n closed-loop poles of the RFDS 
asymptotically approach the number 1 - lt2....*p where 


if Re(Vj^) ^ -a 


-v^ -2a if Re(v^) > -a 


(3.45) 


(b) As p approaches O.the remaining n-p optimal closed- loop poles 
asymptotically approach straight lines which Intersect at 
(-a* 0) and make angles with the negative real axis of magnitude 
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+ Jl 

- n-p 


A ■ 0,1,2,... if n-p Is odd 


(A+ |)ir 
n-p 


A ■ 0,1,2,... “ 1 if n“P Is even 


(3.46) 


These faraway closed-loop poles of RFDS are asymptotically at a distance 
2 ^ 

04q from (-a, 0). Moreover, k is independent of a. 


(c) As p approaches * , then n closed- loop poles approach the 


numbers 


ir^, 1 * 1,2,3,... n where 


TTi - 


-T. - 2a 


if Re(TT^) ^-a 


if RgCTT^) > -a 


(3.47) 


Proof ; When a« 0, Theorem 3.1 reduces uo Theorem 3.11 of [Kw 1]. The 

cases where a > 0 trivially follow from the construction procedure of 

the RFDS root-loci described above. 

« 

The following example gives several root-loci plots obtained using 
the previous construction procedure. 


Example 3.1 Consider the RFDS problem with cost functional given by 


ft 

ft 

ft 

1 [2 1] x(t)] + P u^(t)l 


U -1 J 


dt 


subject to the dynamic constraint 


(3.48) 



The loci 


T -1 s + 2 

The transfer function £ (si - ^ b Is given by ( a4.3) ' ( ' s»4 7' ‘ 

of the closed-loop poles corresponding to RPDS design with values of a 

equal to 0,1, 2. 5 and 3.0 are plotted in Fig. 3.1 a,b,c and d respectively. 

As expected, the RPDS root-loci lie to the left of the O ■ -a line 

in all 4 cases. When the value of a is less than 2.0, one of the closed- 

loop pole is drawn to -2.0 with the remaining pole going off to Infinity 

along the negative real axis. (Thus forming a first order Butterworth 

pattern). In cases where a is larger than 2.0, the finite poles approach 

-2.0a + 2.0 asymptotically. (Property (a) of Theorem 3.1). 

It is interesting to note that -2.0a 4* 2.0 equals 3.0 when a ■ 2. 5. As 

a result, the open-loop pole located at -3.0 remains fixed asp varies (Fig. 3.1c) 

Property (a) of Theorem 3.1 has important design implications. 

Let us consider a LQ state feedback design configuration of the form 

T 

depicted In Fig. 3.2 where the output matrix £ in the figure is identical 

to the square root of the state weighting matrix ^ in (3.48). If • 0.0 
T -1 

and If £ (s £ - A) ^ is minimum phase, thun by Theorem 3.1 the p 

finite closed-loop poles will asymptotically approach the p minimum phase 
T -1 

zeroes of £ (s £ - A) ^ as p approaches 0. These finite modes are 
therefore hidden from the output y(t) . The remaining n-p visible modes 

f 

form a Butterworth pattern which is known to give step response of small 
overshoot. In cases where a > 0, the root- loci are constrained to lie 
to the left of c ■ -ot . By property (a) of Theorem 3.1, the p finite modes 
will be hidden from the output only if the p minimum phase zeroes of 
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FIG. 3.1 RooC-Locl of a Single-Input RPDS with 


Different Values of a 
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T -1 

£. (s 1. * A) b are located to the left of or ■ 

la light of the above observation, one needs to exercise caution 
In using RPDS for state feedback design. If the finite optimal poles 
can be made 'reasonably* unobservable at an acceptable bandwidth, then 
the use of RFDS may not be desirable since It nay Introduce some slow 
modes In the observed output which would otherwise be absent. However, 
there may be situations where excessive bandwidth Is required before some 
of the finite modes are made unobservable In the output. It Is advantageous 
In such cases to consider the use of RFDS so that the finite modes that 
appear In the output are guaranteed to decay at a certain rate. 

The following example demonstrates the relationship between the 
stability factor a and the observed output y(t) as the weight on the 
control becomes vanishingly small. 

Example 3.2 Consider a RPDS state feedback arrangement of the form 
shown In Fig. 3.2. The underlying system dynamics Is specified by the 
differential equation (3.49) and the state variables x(t) are related 
by the observed output y(t) by y(t) ■ [2,1] x(t). The RFDS state feedback 
gain Is obtained by solving the RFDS problem described In Example 3.1 
with a chosen to be 1.0 and 3.0 and the scaling factor p In (3.48) set 
to 10000. 

T 

An Initial state perturbation equal to [10,10] Is applied. The 
resulting time simulation for the observed output y(t) and the control 
u(t) are plotted for the closed- loop RFDS systems with a « 1.0 (Fig. 3. 3 (a)) 
and a - 3.0 (Fig. 3.3(b)). 

Despite the presence of a slow mode near -2 (Recall from Example 3.1 




y(t), u(t) 
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FIG. 3.3 (a) Tine Rasponse of a RPDS where Che Slow Mode 
is Cancelled by Che Syscem Zero 






-56- 


and Fig. 3.1} > the RPDS design with a« 1.0 Is noted to decay rapidly. 

This Is a result of the fact that the closed-loop eigenvalue near -2 Is 
hidden from the output y(t) as a result of pole-zero cancellation. However, 
the same type of cancellation does not occur for the design with ot equal 
to 3.0. In this case, all the closed-loop poles are constrained to lie 
to the left of a > -3.0 line, thus preventing them from being hidden by 
the system zero at -2.0. This Is evident from the slower decay observed 
In this case. 

The presence of the slow mode Is also responsible for the relatively 
long settling time of the feedback signal u(t) observed In both Figs. 3. 3(a) 
and (b) . 


3.4.2 The Multiple-Input Case 

The procedure described for constructing the RPDS root-loci In the 
single- Input case Is equally applicable to the multiple- Input case. The 
asymptotic behavior of the modes that stay finite as P approaches 0 are 
the same for both cases. However, the far-off closed-loop poles In the 
multiple- Input case generally do not form a single Butterworth pattern, 
but group Into several Butterworth configurations of different orders and 
different asymptotic rail. The exact detail of such patterns Is not 
considered here. Interested readers may consult [St 1] and [Kw .2] for a 
thorough treatment of this subject. The corresponding multiple-input RPDS 
results are summarized In Theorem 3.2 and demonstrated In Example 3.3. 

Theorem 3.2 Consider the LQ regulator problem with cost functional (3.38) 
and dynamical constraint (3.39). Suppose that (A, ^) Is controllable and 


T 

that ^ is of the form £ ■ C C where £ has the same rank as and 
T 

where (£ , A) Is observable. Let H(s) be the transfer function matrix 

H(s) • £^(5 I - A)“^ B (3.50) 

Suppose that <|)(s) <■ det(s£- A) and write 


det (H(s)) 


<Ka) 


-k 


where k 0 and where the are the 

facts hold. 


P 

IT (s - V ) 

i-1 

n 

7T (s - TT ) 
i»l ^ 


poles of A. 


(3.51) 


Then the following 


(a) As p approaches 0, p of the optimal closed- loop poles approach 
the values v^, i > 1,2,... p where 


if Re(v^) £ -a 


- 2a if Re(v^) > -a 


(3.52) 


The remaining closed-loop poles go to. infinity and group into 
. several Butterworth configurations of different orders and different 
radii. 

(b) As p approaches , the n closed-.! cop RPOS poles approach the 
numbers tt^, i « l,2,...n where 

if Re(ir^) £ -a (3.53) 



2a 


if Re(TT^) > -a 
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ZI2®£! The results follow trivially from Theorem 3.12 of [Kw IJ and the 
root-loci construction procedure described above. The detail of the proof is 
omitted. 


Consider the following multiple-input RFDS problem 


min J 
u(t) 


lim I 

ti-^ jg 


e^°^^ [x^ £ + o^(t) u.(t)Jdt 


where 


0 -5 2 -2 1 T 


8 -14 0 


and subject to the dynamic constraint 


7 -6 


0 1 


2 x(t) +10 u,(t) 

8 2 0 


-2 0 


It can be readily shown that the system zeroes for [A, £ ] are 

located at -1.0 + j 1.0. The loci of the closed- loop poles corresponding 
to RPDS designs with values of oi equal to 0.0, 0.5 and 2.0 are plotted 
in Figure 3.4 a,b,c respectively. In all three cases, the various branches 
of the loci stay In the half plane to the left of the -a line. More- 
over, two of the four branches asymptotically approach infinity along 
the negative real axis, forming two first-order Butterworth patterns. 

The asymptotic behavior of the two remaining branches depends on the 
value of cx. For ex less than +1.0, these branches eventually arrive at 


- -ya 
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fig. 3.^ Root-Loci of a Multiple- Input RPDS with Different Values of a 
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a s 2.0 


40 -30 - 20 | 


FIG. 3,4 Root-Loci of a Multiple-] 

Values of a 
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the system zeroes -1.0 + J 1.0. If a Is larger chan 1.0 » these branches 
approach Instead -2a '^1 + j 1.0. 


3.4 Regulators with Prescribed Damping Ratio 

A class of LQ regulators are chose with prescribed damping ratio 

(RPDR) . A LQ regulator Is said to have a damping ratio cos 6 if It Is 

a. 

stable and each of Its closed-loop poles a + jw satisfies 1 l5_tan 9 

^ " 1 

for some given 9. Dlagramatlcally (Blgure 1.2), all the closed- loop 
poles are found In a cone centered at Che origin of the complex plane with 
edges extending along an angle of size 26 which Is symmetric about the 
negative real axis. For stability, the value of 6 has to lie between 
0 and ^ . The procedure described In the following theorem solves the 
RFDR problem as a special case of Che RFDS problem. 


Theorem 3.3 Consider the following LQ regulator problem 

min J ■ lim f e^^ £ i(c) +u^(t) pRu,(t)]dt 

u(t) 

± ^ ^ . n > n 


(3.57) 


subject to Che dynamic constraint 

x(l) • a x(t) + B u(t) (3.58) 

Suppose that all Che requirements on A, B, £ and R that guarantee 
the existence of a stable minimizing control law are satisfied and that 
Che scalar P Is an arbitrary positive constant. The factor a is 
selected In the following manner. For each open- loop pole 0^ + jw^ of A, 
we associate a positive scalar a^ which Is defined by 
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a 


i 


0 


if 


T r < tan 0 


-(tan 



otherwise . 


(3.59) 


a Is then chosen to be the largest of the a^'s and used in (3.57). 

•1 T 

Let u(t) ■ -R B K x(t) be the feedback control law obtained 

~ ~ “Cl — 

from the above minimization problem as p approaches , where 
K ■ lim K (p) and K (p) is the unique positive definite solution 

pfoo 

of the algebraic Rlccati equation 


^(p) (A + a I) + (A + a I)^ ^(p) - ^(p) B i r'^ I^(p) + £ - 0 

(3.60) 

-1 T 

The LQ regulator resulting from the application of _u(t) ■ -R 2 . ^ 
to (3.58) has a damping ratio equal to cos 3 . 


Proof; Let + j oj^’s be the poles of the closed- loop system described by 

• “T T 

x(t) ■ (A - B R ^ I^) x(t). Applying the results of the property (c) 
of Theorem 3.2 to the previous minimization problem, we obtain 


a . + jco. 


1 0^ + j 0)^ if < - a 

-Of - j (0^ - 2 a if 1 - a 


Me need to show that for each i, the inequality 



tan 6 


is 


(3.61) 


satisfied. There ar? two different situations that need to be considered 
individually. 


Case (I) 



< tan 0 


It follows from (3.61; that 


for all values of a^. This 
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In turn Implies that 



< 

‘-I 

Ol 


^1 


Case (II) 


> tan 6 


By definition of a, and |a)^| > |a^|tan 9, wa have 


^ tan0 


hi . 

tan6 

2 


Hence ^ -a and ■ -o^ - 2a . 

Applying the definition of a and we obtain 


cr « -c - 2a 
11 


i -”i - 


'i - 2 


It follows from (3.64) and (Property (c) of Theorem 3.2) that 

0). lo). I 

.< tane (: 

ki I 
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3.6 Concluding Remarks 

We have developed In this chapter several techniques for analyzing 
the elgenstructure properties of RFDS. 

Two methods for computing the sensitivity of the closed-loop RFDS 
poles with respect to the stability factor a are introduced. The one 
based on the Hamiltonian system is computationally useful for updating 
the closed-loop poles given small variation of a. 

We also present a two-step procedure for plotting the root-loci 
for single-input BPDS. This procedure also provides the necessary frame- 
work for derivation of asymptotic root-loci properties for multiple-input 
RFDS. Using the property of the RFDS root- loci as the state weightings 
become vanishingly small, we obtain a novel one-step solution to the 
Regulator with Frescribed Damping Ratio problem. 


CHAPTER IV 


ROBUSTNESS PROPERTIES OF REGULATORS WITH A 
PRESCRIBED DEGREE OF STABILITY 


4.1 Introduction 

A critical property of feedback systems is their robustness* i.e. 
their ability to maintain system performance in the face of uncertainties. 
In particular, it is important that a closed'* loop feedback system remains 
stable despite the difference between the model used for design and the 
actual plant in the absence of feedback gain recomputation. Such 
differences commonly arise as a result of unknown and/or unmodelled 
dynamics of a plant. 

So far in this thesis, we have been looking at RPDS design from a 
transient response point of view (i.e. the ability of the feedback system 
to attenuate the initial state perturbation at a prescribed rate). The 
robustness specifications, coimnonly quantified in terms of stability 
margins and noise attenuation requirements, have been absent from the 
RPDS problem formulation. The purpose of this chapter will be to 
discuss the robustness properties, particularly those related to stability 
of RPDS and its dual, the Kalman Bucy filter with a prescribed degree 
of stability (FPDS). 

Many of the results (e.g. Theorems 4.8 and Corollaries 4.2, 4.3 
and 4.4) presented in this chapter and their robustness interpretatlo>is 
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have been created by Anderson and Moore for the single-input RPDS problem. 
The objective of this chapter is to generalize these results for the 
single-input case to the multiple-input systems. 

Required background for robustness analysis of linear time-invariant 
MIMO feedback systems is briefly reviewed in Section 4.2. The use of 
minimum singular values of the return difference and Inverse return 
difference matrices as measures of MIMO feedback system's ability to 
tolerate modelling uncertainties is emphasized. 

In Section 4.3, an appropriate framework for robustness analysis 
of RPDS is introduced. It is shown that the RFOS state feedback loop 
can be redrawn as a unity negative feedback system, which is the cannonical 
feedback structure assumed in the current works on robustness theory. 

Baaed on such a framework, the properties of the return difference 
matrices and the inverse return difference matrices of RPDS and their 
corresponding robustness implications are studied in Sections 4.4 and 
4.5 respectively. RPDS is shown to possess excellent stability margins 
with respect to the stability and degree of stability properties. 

However, the ability of RPDS to tolerate modelling uncertainties only 
Improves with increasing value of a under very specific context. 

Section 4.6 discusses the issue of roll-off requirement at high frequencies. 
This problem is of importance from the robustness point of view because 
the quality of a nominal design model inevitably deteriorates at high 
frequencies as a result of unmodelled and/or unknown dynamics. ^ 
explicit relation between the cross-over frequency and the choice of a 
is derived. The robustness results obtained in Sections 4.4, 4.5 and 
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4.6 are briefly summarized in Section 4.7. 

This chapter closes with a discussion on the dual robustness 
results of FFDS in Section 4.8. 

4.2 Robustness Analysis of Linear Time- Invariant MDfO Systems 

The robustness analysis of feedback systema requires the determination 
of regions about the nominal model for which a particular system property 
is preserved. We shall focus only on the robustness with respect to 
the closed- loop stability property. 

The importance of obtaining robustly stable feedback control systems 
has long been recognized ([Bo 1], [Ho 1]}. In classical frequency domain 
techniques for single-input single-output (SISO) control system design, 
the robustness issue is naturally handled. The various graphical means 
(e.g. Bode plots, Nyquist diagrams and Nichols charts) for displaying the 
system model in terms of its frequency response allow the control syst- a 
designers to determine by inspection the minimum change in the frequency 
response of the model dynamics that leads to instability. These changes 
are commonly quantified in terms of the gain and phase margins.^ 

Extension of these SISO robustness measures to the MCiO case is by 
no means straight forward. A satisfactory notion of stability margins 
for a multi-loop feedback system must be able to characterize the 
ability of a system to tolerate gain and phase variations in all of 
its loops simultaneously. It is not until recently that an appropriate 
framework for this purpose has become available. The basic work in this 

1 


Readers are referred to section 3.2 of [Le 1] for a comprehensive 
discussion uf these measures. 
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area la due co Sefonov [Sa 2J. Hla work waa fonaulated in che cl&e 
domain and uaed acme baalc concepca of functional analyala* aa in modem 
Input-outpuC formulation of atablllty theory. The approach adopted here 
for analyzing robuatneaa propertlea of time-invariant RPDS la baaed on the 
frequency domain formulation developed by Doyle [Do 4] and Lehtomakl [La 1]. 

A review of MIMO ayatem robuatneaa reaulta that are relevant to 
the aubaequent development of chla chapter are preacnted In aubaectlona 
4.2.1 to 4.2.4. 

4.2.1 Characterization of Model Error 

We ahall baae our dlacuaalon for the remainder of thla chapter on 
the unity negative feedback ayatem depicted In Fig. 4.1. T(a) la aaaumed 
to be a atrlctly proper rational tranafer matrix (l.e. che atate apace 
realization haa no feedthrough term) t^lch repreaenta the plant and any 
compenaatlon Chat la uaed. The perturbed veralon of T(a) la denoted by 
T^(a} where ^(a) la again a proper rational tranafer matrix. Before 
Introducing che various robuatneaa testa for HIMO feedback ayatems, we 
need to define aome appropriate measurea of deviation between T(a) and 
T(a). There are many waya that one can repreaent Che model uncertainties. 

It la Important to point out chat different typea of model error repre- 
aencation will emphaalze different aspects of the differences between 
T(s) and T(a) and will thus give essentially different assessment of 
the robustness properties under certain circumstances.^ 

I ■ " ^ 

For a concrete illustration of this point) readers are referred to 

section 4.5 of this chapter. 




FIG. 4.1 Unity Negative Feedback System with 
Loop Transfer Matrix Given by T(s) 
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Following Che notation of [Le Ij , we let the matrix ^ (s) denote 
Che particular modelling error under consideration wheve the index i is 
used to distinguish between different types of error models* The two 
most Intuitive chazacterlzatlons of model error are given by 

Z^is) - i(s) -T(s) (4.1) 

which is the absolute error inX(s) and 

E 2 (s) - T*^(s) (T(s) - T(s)) (4.2) 

which is the relative error inr(s). It is convenient for the subsequent 
development to define Che multiplicative perturbation matrix L(s) by 

T(s) » T(s) L(s) (4.3) 

Note that L(s) has a nominal value of Based on the above definition, 
we can reexpress (4.2) as 

£^(3) * L(s) - 1 (4.4) 

A feedback representation of the perturbed system using X(^) L(s) 
is depicted in Fig. 4.2. The two error measures introduced above are 
obvious multivariable generalization of the error measures c(s) - t(s) 
and (c(s) - t(s))/t(s) which are commonly employed in classical stability 
analysis using the Nyquisc diagrams. 

In view of the use of the inverse transfer function g ^(s) in 
stability results employing the Inverse Nyquist diagram [Ro 1], it 
seems natural to introduce an alternative definition of the absolute and 
relative error between the ncridnal and the perturbed system in terms of 
T ^(s) and T ^(s). With this tn>e of error representation, the absolute 
error is given by 
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FIG. 4.2 Feedback System with Multiplicative 

Representation of Uncertainties in T (s) 
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E^Cs) - t"^(s) - t"^(s) (4.5) 

and the relative error by 

E^(s) - (i"^(s) - t"^(s)) T(s) (4.6) 

- l'^(s) - ^ (4.7) 

Several comments are in order here. First, the absolute type of errors 
defined above are additive in nature whereas the relative type of errors 
are multiplicative in nature. While both types of errors can be used to 
derive robustness theorems, the notions of gain and phase margin are 
associated only with the relative type of errors. Second, the magnitude 
of the absolute type of errors is affected by both the modelling un- 
certainties and the compensator gains. This in turn makes it very 
difficult to assess the Improvement of uncertainty tolerance due to 
compensator adjustments. Problems of this nature do not occur with 
relative type of errors, since the scaling effect of the compensator 
gain is naturally handled by normalization (i.e. by multiplication of 
^(s) and T ^(s)). In this thesis, we shall only work with robustness 
theorems that are based on relative type of errors. 

4.2.2 Multivariable Nyquist Theorem and Its Generalization 

All robustness results that we are going to present in the sequel 
are based on the multivariable Nyqulst theorem. This theorem can be 
stated as follows. 

Theorem 4.1 Tlte MIMO system in Fig. 4.1 is closed-loop stable if and 
only if the image of the Nyquist contour (Fig. 4.3) under the map 
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jw 



FIG. 4.3 Nyqulst Contour Enclosing all the Unstable 
Open-Loop Poles of ^(s) 
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det (I. + X(.8)) encircles the origin P times in the counterclockwise 
direction where P is the number of unstable poles of X(s) that are 
encircled by the Nyqulst contour as the radius R of the half circle 
becomes sufficiently large. 

Remark This theorem is a natural generalization of the familiar 

Nyquist criterion with det(^ + T(s)) taking the place of 1 + t(s). It 
is derived using the relationship 

det(I+T(s)) - 

where (s) and (s) are respectively the open-loop and closed-loop 
OXs CL 

characteristic polynomials of the underlying system, and the Principle 
of Argument of complex variable theory. 

Remark The Nyquist diagram of det(i+ X^)) is commonly called the 
multivariable Nyquist diagram. Despite the similarity in form between 
the above theorem and the Nyquist criterion for SISO system, it is not 
possible to infer robustness properties of a MIMO system by inspection 
of its multivariable Nyquist diagram. In other words, the distance 
between the multivariable Nyquist diagram to the origin of the complex 
plane does not constitute an appropriate measure of stability margin. 

The reason for this will become clear in the next section. 

The multivariable Nyquist theorem provides a procedure for checking 
the presence of closed-loop poles Inside the domain enclosed by the 
Nyquist contour. Generalization of such procedure to regions enclosed 
by other contours in the complex plane is straightforward. 
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This can be accomplished by plotting the Image of the contour In question 
under the nap det(I +T(s)) and counting the resulting number of en- 
circlements about the origin. If this number turns out to be Identical 
to the number of open- loop poles enclosed by the counCour of Interest, 
we can then conclude the absence of closed-loop poles Inside the region 
enclosed by such contour. 

Remark A rigorous Justification of the above discussion follows 
directly from the application of the Principle of Argument and (4.S). 

A c( ntour of significance fur the euialysls of feedback systems Is 
the a - Nyqulst contour. This contour Is obtained from the conventional 
Nyqulst contour by shifting the jw-axis to the line a ■ -a (see Fig. 4. 4) 
Applying the procedure discussed above to this contour results In a 
graphical test for checking the degree of stability. A precise 
description of this result is given In the following theorem. 

Theorem 4.2 The MIMO system in Fig. 4.1 has a degree of stability o. 

If and only if the Image of the <x - Nyqulst contour (Fig. 4.4) under the 
map det(I +T(s)) encircles the origin P times in the counterclock- 
wise direction where P is the number of poles of chat are enclosed 

by the a - Nyqulst contour as the radius R of the. half circle becomes 
sufficiently large. 

Remark The SISO version of the above theorem was employed by Anderson 
and Moore [An 1] for interpretation of RPDS robustness properties. 
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4.2.3 Robustness Theorema for MIMO Systems 

In this sectloQt we review several theorems that are key to the 
study of RPDS robustness properties. The derivation of these results 
from the basic theorems (Theorems 4.1 and 4.2) are sketched In a tutorial 
fashion. For convenience, we will assume that In all the remaining 
theorems and corollaries, the Nyqulst contour (a»Nyqulst contour) Is 
chosen with R sufficiently large so that Theorem 4.1 (Theorem 4.2) 
may be applied. 

Suppose that the nominal closed-loop system In Fig. 4.1 Is stable 
and that T(s) and T(s) have the same number of unstable poles. Then It 
follows from Theorem 4.1 that the perturbed system Is stable If and only 
If the image of the Nyqulst contour under det(I -f T^(s)) and det(I 
have the same number of encirclements about the origin. 

Consider the c'ase where T(s) Is a continuous deformation of Ks) . 

For a sufficiently small deformation, the number of encirclements of 
the multivariable Nyqulst diagram about the origin stays unchanged. 
Consequently, the perturbed system will remain stable. However, If 
the deformation makes the multivariable Nyqulst diagram crossing the 
origin, the number of encirclements about the origin will change as a 
result and the perturbed system becomes unstable. The point where 
crossing of the origin takes place Is characterized mathematically by 

det(I + T^(joi)) - 0 (4.9) 

for some value of O) > 0 where T (s) is the deformation ofX(s) that 

— — o 

touches the origin for the first time. This point clearly marks the 
borderline between stability and Instability. It seems natural from 
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the above discussion to characterize the stability margin of a system at 
each frequency (ti by the nearness of the nominal return difference 
matrix (I T(jai)) to singularity. 

The determinant detCl appears to be an obvious candidate 

for such measure since the determinant of a matrix is zero If and only 
If It Is singular. Moreover, the size of det (I +^(s)) at each 
frequency o) can be readily obtained by Inspection of the multivariable 
Myqulst diagram. Unfortunately, the determinant of a nonsingular matrix 
turns out to be an unreliable measure of closeness to singularity. It 
Is well known that a matrix with a reasonably large determinant can be 
made singular by a small perturbation (compared with the determinant 
of Its elements) . The standard measure for closeness to singularity of a 
matrix Is given In terms of the matrix norms. When the 2-norm Is used 
(as Is the case of this thesis), the distance from singularity for a 
given matrix A Is quantified by the minimum singular value a min (A). ^ 

The precise use of oaiin (A) In detection of closeness to singularity 
Is given In the following theorem. 

Theorem 4.3 Given a nonsingular complex matrix A, a solution to 

the following minimization problem 

1 1 1 1 1 • ^ singular 

is given by 

1 1 E 1 1 ■ amiii (A) 

and u 

E • amin (A) U ^ 

^ Readers unfamiliar with singular value may refer to [Str 1] for an 
overview of their properties 


(4.10) 

(4.11) 

(4.12) 
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where U and V are the right and left singular vectors of A corresponding 
to Onln(A) . 

In other words, If A ^ Is singular, then 

Il 5 .il (^) ^ Omln (A) (4.13) 

The following corollary to Theorem 4.3 Is Che key to the derivation of 
robustness theorems that are based on Che multivariable Nyqulst theorem. 

Corollary 4.1 Let A be a nonsingular complex i&tiTlx and 

omax (^) < omln(A) (4.14) 

then A + E Is nonsingular. 

Remark The condition (4.14) Is sufficient but not necessary, because 
only magnitude Information of A and ^ are employed. 

Before applying the previous corollary to derive Che robustness 
theorems for feedback systems, note that we can rewrite the perturbed 

return difference matrix 

I+i(s) • I+T(s) L(s) (4.15) 

In two different ways 

I +I(s) J.(s) - (l'^(s) +T(s)) L(s) 

- [(I + T(s)) (l"^(s) -p] L(s) (4.16) 

and 

I + T(s) L(s) ■ T(s) (2 ^(s) + L(s)) 

- T(s) [(I +l"^(s)) + (L(s) - T)] 


(4.17) 



-80- 


Dlrect application of Corollazy 4.1 to (4.15) reaults in the following 
condition that guarantees the nonsingularitjr of (I^ T(s) L(s}). 

anax(L ^(s) “ p ^ anin(I + 1 ( 3 )) (4.18) 

In view of the discussion at the beginning of this section, it is clear 
that any perturbation matrix L(s) satisfying condition (4.18) for all 
values of s along the i^yquist contour will not destabilize the closed- 
loop system. This result is formally stated in the following theorem. 

Theorem 4.4 If the MIMO feedback system in Fig. 4.1 is stable, then 
the perturbed system in Fig. 4.2 is stable if 

(1) T(s) and T(s) *• T(s) L(s) have the same number of unstable 


poles (4.19) 

(ii) T(s) has no pole along the ja>-axls^ (4.20) 

(ill) ^(jcu) has no eigenvalue at 0 or on the negative real 

axis for ^ 0 (4.21) 

(iv) amax(L ^(jw) - I) ■ amax(^(jw)) 

< amin (I + T(j 0 )) (4.22) 

where E^(s) - ^"^(s) - t"^(s )) T(s) (4.23) 


Another useful nonsingularity characterization of (I T (s) L(s))) 
is provided by the condition 

omax(L(s) - I) < amin(I +T ^(s)) (4.24) 

This is readily obtained by applying Corollary 4,1 to the alternative 

I 

This condition can be relaxed by Indentiog the Nyqulst contour along 
the ju-axis at locations of open-loop imaginary poles 
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expresslon of i.+ T(s) L(a) as given in (A. 17). The related robustness 
theorem is stated below. 

Theorem A. 5 If the MIMO feedback system in Fig. A.l is stable, then 
the perturbed system in Fig. A. 2 is stable if 

(i) ^(s) and 2(s) Js(s) number of unstable 

poles (A. 25) 

(ii) ^(s) has no pole on the Jci)-<xis (A. 26) 

(ill) oaax(L(Ju) “ 1.) ■ omaxCE^CjCi))) < Omin(^ + X 

for all wlO (A. 27) 

where ^ (s) - f^(s)(T(s) - 1(8)) (A. 28) 

Remark Strictly speaking, the argisnents given above are not enough 
to prove the robustness theorems in this section. A major part that 
has been omitted la the embedding argument which ensures that the 
perturbed multivariable Nyquist diagram associated with X '^l.(J<^) can 
be reached through a continuous deformation of the original Nyquist 
diagram which preserves the number of encirclement of the origin. The 
results presented in Theorems A. A and A. 5 are however correct. Readers 
interested in more technical details should consult [Le 1]. 

We shall make several brief comments on the relationship between 
Theorems A. A and A. 3. Observe from inequalities (A. 22) and (A. 27) 
that as omin(X + X (jw)) amin(X ^ l”^(J<**)) increases, bounds on the 


respective kind of model error become less stringent. For purpose of 
improving the system's ability to tolerate uncertainties, it would 
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seea desirable to make both the quantities Cmln(I T(jp}) and 
amlnd + T ^(Jcu)) as large as possible. This Is however Impossible 
because the return difference matrix (I and the Inverse return 

difference matrix (I 4> T~^(Joi)) are related by the following matrix 
Identity 

(I + (I +KJ“))"^ -i (^*29) 

which prevents us from making both amin(_I + T(ju)) and omin(I_ + T ^(jw)) 
large Independently. 

In the low frequency regions where Omln(T(ju>)) is large (l.e. large 
loop gain In all feedback loops), the quantity ominCI^ T(ju)) is also 
large. This In turn constraints omln(^ to assusw values 
close to unity. Cor.sequently, Theorem 4.4 will tend to give a less 
conservative test with respect to the multiplicative errors (l.e. errors 
of the form T(s) • T(s) L(s)} than will Theorem 4.3 at low frequencies. 

Likewise, when Omax(T(jo))) Is small (as Is the case in the high 
frequency region). Theorem 4.4 will provide a less conservative test 
with respect to the multiplicative errors. This Is a consequence of 
o»in(_I + T”^(ju)) being large and amin(i + T_(j(i))) being near to unity. 

Using det(I + T (s-a) L(s-a)) (l.e. the Image of the (»-Nyqulat 
contour under det(I ■♦•1(8)) in place of det (I ■♦• T (s) L(s)), and 
repeating the argument (based on the deformation of the multivariable 
Nyqulst diagram) employed In the derivation of Theorems 4.2 and 4.3 
leads directly to the following theorem on robustness with respect to 
degree of stability. 
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Theoree 4.6 If the MIMO feedback system In Fig, 4.1 has a degree of 
stability a» then the perturbed system In Fig. 4.2 possesses the same 
degree of stability If 

(1) TCs*a) and T(s-a) • T(s-a) L(s-oi} have the same number of 

unstable poles (4.30) 

(11) T(s-a) has no pole along the ju-axls (4.31) 

(111) LCjojHX) has no eigenvalue at 0 or on the negative 

real axis for oi ^ 0 (4.32) 

(Iv) amax(L ^(ju^) - I) ■ omax ^(jw-c) < mlnC^ + T(jO)-a)) 

~ ^ “ (4.33) 

for all 01 ^ 0 

where E^(s) - d“^(8) - r^(s)) T(s) (4.34) 


Theorem 4. 7 If the MIMO feedback system In Fig. 4.1 has a degree 

of stability a , then the perturbed system In Fig. 4.2 has the same 
degree of stability If 


(1) X(8“0i) and T (s-a) » T(s-a) L(s-ct) have the same 

number of unstable poles (4.35) 

(11) T(s-a) has no pole on the jcu-axls (4.36) 

(111) ffmax(L(jO)-a) - p « amax(E 2 (j^^-ot)) 

< omln(_I + if^(jco-a)) for alloi^ 0 (4.37) 


where E^(s) 


t“^(s) (TCs) - T(s)) 


(4.38) 


p^mark Our remark on the relative effectiveness of the robustness 
tests prescribed by Theorems 4.4 and 4.5 at different frequency regions 
also applies to Theorems 4.6 and 4.7. 
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4.2.4 Multi- loop Stability Margins 

A natural extension of the gain and phase margin measure to MIMO 
systems can be obtained by using Theorems 4.4 and 4.5, It turns out 
that these stability margins can be conveniently characterized In terms 
of the lower bound of amlnC^ -t* ^(Ju)) and amln(^ T for 

CO ^ 0. The definition of the multi- loop margins given below are due to 
Lehtomakl [Le Ij . Diagonal (l.e. noninteracting) perturbation of the 
multiplicative type Is assumed throughout (see Fig, 4.5). 

Definition 4.1 The multi- loop gain margin is the pair of real numbers Cj^ 
and C 2 defining the largest Interval (c^fC^) such that when f>^(s), 

1 • l,2,...m In Fig. 4.5 are all real, and satisfy the lnequed.ltles 

Cl < If < C 2 1 ■ l,2,...m (4.39) 

the closed- loop system remains stable. 

Definition 4.2 The multi- loop phase margin Is the pair of real numbers 
6 and -6 defining the largest interval (-6. ,8.) such that when 
l^(jco) , 1 - l,2,...m in Fig. 4.5 are of the form e where 4>^'s are 
real and satisfy the inequality 

”®1 ^ ^ ^ * l,2....m (4.40) 

the closed- loop system remains stable. 

Following the notation in [Le IJ , we denote the multi-loop gain 
margin of (4.39) by 

GM - (c^.c^) 


and the multi- loop phase margin of (4.40) by 
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Perturbed Plant 
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PM - (-0j^,0j^) 

In a similar fashion, we can define multi-loop gain and phase margins 
chat guarantee a prescribed degree of stability for the closed- loop 
system. These bounds ace useful for characterization of a system's 
ability to maintain a specified speed of response in Che face of un- 
certainties . 

Definition 4. 3 The multi- loop a - gain margin (denoted by QIq) is the 

pair of real numbers c^ and C 2 * defining the largest interval (c^t C 2 ) 
such that when 1 ■ l,2,...m in Fig. 4.5 are of the form 

where are real and satisfy the Inequality 

Cl < ^i < C 2 , i ■ l,2,...m (4.41) 

Che closed-loop system retains a degree of stability cx . 

Definition 4.4 The multi- loop a-phase margin (denoted by PMa ) is 
Che pair of real numbers -0^ and 0^, defining Che largest interval 
(-02,0|^) such that when i • l,2,...m in Fig. 4.5 are of the 

form e where <|i^'s are real and satisfy Che inequality 

- 01 < (|i^ < 01 , i - 1,2,..,.» (4.42) 

Che closed- loop system retains a degree of stability a , 

The interpretation of Che two types of gain and phase margins 
defined above require some explanation. First of all, these margins 
are valid for perturbations (either pure gain or pure phase changes) 
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applied simultaneously at each Input channel. This does not exclude 
the possibility of having pure gain changes at some channels and pure 
phase changes at others. Thus, they differ from those stability margins 
obtained by breaking one loop at a time. However, the word 'simultaneously' 
does not suggest that one can change the gains and phases of each Input 
channel simultaneously Inside the limit prescribed by (4.39) and (4.40) 
((4.41) and (4.42)j without causing Instability (the degree of stability 
to go below a). Secondly, the multi- loop margins only ''.over a limited 
class of perturbations. In particular, they are based on the assumption 
of diagonal L(s) . 

The following corollaries to Theorems 4.4 and 4.5 provide a 
characterization of multi-loop Q1 and PM in terms of the lower bound for 
omln(I + T(ju)) and amln(I_ + T ^(jw)). 

Corollary 4.1 If the nominal closed-loop system In Fig. 4.1 Is stable 
with no open-loop pole on the ju-axls and 

omLnil + T(j(a)) > 8 , 0 < 8 £ 1, for all O) ^ 0 (4.43) 

then the multi- loop gain and phase margins are characterized by 

and 

PM D (-2 sln"^ I , 2 sln'^ | ) (4.45) 

respectively. 

Remark The case for 8 ^ 1 Is not considered because It Is In- 

consistent with the strict properness assumption of T(s). 
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Corollary 4.3 If the nomloal closed-loop system In Fig. 4.1 is stable 
with no open- loop pole on the Ju-azls, and 

<Jmln(i + ^^(jw))> 8 for all w ^ 0 (4.46) 

then the oniltl-loop gain and phase margins are characterized by 

GMD (1-8, 1 + 6) (4.47) 

and 

PM D (-2 sin I . 2 sln"^ |) (4.48) 

respectively. 

Likewise, we can characterize the multi- loop margins with respect 
to degree of stability In terms of the bound for amln(^ + T (o-Jco)) and 
Cmln(^+X ^((^J^*^)) respectively. The result Is summarised In the 
following corollaries to Theorems 4.5 and 4.6. 

Corollary 4.4 If the nominal closed- loop system In Fig. 4.1 has 
a degree of stability a with* .lo open- loop pole on C7 « -a, and 

Omin(I + T(ja)-a)) > 6 , 0 < 5 <. 1 for all w ^ 0 (4.49) 

then the multi-loop gain and phase margins with respect to the degree of 
stability cx are characterized by 

and 

PM(^ D (-2 sin I , 2 sin"^ | ) (4.51) 

respectively. 

Corollary 4.5 If the nominal closed- loop system in Fig. 4.1 has a degree of 
of stability a with no open-loop pole on O ■ -a, and 
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amln(I + ^ 8 for all w ^ 0 (4.52) 

then the multi-loop gain and phase margins with respect to the degree of 
stability a Is characterized by 

GM(x C (1 - 6 . 1 + 6) (4.53) 

and 

PM^ C (-2 sln"^ I , 2 sln"^ | ) (4.54) 

respectively. 

It la Important to emphasize that the bounds on the multi-loop gain 
and phase sargln characterization provided by the above corollaries are 
In general very conservative. This Is a result of the fact that only 
magnitude but no structural information on ^ + T(jco) and ^ + ^^(Jca) 
are employed In the derivation of these bounds. These stability margins 
are commonly known as the guaranteed gain and phase margins. 

4.3 Formulation of the RPDS Robustness Problem 

Consider ' LTI system described by the dynamical equations 


^(t) 

■ A x(t) + ^ u(t) 

(4.55) 

2.(t) 

- C x(t) 

(4.56) 

u(t) 

- -Sa 

(4.57) 


be a state feedback control law for the above system. The resulting state 
feedback configuration Is shown In Fig. 4.6. 

Since the stability of a RPDS Is characterized by the eigenvalues 
of the closed- loop matrix ( A - B , It Is strictly a property of the 
state feedback loop and Is Independent of the output matrix^. To 
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FIG. 4.7 Feedback Configuration for Robustness Analysis of RPDS 
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perfora robustness analysis, we can redraw the BPOS state feedback loop 

In a unity negative feedback form (see Fig. 4.7) that Is compatible with 
the robustness theorems considered In the last section. The RPDS feed- 
back system Is therefore akin to that found In classical situation where 
unity negative feedback Is applied around a system with loop-transfer 
matrix given by Ga (si. - 

Now, It becomes clear that we can characterize the robustness 
properties of RPDS In terms of the two qt^antltles Omln(I + T^ (ig))) 
and o mln(.I (jc*i)) where we define 

A 

Ta(jco) - G^(jO)I - A) ^ B (4.58) 

4.4 Properties of RPDS Return Difference Matrices and Related Robustness 
Results 

Several well known feedback properties of LQ regulators are characte- 
rized In terms of its return difference matrices by the Kalman Frequency 
Domain Equallty([Ka IJ , [Le Ij). Since RPDS represents a special case of 
LQ regulator, it seems natural to begin our study of Its return difference 
by introducing two modified versions of Kalman Frequency Domain Equality 
that are obtained directly from the RPDS Rlccatl equation. The derivation 
of these results are given in section 4.4.1. The three subsections that 
follow discuss three important consequences of these equalities and their 
respective robustness Interpretations. 

4.4.1 Two Fundamental Equalities 


There are two versions of Kalman Frequency Domain Equality for 
RPDS. Each of them results from a different arrangement of the RPDS 
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algebralc Rlccatl equation. These results can be stated as follows. 


Theoren 4.8 Let the matrix be the unique positive definite 

solution to the algebraic Riccati equation 



(A +ai)^ K^+ K^A +ajy + a “ la® la “ ^ 

(4.59) 

with 


(i) ai 0 and R> 0 

(4.60) 


(ii) (A, B) controllable 

(4.61) 


(ill) A) observable 

(4.62) 

then 


(I + Tj,(-a-a))^ R (I +'T^j(-a+s)) - R + H(-s,s) 

(4.63) 

and 


(i + Ia(-s))^ l(i + Ia(*)) • R + H^(-8 .j) 

(4.64) 

i^ere 

H(5,s) - B^(Ci - Ol - A^)"^ ^(si - ca - A)"^ B 

(4.65) 

and 


^(5,s) - ^^(CI - J)‘^ (2a + 2) <al - A)"^ B 

(4.66) 


Proof ; Direct manipulation of (4.59) yields 

(-si - A - ol)^ K + j^(si - -2 

(4.67) 

upon prenultiplying and postmultiplying (4. 67) by I (-si. - 1,1 

and [(sj[ - ^ B ] respectively, we obtain 

R T^jCs-a) + (s-o) R + T^ (-s-a)R^ (a “«) • H(-«.a) (4.68) 

adding R to both sides of (4.68) gives (4.63). 

To derive (4.64), note that we can rewrite (4.59) as follows 
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A + (a+ 2a 1^) - ^ B R*^ - 0 (4.69) 

by regrouping the two tern involving . Repeating the sane algebraic 
manipulations as before with(4.69)we obtain (4.64) and this completes the 
proof. 

Remark As expected, both equalities (4*63) and (4.64) reduce to the 
fastlllar Kalman Frequency Equality for LQ regulators when a *0. 
Furthermore, If 

det(JojI. - A) #0 for all o) ^ 0 (4.70) 

we can rewrite equality (4.64) as 

(I +UJw))^ K (1 i j") (^*71) 

for values of s on the Jto-axis. 

Similarly, If 

det(Jw^ - A - o4) f* 0 for allco^O (4.72) 

we can rewrite equality (4.63) as 

Ci + Ia^-« + + Jw)) - R + H(-jco, Jw) (4.73) 

for values of s on the Ja>-axis. 

4.4.2 Common properties with IX’ Regulators 

It was pointed out In our remark to Corollary 2.5 that RPDS Is 
a special case of LQ regulators. As a result. It possesses all the 
robustness properties of the latter. A formal verification of this 
fact Is provided la this section. The following corollary to Theorem 4.8 
Is key to the derivation of all the subsequent results on LQ-related 
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properties of RFDS. 

Corollery 4.6 Let the Mtrlx K be the unique positive definite solution 

“'o 

of the elgebrsic Rlccatl equation (A, 59) (with ell conditions on A, 

£ end R satisfied) , and dat(Jca^ - A) # 0 for all cd >0. 

Then 

froof ; The inequality (4.74) holds if and only if Hpj(-jcii,J<i)) > 0 

for u ^0. The positive definiteness of j^(-Jcu,J(u) in turn follows from 
that of Ky and this completes Che proof. 

Remark The result described in this corollary is the well known Kalman 
Inequality for LQ regulators. In the single-input case (4.74) reduced to 
the following scalar inequality [Ka 1] 

ll + tQj(Jo»)l >1 foru^O (4.75) 

By inspection of the NyquisC diagram corresponding to (4.75) (see Fig. 4. 8) 
it is straightforward to observe that a single- input RPDS state feadback 
regulator has a guaranteed Of of (j, «) and a guaranteed FM of <-60°»60°). 

In Che multiple- input case, the inequality (4.74) does not provide 
a bound on the quantity amln(^ 4 j^(j<t))) for arbitrary choices of R. 

In cases where R is chosen to be a scaled identity Ci.e. R ■ c^ for 
some c > 0), then (4.74) reduces to 

(I + T^(jw»*^ a * Tjj(jw)) > I for all « > 0 (4.76) 

which in Cum provides a lower bound 


amln(i ♦ TaO)) > 1 


tax all (I) > 0 


(4.77) 
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for the minioua singular value of ^ 4* Tg(j(4. It follows Immediately 
from Corollaries 4.2 and 4.3 and (4.74) that RPDS designs with R chosen 
to be a scaled Identity matrix have a guaranteed multi- loop gain margin 
given by (0.5,o>) and a guaranteed phase margin given by (-60*^, 60°). 

It turns out that these guaranteed stability margins also apply 
to RPDS designs with R chosen to be diagonal. The derivation of this 
more general results Is due to Lehtomakl (see [Le 1], [Le 2]) and Is 
Included here for completeness. The proof Is motivated by the observation 
that (4.74) can be rewritten as 

^ H 

(i + Sx^jw)) ^ - for all 01^0 (4.78) 

where we define 

^(joj) - R^^^ (4.79) 

Inequality (4.78) then provides a bound 

amln(^ + ^(ju)) > 1 for all (4.80) 

on the minimum singular value of I. + Tq(J(o). If ^(Jco) Is used Instead 

stability test based on Theorem iA , It Is necessary to 
manipulate Fig. 4.7 into the equivalent form depicted In Fig. 4.9^ 

Using Theorem 4.4 and (4.80) leads directly to the following result 

Theorem 4.9 Given a RPDS system with loop-transfer matrix (s) » the 
perturbed regulator with loop-transfer matrix (a) * it^s) Is 

closed-loop stable If 

^ This ensures the equivalence (from stability point of view) of uhe feed- 
back loops In Fig. 4.7 and Fig. 4.9. 
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(1) T (s) and T (s) have the same number of unstable poles (4.81) 

•“Ot “Ot 


(1^) joi-aacis 

(til) omaxCR^^^ l" < 1 


(4.82) 

for all (0 > 0 (4.83) 


• 1/2 "“1 ** 1/2 

Observe that the condition amas(R L (ju) R *■ "i) ^1 In 


(4*83) can be written as (see ISa 3J) 


R LCjw) + L CJW) R - R 10 


for all o» > 0 (4.84) 


If R and L(jco) are diagonal, (4.84) slnrjllfles to the following conditions 


Re(l^(jco) 2 l| 


for all <u ^ 0 
and 1 * 1,2, . . .m 


(4.85) 


The corresponding guaranteed gain and phase margins properties can be 

readily obtained from (4.85) by assuming 1 . (jw) to be a real scalar and 
j<l>l 

^^(j“) ■ ® » respectively, for 1* l,2,...m. 

4.4.3 The Effect of the Stability Factor a on RPDS Return Difference 
Matrices 

The effect of the stability factor a on the matrix 
(_I +T(x(-s))^ R(I +Xa(s)) is provided by the following corollary to 
Theorem 4.4. 


Corollary 4.7 Let the matrices and ^ be the unique positive 

** 1 2 

definite solutions of the RPDS algebraic Rlccatl equations 

(A + a.I)’^ K + K (A + 0 .. 1 ) + a - ^ 1 ^ ^ 

- 1 - 1 


**It can be trivially shown that condition (ill) in (4.79) automatically 
guarantees condition (111) of Theorem 4.4. 
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aud 


(A+a2i)^ 


K 

-TX2 


+ K (A + a-»I) + Q - K B r’^ 

1X2 ~ 2- ^ -Oj 



respectively with 

. (i) a> 0, R > 0 

(11) B) controllable 

(111) (2^^^, A) observable 

(Iv) det(jul - A) 0 for all w > 0 

then 

(I + ^ (ju)® R (I (jw)) 

> (i + T(j (jco))® R (i + (jw)) for all 0 ) ^ 0 


If > 0 I 2 — ® 


Proof ; Direct application of Theorem 4.4 to the algebraic Rlccatl 
(4.86) and (4.87) yields 

R (I +T^^(j«)) - |^^(>j<0,jcu) +R 

and 

((i + X. (jw))® R i + S. 

“®2 ” 2 2 

respectively. It follows from Lemma 3.1 and condition (Iv) In the 
theorem statement that 

Si (-Jw.Jw) > (-jw.jc^) for all o» > ( 

1 2 

Adding R to both sides of (4.95) yields 

R + H (-joj.ju) > R + H (-Joj.jo)) 

-txi -a2 


(4.87) 

(4.88) 

(4.89) 

(4.90) 

(4.91) 

(4.92) 

equations 

(4.93) 

(4.94) 

I (4.95) 


for all 0 ) > 0 


(4.96) 
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and the theorem Is proved. 

For single-input BFDS (4.92) reduces to the following scalar In- 
equality 

U + t^ (jcu)l > [l + t^ (jaj)i for all <o ^ 0 (4.97) 

Geometrically, (4.97) means that for each value of frequency (o , the point 

t (Joi) In the complex plane Is always farther away from the critical 
1 

point (-1,0) than the point t (Joi) (see Fig. 4.10)). 

U2 

In the multiple- Input case, the Inequality (4.92) does not In general 
+ (Jw))® (I+X,(jo))) > (I+T (jo)))“ (I+T (jco)) for 

-Oi -Oj^ -^2 -02 

> a^. However, If R Is chosen to be ^ scaled Identity matrix, (4.92) 
then simplifies to 

(I + X, + (»)® (I+X, OJ) (^‘98) 

■^1 "®1 "“2 "“2 

for all <») ^ 0 

The following Inequality on the minimum singular value of RFOS return 
difference can be readily derived from (4.98) by using ':he properties 
of singular values 

cmln(i + ^ (jw)) > amln(I + T (joj)), for all u ^ 0 (4.99) 

1 ^2 

Inequalities (4.98) and (4.99) allow us to assess the effect of o on the 
robustness properties of RPDS design with R chosen to be a scaled Identity 
matrix. It follows directly from condition (4.21) In Theorem 4.4 that 
the tolerance of modelling error of the form E^(s) ■ (T ^(j(s) “ T \(s)) T^(s) 
Improves with Increasing values the stability factor a. Unfortunately, the 
use of the matrix function T (s) is not cornnon, thus making the above 
conclusion somewhat difficult to Interprets. We shall explore further 
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Diagrams for RPDS Loop Transfer 
ns t (s) and t (s) with a. > a. 
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coasequences of this result by restricting ourselves to perturbation 
matrices L(s) that are diagonal* Then condition (.4*22) in Theorem 4,4 
simplifies to 

- 1 |<6(J0)) - omin(I + T^(jo»)) 

for all 0 ) ^ 0 (4.100) 

1 * 1,2*. ...m 

The allovable values for as defined by (4.100) are depleted In 

Fig. 4.11. It Is clear from the diagram that the area of allowable 

(shaded) region for l^(Jco) Increases with B (jeo) which Is Itself and 

Increasing function of a. This represents a very special type of 

robustness Improvement for RFDS made possible by Increasing the value of a. 

Attempts the find a characterization of robustness variations with a for 

error models defined In terms of T (s) have not been successful. This Is so 

“Cl 

because robustness theorems relating amln(^ ^ (jeo)) and error models 
defined In terms of I( 2 (s) are only available for the absolute type of 
error representation. We have already pointed out In section 4.2 the 
difficulty In accounting for the effect of compensator with this type of 
error model. Indeed, an Increase In the value of a will Increase both 
the error and the quantity <7mln(^<f ]^(ju)) , thus making It Impossible to 
ascertain the effect of a on robustness. 

We shall conclude this section by discussing the effect of a on 
the guaranteed and actual stability margins of RFDS. Since X^(s) Is a 
strictly proper, rational matrix function. It follows that 

lim amin(_I +T,^(jw)) for all values of a > 0 (4.101) 

This In turn Implies that we cannot find a lower bound larger than 1 for 
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FIG. 4.11 Set of Allowable Values (shaded) for &^(s) Corresponding 

to aoin(I, + (jco)) ■ (4.12a) and aaln(Jt + (Jw)) • 6 

(4.12b) 
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<mLn(l+ ^(Jco)) that holds for all values of <u. As a result, the 
guaranteed gain and phase oarglns for RPDS with R chosen to be a scaled 
Identity matrix Is unaffected by the choice of a. It can be shown that 
this result also holds for RPDS design when R Is chosen to be diagonal. 

The effect of a on the actual stability margins of ?^S turns out 
to be highly system dependent. As expected, the upward gain margin Is 
always « (since actual gain margin can be no worse than guaranteed gain 
margin). The downward gale margin and phase margin may however either 
Improve or deteriorate with Increasing a. The following simple example 
Illustrates this point. 


Example 4.1 Consider the following RPDS problem 


min 11m 
u(t) 



e^*" (x^(t) + u^(t))dt 


x(t) • ax(t) + u(t) 


(4.102) 

(4.103) 


If a ■ -1, the RPDS feedback gain Is given by 

^ - (a - 1) +((a - 1)^ + 1]^^^ (4.104) 

and Che corresponding gain and phase margin are given by 

GM • ■(- gx"^,*) (4.105) 

and , /» 

-1 9 1/9 -1 2 

PM - (-TT + tan ^(ga - 1) ^ - tan ^(g® - 1^ ^ 

(4.106) 

respectively. It follows from the fact that g^j Is an increasing function 
of a that both the downward gain margin -g^ ^ and the phase margin given 



105- 


by (4.106) deteriorate with Increasing a (see Fig. 4.12a). If we have 
a • 1 Instead, (i.e. the open- loop system is unstable) then the gain and 
phase margin of the resulting design become 

- (ga^.“) (4.107) 

and 

PM - (-tan"^(ga^ - 1)^^^, tan"^(ga^ - 1)^^^) (4.108) 

respectively, where g^ is now given by 

g^ - (a + 1 ) + [(a + 1 )^ + l]^/2 (4.109) 

Note that in this case, both the downward gain margin and the phase margin 
Increase monotonlcally with a (see Fig. 4.12b). 

4.4.4 Properties of the RPDS Return Difference Matrices on the 
q-Nyquist Contour 

Thus far, we have not been using the version of the BPDS Frequency 
Domain Equality given by (4, 64). a useful matrix inequality similar to 
(4. 74) can be derived from this equation. 

Corollary 4.8 Let the matrix ^ be the unique positive definite 

solution of the algebraic Riccatl equation (4.59) with the respective 
conditions on A, ^ , 2. &nd R satisfied and det(Jo) ^ ~ A -q^) ^ 0 for 
all (0 >0, then 

(I ■‘"Iq^"® ■** w > 0 (4.110) 

Proof ; The inequality (4.110) holds if and only if H(-jo), jw) ^ 0» all 
cj > 0. The positive semi-deflnlteness of H(”jw,jw) in turn follows from 
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2 > Ot A) obsarvAbl* and dat(Ju _I “ A " al) ^ At ^ 0« and 

thla conplatas eha proof. 

Noca thac InaquallClaa (4.74) and (4.110) bacoaw IdanClcal whan 
a ■ 0. For tha caaa a 0» wa can ragard (4.110) aa arising from (4.74) 
by raplacing ^ Ju by 4 J u-a . Indaad tha valuaa of tha matrix function 
for 0 ) ^ 0 on tha laft hand slda of (4.110) can ba obtalnad by avaluatlng 
tha RPDS ratum dlffaranca matrix along the a-Nyqulst contour shown In 
Fig. (4.4). For ^a slngla-lnput case, (4.110) slmpllflaa to tha following 
scalar Inequality 

|l + ■** ^ 1 Cl) ^ 0 (4.111) 

which admits tha same geometrical Interpretation as (4. 75) in that every 
point on tha Nyqulst diagram of ta(s) evaluated along the a-Nyqulst 
contour la away from the critical point by a distance of at least unity 
(sea Fig. (4.8)). 

Recognizing tha similarities between Inequalities (4.74) and (4.110) 

It seems natural to use the closed- loop system depicted In Fig. 4.9 In 
order to Interprete robustness Implications of (4.110). The following 
result can be readily obtained using inequality (4.110) and Theorem 4.6. 

Theorem 4.10 Given a RPDS with loop transfer matrix T^(s) and a 

prescribed degree of stability a , the perturbed regulator with loop 
transfer matrix T^(s) *T^s) L(s) will retain a degree of stability 
If 

(1) ^(s-a) and T ^^s-a) have the same number of unstable 

poles 

T^(s-a) has no pole on the Jb>-axls 


( 11 ) 


(4.112) 

(4.113) 
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1/2 -1 -1 /2 

(ill) Omax(n ' i, ” «) i " - I) < 1 for all « ^ 0 

(4.114) 

Note that condltloa amax(R^^^ l"^(J« - a) - p < 1 in (4.114) can 

be rewritten aa 

R L(J(a - a) L®(jw - a) R > R for all w ^ 0 (4.115) 

If both R and L(e) are choaen to be vllagonal. then (4.115) above slmpllflaa 
to the following eat of Inequalltlea 

Ra(l^(jw - o)) > I for all w > 0 (4.U6) 

and 1 * 1.2. ....a 

For caae of representing a pure gain change. It la clear from 

(4.116) that the pertxirbed regulator will retain a degree of stability a 

1 ^*1 

If each satisfies ^ 2 • Similarly, for can conclude 

from (4.116) that Che perturbed regulator will retain a dagree of stabi- 
lity a for < 60°. 1 ■ 1.2.....m. It Is interesting to note that 
these guaranteed margins with respect to degree of stability are identical 
to the guaranteed stability margins derived in Section 4.4.2. 

4.5 Properties of the RPDS Inverse Return Difference Matrices and 
Related Robustness Results 

We have noted in section 4.1 Chat useful robustness ^laracterizations 
(including the guaranteed gain and phase margins) for MIMO systems can be 
stated in terms of the minimum singular value of the Inverse return 
difference matrix. It was also mentioned briefly in section 4.2 that 
this quantity is a useful cosipleoent to the minimum singular value of 
the return difference matrix. The purpose of this section is to examine 





- 110 - 


i 


several Important properties of SPDS Inverse return difference matrix 
and the related robustness interpretations. Results in this section 
are presented in a fashion which closely parallels that of section 4.4. 


4.5.1 Properties in CoiMon with LQ Regulscors 

Onllke the case of the return difference matrices, no useful 
frequency domain characterisation of inverse return difference matrices 
is known for LQ regulators. However, when the control weighting matrix 
is chosen to be a scaled identity matrix, then a lower bound for the 
minimum singular value of the LQ inverse return difference matrix can be 
obtained in terms of the lower bound for the minimum singular value of the 
LQ return difference matrix derived in the last section. A precise 
statement of this result is given by the following theorem. 


Theorem 4J.1 Let the matrix ^ be the unique positive definite 
solution of the algebraic Rlccatl equation (4.59) with conditions on 
A, B, 2 K satisfied. Also assunm that 


and that 


Then 


det(ju I, - A) }* 0 for all w ^ 0 (4.117) 

^ (Ju) is invertible at all co > 0 (4.118) 

a -* 

, amin(I + T^(jw)) ‘ 

omin (I, + (jw)) 1 i + amin(I + Ta(jo)) (4.119) 

^ for all oj > 0 (4.120) 


Proof; The relation (4.119) follows from direct application of an 

inequality due to Nusman and Sandell [Nu 1]. To derive (4.120) from 
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(4.119), we make use of the Kalman's Inequality (4.80) derived In the 
leist section and this completes the proof. 

It follows from Corollary 4.1 and the above theorem that RPDS design 
with R chosen to be a scaled Identity matrix has a guaranteed gain margin 
given by (-j, j) and a guaranteed phase margin given by (-30°, 30°) . Both 
of these margins turn out to be mere conservative than those derived using the 
RFDS return difference matrix. 

It Is also possible to obtain these guaranteed stability margins 
for RPDS designs with R chosen to be a diagonal matrix. The derivation 
Is analogous to that of Theorem 4.5 and requires the use of the equivalent 
system depicted In Fig. 4.9. Based on the equivalent system, we can 
rewrite (4.1J9) and (4. 120) In Theorem (4.11) as 

. amln(I +fa(j«>»)) 

Omln(i+T^ (jw))^ 1 + amln(I +T^(jO>)) 

Using Inequality (4.122) together with Theorem 4.5 leads directly 
following result . Invertlblllty of ^(J^) for o) ^ 0 Is assumed. 


(4.121) 

(4.122) 
to the 


Theorem 4.12 Given a RPDS with loop-transf er matrix Tf, (s) . the 

loop-transfer matrix ^s) • (s) L(s) is stable If 


(1) T (s) and T (s) have the same number of unstable poles (4.123) 

"“Cl 

(11) ^(a) has no pole on the jw-axls (4.124) 

(111) amax (R^^^ L(ja») r"^^^ - P < f for all u > 0 (4.125) 


When R and L(s) are chosen to be diagonal, condition (4.125) In the above 


theorem reduces to 
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i - ll < 7 for all a» > 0 (4.126) 

Slid ^ 

We can now employ (4.126) to establish the guaranteed gain and phase 
ma r gi n s by assuming i^(s) to be real and i^(s} « e'^^ , respectively, 
for i • l,2,...,m. 

It is Important to emphasize that the guaranteed stability margins 
derived above and those In section 4.2 are margins that apply to all 
RPDS designs with diagonal R. The remarkable generality of these results 
In turn account for their conservative nature. Less conservative margins 
can be obtained If the actual values of the design parameters A, R 
and ^ are taken Into account In the robustness test. For a given RFDS 
design, the quantity defined by 

■ min amin(^ + ^(joj)) (4.127) 

Is in fact often greater than y. In cases when this Is true, we can 
replace condition (4.125) la Theorem 4.12 with a less stringent bound 


amax(R^^^ L(jw) - I) < B 


for all 0 ) > 0 


If R and L(s) are both chosen to be diagonal, then the Inequality 

(4.125) implies an improved guaranteed gain margin of (1 - B , 1 B ) 

-1 ®o -1 ®o ° 

and a guaranteed phase margin of (-2 sin y-, 2 sin y-) . 

It Is also noteworthy that guaranteed stability margins obtained 
from robustness tests based on amln(I^ *^2<j(J^)) cannot be further tightened 
by taking Into considerations the specific values of the design parameters. 
The quantity amln(X **’Tq(J<i))) Is always lower bounded by unity as a result 
of Tjj(s) being strictly proper. 


! 

i 

i 



-113- 


Conblnlng the guaranteed stability margins obtained In this section 
and the last, we arrive at the following Improved guaranteed margins that 
are valid for BPDS designs with R chosen to be diagonal. 

GM C Cl - (4.129) 

.1 

PM C (min(-2 sin ^ , 60°) . (4.130) 

(max(2sin“^ ^ , 60°) (4.130) 

(4.129) and (4.130) above represent the tightest guaranteed margins one 
can obtain using Theorems 4.4 and 4.5. This Is an example where appropriate 
combination of robustness tests can lead to a reduction In conservatism. 

4.5.2 The Effect of the Stability Factor o on the RPDS Inverse 
Return Difference Matrices 

If 

det(j0)T - A) # 0 for o» ^0 (4.131) 

and 

det(T^(j<o)) 0 for w>o, (4.132) 

we can rewrite the RPDS frequency domain equality (4.63) In the following 
form 

(I + Ta"^(io»))® R(I + T^^(jw)) 

[R + H^(-jw,ja»)] Tj^"^Cj&») for cu > 0 (4.133) 

Recall from the derivation of Corollary 4.2 that the effect of a on 

B 

(i+ ^ Ow)) R(l. + ^(jw)) is readily obtained by considering the 
effect of a on ^(-jai,ja!>. obtain a similar characterization for 
(I + T, ^(-Jw)) R(I + X,, ^(jcu)), we need to consider the effect of a on 

— ^ “Tjt 
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(-8.3) - T ’•'^(-s) CR + ^ (s) ) T “^(s) (4.134) 

• K (-81 - B]"'^ (R + B^(-sI - (£ + 2kxK )(al - A)"^ B)J 

» "tX ” 

K (8l - A)"^ B] (4.135) 


;rttere ^ is the positive definite solution of the RPDS algebraic Rlccati 

equation C4.59)> The coefficients of both rational matrices T (s) and 

TX 

^(-s.s) are dependent on the stability factor a through the matrix , 

which is known to be an increasing function of a (Lemma 3.1). Unfortunately, 

^ enters the matrix H^(-s,s) in a highly nontrivial fashion and thus 

makes it difficult to precisely characterize the effect of (x on H (s) . 

~o 

When R is chosen to be a scaled identity matrix Bl. we can approximate 
the RPDS Inverse return difference with matrix with 


I +^^(Jo»))“ (I + 4^ao»)) 


/ J,KB .\/ B'KB ,\ 


T 

B K B 
>6 ^ 


I + cu^ (B^ K B)"^ 
— — -a ~ 


(4.136) 


at values of to where jO)_I dominates the matrix (jco^ - A) . Using the 

incremental property of K (Lemma 3.1), we can readily show that 

—a 

2 2 T -2 

^ 4 - 3 (u (B ^ ^) 2 monotonlcally decreasing function of a (l.e. 

^+8o»(^^B) <I +0^(o^(B^ K B) ^ for a. > a, and for all 
w ^ 0). It follows from this observation that (^ + (8*“* 

amin(^ ^^(J<o)))ls a decreasing function of a in the high frequency 
region. Attempts to characterize (1,+ T^(jaj)) (1. + ^ (J«)) at frequencies 

where the approximation (4.136) falls to hold has not been successful. 
However, it is our conjecture that (I + t”^(Jco))® (I + T (jw)) will 

— —a — —a 
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generall7 behave as a decreasing function of a at all w > Q for R 

chosen to be a scaled identity siatrix. If this is Indeed the case, 

then two important design implications follow. First, the quantity 

defined in C4.127) also becomes a decreasing function of a As 

a result, the guaranteed downward gain margin given by 1 - 0^ will 

2 

deteriorate with increasing a. Second, the maximum RPDS system 
bandwidth Umax defined as the highest value of u (see [Ch Ij section 
7.7) at which 

omin(I + (jo»)) -1 (4.137) 

will increase monotonlcally with a . This can be readily verified by 
inspection of the amln(^ + (ju)) plots of a RPDS design (Fig. 4.13) 

with the minimum singular value of the inverse return difference matrix 
being a decreasing function of a. 

4.5.3 Properties of the RPDS Inverse Return Difference Matrices on 
the g-Nyquist Contour » 

The results presented in this section are stronger versions of 

those considered in section 4.5.1. All the margins derived here apply 

to the stability as well as the degree of stability property of the 

underlying RPDS. To begin with, we shall consider RPDS designs where 

R is chosen to be a scaled identity matrix and obtain a lower bound 

^ See section 5.4 for an example where amln(J ^^(Ju))is a decreasing 
function of a for all u ^ 0. 

2 

Recall from Example 4.1 that the same result may not hold for the 
actual stability margins. 


Magnitude (dB) 
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FIG. 4.13 Singular Value Plots for RFDS Design; 
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for the quantity omln(I + ^^(jw-a)) that is valid for all o» > 0. A 
precise statement of the result is given by the following theorem. 


Theorem 4.13 Let the matrix ^ be the unique positive definite 

solution of the algebraic Rlccati equation C4.59) with conditions on A, 

^ satisfied. Also assume that 

det((jM -oil - A) 0 for all <o ^ 0 (4.138) 

and that 


T^(-a + j(o) is invertible at all a» ^ 0 


(4.139) 


Then 


omin(I, + ^^(-a + jcu)) 


amin(JC + (-« + j^^)) 

— 1 + amin(X + (-a + jw) ) 


for all 01 > 0 (4.140) 


Proof ; We shall prove this result from a standpoint different from 

that adopted in Section 4.4.4. Recall from the proof of Theorem 2.5 
that we Introduced a time- invariant LQ problem on our way to obtain a 
solution for the RFDS problem. It was also shown that the desired RFDS 
feedback gain is Identical to that obtained from the related LQ problem 
where the objective is to minimize 


11m / 

t I 

1 •'t 


^ 1 *(t) + vJ(t) R u(t) Jdt 


(4.141) 


aubject to the dynamic constraint 

x(t) “ (A+ap x(t) + B u(t) 


(4.142) 
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The controllability and the coat observability of this LQ problem 

follows from those of the original RFDS problem and the RPDS state feed- 

-1 T 

back gain Is given by • R 1. where Is the positive definite 
solution of the RPDS Rlccatl equation (4.59) . Applying the result of 
Theorem 4.10 (which also holds for LQ regulator problems with R chosen 
to be a scaled Identity matrix) to the solution of the above LQ problem 
leads directly to (4.140) and this costpletes the proof. 

Based on the equivalent system given In Fig. 4.9, the following 
Inequality that holds for arbitrary R (l.e. R not necessarily be diagonal) 
can be readily obtained 

omln(i + R ^^(j« - a) 2 for all (o ^ 0 (4.143) 

The proof of (4.143) Is similar to that given In Theorem 4.13 with 
^(s) replaced by R^^\j(s-a) result below la a direct 

consequence of (4.143) and Theorem 4.13. Invertlblllty of ^(Jomx) for 
cu ^ 0 Is assumed. 

Theorem 4.14 Given a rational transfer matrix ^(s) of a RPDS 

system with a prescribed degree of stability a, the perturbed system with loop 
transfer matrix given by ^(s) ■ ^ degree of stability 

a if 

(1) ^(s-a) and T(s -a) have the same number of unstable poles 

(4.144) 

(11) Xv(>~o) has no pole on the jcu-axls (4.145) 

(111) amax(R^/2 ^-1/2 _ p < 1 for all w ^ 0 (4.146) 

If both R and L are diagonal, then (4.146) reduces to the 
following set of scalar Inequalities 
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lAj^Cjw - a) - l| < Y for all u ^ 0 

and 1 • l,2,...m 


(4.147) 


For case of representing a pure gain change. It is clear from 
condition (4.147) that the perturbed system will retain a degree of 




1 3 ■'^1 

stability o for satisfying J ^ ^ T* Similarly, for “ ® 


the condition (4.147) implies that the perturbed system will not lose 

the degree of stability for |(^^| < 30°. A RPDS design with R chosen 

1 3 

to be diagonal therefore possesses a guaranteed GM^ of (yty ^ ® 

guaranteed FM^ of (-30°, 30°). These margins are found to be more 
conservative than those derived from (4.116). 

To characterize the margin for a given RPDS design, we place 
condition amax(R^^^ t,(Jw ” <*) ” i) ^ 2' (4.142) by 


amax(R^^^ L(jo) - a) - I) < 8a 


(4.148) 


where 8a Is defined to be 


-- 1 . 


8a ■ min amln(I + T_*(1ai - a)) 

o»0 “ 


(4.149) 


This quantity 8a turns out to be greater than j for most RPDS designs, 

and thus enables us to obtain an Improved guaranteed Of of 

-1 ®o -l®o 

(1 - 8a, 1 8a) and an Improved guaranteed PM^ of (-2 sin -j- • ^ ) ' 

Combining the guaranteed CM and PM obtained here with those derived 

a a 

In section 4.4.5, we arrive at the following guaranteed CM and PM that 

a a 

hold for a given RPDS design with diagonal R 


CM D (1 -go* «) 
a 


(4.150) 


PHx D (min(-60°, -2sin‘^ |SL) (60®, 2sin"^ |£.) 


(4.151) 
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These are the best margins that one may derive using Theorems 4.10 and 
4.14. 


4.6 Roll-Off Requirements at High Freouencles 

In our discussion of RPDS robustness properties thus far, we 
have neglected the important design consideration of having the nominal 
loop gain sufficiently attenuated over the frequency range where the 
magnitude of the multiplicative perturbation L(Ju) becomes large compared 
to unity. It Is a physical fact that the quality of a nominal design model 
Inevitably deteriorates at high frequencies as a result of unmodelled 
and/or unknown dynamics of various types. In the multiplicative form of 
perturbation L(s)^ this means that cmax(L(JO))) will assume value close to 
1.0 at low frequencies but will increase to 2 and beyond at high frequencies. 
In the face of such uncertainties, the seemingly excellent guaranteed RFDS 

stability margin Is clearl*' Inadequate. In the high frequency region, 
the + 60^ phase margin is of no value since a neglected time delay 
will ultimately produce phase error In excess of 180°. Moreover, the 
return difference inequality (4.74) from which the guaranteed RPDS 
stability margins are derived does not hold at high frequencies since 
the frequency response of any physically realizable system must have a 
roll-off at a rate greater than or equal to a . The roll- 
off requirement of RPDS (and for LQ regulator In general) In the face 
of uncertainties can be derived from the basic Inequality (4.27) In 
Theorem 4.5. 

At frequencies where all the feedback loops of ^(Ju) are rolled 
off, the quantity amax(T (Ju)) ■ amin~^(T~^(Ja>)) becomes small. As a 

“tX 
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result, we can approximate the minimum singular value of ^ es 

o mln(l + (j w) ) ■ omin(T^^ (jw) ) 

■ omax"^ (^(jw)) (4.152) 

Substituting this approximation into (4.27) gives 

onax(L(ja») - ^ amax(T^(Jw)) < 1 (4.153) 

It is now clear from (4.153) that the loop tranafer matrix Tq^(Ju) has 

to attenuate faster than omax ^(L(Ju) - ‘^) in order to satisfy the stability 

requirements at frequencies where the physical process is ill-represented 

by the model. But the roll-off rate of Tq(J(») is limited to <u ^ as 

evidenced by the following approximation 

T 

B*K B 

T CJw) - (4.154) 

nx JW 

which holds for values of O) sufficiently large. Substituting (4.154) 
into (4.153) we obtain 

caax(L(Jw) - 1) (jmax (i^ ^ (4.155) 

It follows from (4.155) that the cross-over frequency of TQ(j(«>) has to 
be located well below the frequency where omax(L(Jo)) - ^starts growing 
large. 

Equation (4. 155) also makes explicit the relationship between the 
roll-off frequency of ^(Jo)) and the stability factor a . If we define 
the cross-over frequency (a amax of T^Cju) to be the frequency 

when 

Oinax(Tft(j<»i)) ■ 1 (4.156) 

T 

it can then be readily shown that the quantity 0max(^ ^ B) provides 
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an upparbound foroioiaax. Moraovar, wa can deduct fron Lena 3.1 that 
Oiaax(B^ i* *n Increaalng function of a. The lapllcatlon of the roll- 
off raqulramant (4.155) on the choice of atablllty factor a nov baconea 
clear. While a large value of a will Inprove the syaten'a apeed of 
reaponae, it aay alao extend the bandwidth Into reglona ^ere LCjJU) la 
large. In picking a value for the atablllty factor a , the dealgner 
needa to find a compromlae between the requlreawnta for atablllty In 
face of uncertalntlea and apeed of perforaance. 

4.7 Sunarv of RPDS Robuetneae Prooertlea 

We ahall auaaarlze in thla aectlon the varloua robuatneaa 
propertlea of RPDS dlacuaaed In the two prevloua aectlona. 

(1) Guaranteed gain and phaae ear;^na for atablllty 

If the matrix R la choaen to be diagonal* then a RPDS 
dealgn poaaeaaes the following guaranteed atablllty marglna 

GMD(0j,,-) (4.157) 

PM D (nln(-60®, -2aln'^ ~ ), (4.158) 

aax(60®, 2aln”^ ^)) 

where 8^ la given by (4.127) 

(11) Guaranteed gain and phaae marglna for the prescribed degree 
of atablllty. 

If the matrix R la chosen to be diagonal* then a RPDS design 
possesses the following margins with respect to the 
prescribed degree of stability 
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D (8a . - ) C4.159) 

D (jain(r2«in ^ , -60®) 

8 ^ 

m«x(2«ln"^ ,63®)) (4.160) 

(ill) Eff«ct of ch« scabillty factor a on the robustnosa propertlca 
of BPDS 

The RPDS's ability to tolarata aodalllng error (quantified In 

tenne of the magnitude of the modelling error In question) Improves with 

increasing value of the stability factor a only for some special choices 

of error representation. Among the four types of modelling error presented 

In this thesis, we only obtain an Improved tolerance for 

^(s) ■ (T ^(s) - T ^(s)) T(s) with Increasing a . The guaranteed RPOS 

margins summarized In (1) and (11) above may deteriorate with increasing 

a If 8 And 6 turn out to be a monotonlcally decreasing functions of a 
0 a 

(sec section 5.4 for an example of being a decreasing function of a). 

It was shown In section 4.4.4 that the behavior of the actual RPDS 
gain and phase margins with changing value of a Is highly system dependent* 

(Iv) High frequency roll-off requirement 
T 

The bound amax(B ^ B) for the maximum cross-over frequency of 
RPDS Is an Increasing function of a (see section 4.b>. This In turn Imposes 
an upper limit on the value of a that we nay coqiloy In RPDS design* 

It Is well known to control system designers that too large a cross-over 
frequency may result In excitation of the umsodelled and/or unknown dynamics 
at hJ.gh frequencies which la undesirable from a stability robustness point 
of view. 
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It Is important to emphasize that the conclusions summarized 
in this section are only valid with respect to the types of model error 
under investigation. The relative error in T^^Cs) is particularly useful 
for our purpose, because its associated robustness theorems Clheorems 4.4 
acd 4.6) allow simple derivation of BFDS multi- loop margins with respect 
to stability and degree of stability. Further study on the ability of 
SPDS to tolerate model errors other than those considered here (see 
[Le i.i Section 3.9 is needed for a more complete understanding of 
its robustness properties. 

4.8 Robustness Properties of Kalman Bucy Filter with a Prescribed 
Degree of Stability (FTPS) 

As we have noted in Chapter II, KBF with a prescribed degree of 
stability (FFDS) is the mathematical d\ial of RFDS. Dual robustness 
results are therefore obtainable for such designs. These robustness 
properties ensure the nondivergence of the filter under variation in the 
nominal model of the plant which is to be estimated. In section 4.8.1 
we set up the framework for robustness analysis of FFDS. The robust- 
ness properties for this class of filter are discussed in sections 
4.8.2. to 4.8.4. 

4.8.1 Formulation of the Robustness Problem for Kai" «" Filters 

with a Prescribed Degree of Stability 

The basic FPDS problem considered here is identical' to that 
discussed in Theorem 2.5. The underlying linear system is given by 


x(t) 


A x(t) + i(t) 


(4.161) 
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i(t) -C x(t) + iCt) C4.162) 

where ^(t) and ^(t) are zero mean white noise processes with spectral 

intensity matrices H and 6 respectively. We also assume the observ- 

1/2 

ability of (£, A) and controllability of (A, H ) . The state estimate 
x(t) is then specified by 

xCt) - A x(t) + v(t) (4. 

v_(t) " 2.(t) “ C. x(t) (4.164) 

where 

^ " hx (4.165) 

and ^ is the unique positive definite solution of the Riccatl equation 

(A + aI)Ecj + la + 1 - C- la " ^ (4.166) 

A useful method for describing the state estimate dynamics of 
FPDS (in fact for KBF in general) is given by the following set of 
feedback equations 


e(t) “ A eit) + w(t) 

(4.167) 

s.(t) « C e{t) - e.(t) 

(4.168) 

w(t) - -l^a(t) -i(t) 

(4.169) 

where s(t) ■ i(t) - x^t) 

(4.170) 

and 


e(t) ■ x(.t) - x(t) 

(4.171) 


A block diagram representation of these equations is given in 
Fig. 4.14. By Ignoring the noise sources, this can be rearranged into 

T 

t a unity feedback of the type considered in Section 4.2. It is now readily 

i 

i 
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Opan-loop Error Oynamict 



Fittir Coin 


FIG. 4.14 Feedback Representation of the Error Dynamics of FPDS 


Perturbed Error Oynamics ^(s) 


r 1 



FIG. 4.15 Feedback Configuration for Robustness Axxalysls 

of FPDS Error Dynamics 
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apparent that 

F^(a) A C(al - A)“^ (4.172) 

the loop transfer matrix of the error dynamics loop of KBF with prescribed 
degree of stability Is the dual of T^(s) for RFDS. The discrepancy 
between the linear model employed In the design of the filter gain and 
the actual system dynamics are modelled as multiplicative perturbations 
L(s) inserted In the closed-loop after C_ (see Fig. 4.15). 

4.8.2 Common Robustness Propei-uies with KBF 

Based on the feedback representation of the FFOS error-dynamics, 

we are ready to characterize the stability margins of KBF with a prescribed 

degree of stability. The nominal error-dynamics model corresponds to 

the Ideal situation In which the linear design model represented by (4.161) 

and (4.162) are exact. Since all the results developed In this section 

are mathematically dual to the BFDS results obtained In sections 4.4 

and 4.5, the respective derlvatloi^ will be omitted. The main emphasis 

will be on the Interpretation of robustness results In the estimation 

* 

context. All the FFOS robustness properties presented In this section 
are known properties of KBF, and are Included here for completeness. 

To begin with, we state two frequency domain equalities for FFOS. 
These are the exact dual of (4.63) and (4.64) In Theorem 4.8 and are 
basic to the derivations of the robustness results given In this section. 

Theorem 4.15 Let the matrix ^ be the unique positive definite 


solution of the algebraic Rlccatl equation 
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+ +H -Z c'^e’^c^ - JL (.4.173) 

with 

(i) H > 0 and 0 > 0 (4.174) 

Cii) (A, controllable (4.175) 

(ill) (£, A) observable (4.176) 

then 

(£ + ^(-a+s)) 0(i + ^ (-a-a))*^ -0+M(s,-s) . (4.177) 

and 

<i+4(s)) 0(1 •^ 4 (-s))'^ *0 + M^(a,-s) (4.178) 

where we define 

M(5.a) • C(51 - al~ £(sl - a£ - a)'^ (4.179) 

and 

M ( 5 ,s) - CO; I - A)“^ ( 2 aT +5) (si - A)'"^ ^ (4.180) 


Remark The equalities (4.179) and (4.180) are derived from the 

FPDS algebraic Rlccatl equations using manipulations similar to those 
employed In the proof of Theorem 4.8. 

Remark The two frequency domain equalities stated in Theorem 4.14 

correspond to the two Interpretations of FPDS given In section 2.4. 

Using the positive deflnltness of ^ and equation (4.180) we 
can readily obtain the following corollary to Theorem 4.15. 

Corollary 4.9 If £jls the unique positive definite solution of 

the FPDS algebraic Rlccatl equation (4.173), with the respective require- 
ments on A, £, H, and 0 being satisfied and If det(Jo>I - A) 0 for 
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all (I) > 0 , we have 

(1+4 (Jw'). e(i + > 0 for all « > 0 (4.181) 

Combining inequality (4.181) with Theorem 4.4 readily leads to the 
following robustness theorem which holds for KBF in general. 

Theorem 4.16 Given a FFDS with loop transfer matrix of error 

dynamics given by l^(a)> the mismatched error d 3 mamlcs with loop transfer 
matrix Fa(®) ■ t,(®) Za(®) ^ closed-loop stable if 

(i) Fg(s) and Fq^s) have the same number of unstable poles (4.182) 

(11) FjjCa) has no pole on the jco-axis (4.183) 

Clii) omax (0"^^^ L"®(ja)) - i) 1 1 for allw^O (4.184) 

The condition (4.184) describes the inherent robustness properties 
of the FFOS design procedure. It says that every FFDS design can tolerate 
at least multiplicative perturbation L(s) satisfying the bound (4.184) 

If the observation noise at each output channel are uncorrelated (l.e. 

0^ is diagonal) and if the model mismatch can be represented by a 
diagoual multiplicative perturbation of the error dynamics, condition 
(4.184) can then be Interpreted in terms of the gain and phase margins 
of each output channel in the feedback representation of error dynamics 
(Fig 4.15). 

The derivation follows from recognizing the equivalence between 

omax(e"^^^ lT®(ju) - i) i 1 for all w ^ 0 (4.185) 

I 

It can be readily shown that Theorems 4.4 to 4.7 also apply to model 
error described by T(s) ■ L(s) T(s) 
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and 

0L®(Jci)^ + L(Jw) 0-9^0 for all (u ^ 0 (4.186) 

When £and L(s) are both diagonal, (4«182) further simpllfiea to 

2Re(jL(jw)) > 1 for all oi > 0 

^ ” i - 1,2 (4.187) 

If f>^(Joid Is real, then (4.187) becomes 

\> Y (4.188) 

Alternately, if llj^(j(o) • e , then condition (4.187) becomes 

< 60° (4.189) 

The conditions (4.188) and (4.189) can be Interpreted aa Implying that 
FFDS design employing uncorrelated observation noise leads to a guaranteed 
gain margin of (-|-, «) and guaranteed phase margin of (-60°, 60°) In each 
output channel of the error dynamics feedback system (Fig. 4.14). 

These margins are relative to the Ideal situation that (4.161) and (4.162) 
are exact. 

It Is Important to stress that the guaranteed margins thus derived 

holds for every FFDS design using diagonal £ . The generality of this 

result In turn accounts for Its conservatism. Less conservative margins 

can be obtained for a given FFDS design If we combine the bounds derived 

above with those derived using Theorem 4.5 and Corollary 4.3. The resulting 

guaranteed gain and phase margins obtained from the latter are given by 

Y Y 

(1 - Yq» 1 Yq) and(-2sin“^ 2sin”^ ^ ) respectively where Y^ Is 
defined to be 

Y - min omind + 0"^^^ (4.190) 

° 01^ - - -a 
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and 0 Is diagonal. 

Since the value of usually lies between and 1, the down** 
ward gain oargln 1 ~ results from application of Corollary 4.3 Is 

therefore always less than or equal to -j. However, the related guaranteed 

—1 Y —1 Y 

phase margin (-2sln 2sln~'^ •^) will extend beyond the Interval 
(-60°, 60°) only when > 1. 

It Is clear from the above discussion that the following guaranteed 
01 and FH hold for error dynamics of a given FFDS design with 6 chosen 
to be diagonal 


GM D (1 - Yg. •) (4.191) 

Y Y 

PM D (min (-60°, 2sln"^ ^) , max (60°, 2sln“^ ^)) (4.192) 

(4.191) and (4.192) Indeed yield the least conservative guaranteed margins 
one can obtain using Theorem 4.4 and 4.5. 

4.8.3 The Effect of the Stability Factor a on the Robustness Properties 
of RPDS 

Like the case of RPDS, the effect of a on the robustness properties 

of the FPDS error-dynamics can be characterized In terms of the effect 

of a on the matrix function I + (*) ^ ^(s) . 

— -a “ “U 

The result described In the following corollary to Theorem 4.14 
makes precise the behavior of 1, -f (s) as a varies. It Is the exact 
dual of Corollary 4.2. 

Corollary 4.10 Let the matrix^ and be the unique positive 

solutions of FPDS algebraic Rlccatl equations 
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(A+ctI)r (A +01)’^ + 2 - -0 

-®1 ~®1 - -Oj^- — “ 

(4.193) 

and 


(A + aI)^ + =-^ - 0 

=Oi -»2 - -ttj - -»2 

(4.194) 

where 


(1) 5 ^ 0, e > 0 

(4.195) 


(11) (£. ^ observable 

(4.196) 


(111) ^ observable 

(4.197) 

and 

(Iv) det(jwl - A) j* 0 for O) ^ 0 

(4.198) 

Then one 

has (I + F (j«))0 a+ Fj OT)® 
^1 1 



> (I + F (Jw))0 (I + F„ (jw))® for all <o > 0 

“ “®2 ” ” ^2 

if <*2 ^ ®2 — ® 

(4.199) 


When £ is a scaled IdenclCy matrix (l.e. the noises at each output channel 
have the same Intensity and are uncorrelated) » the following Inequality 
on the minimum singular value of FPDS can be readily derived from (4.199) 
by using the properties of singular values. 

omin(i+ > ominCi + F«®(Jw)) for all w > 0 (4.200) 

1 2 ■“ 

This Inequality provides us with a useful way of assessing the effect of 
a on the robustness properties of the FPDS designs. Using condition (4.22) 
of Theorem 4.4 together with condition (4.198) leads directly to the 

conclusion that FPDS tolerance of model mismatch represented by^ 

E^(s) - (F^^(-8) - F^Vs)) l^(-s) (4.201) 

^ The singular values of E.(s) given in (4.201) will be different from 
those of defined by E^(s) ■ (F^^(s) ” 
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Isproves with ot . As we have noted before, error neasare of this type Is 
unintuitive and difficult to Interprets. 

Like the case of RFDS. a precise characterization of (jcu) 

with respect to change of a Is not available. However, when 6 Is an Identity 
matrix. It can again be shown that <7mln(^ + ^ Is a decreasing 

function of a for sufficiently large value of at. 

4.8.4 Robustness Properties with Respect to the Degree of Stability a 
We shall first examine the robustness Interpretations of eq\iallty 
(4.179) for FFDS. The following corollary to Theorem 4.15 that characte- 
rizes the behavior of the FPDS return difference matrix on the a-Nyqulst 
contour Is the exact dual of Corollary 4.8 

Cor ollary 4.11 Let ^ be the unique positive definite solution 

of the algebraic Rlccatl equation (4.163), with the respective requirements 
on A, £, £ and ^ being satisfied and det(JO)l. - a^) ^ 0 for all o> ^ 0. 

Then 

(I + p;j(-a + j(o))0 (I + Iij(-a + jw))® > e (4.202) 

Working with the equivalent feedback representation of the FPDS error 
dynamics given by (4.167) to (4.171), the theorem given below Is a direct 
consequence of (4.202) and Theorem 4.5. 

Theorem 4.17 Given a FFDS with a prescribed degree cf atablllCy a and a loop 
transfer matrix P^(s) . The mismatched error-dynamics with loop transfer 
matrix |^(s) - L(s) ^(s) has a degree of stability a If 
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(i) F(s - a) and F(s - a) have the sene number of unetable 

poles (4.203) 

(11) ^(s - a) has no pole on the ju-oxls (4.204) 

(111) onax(9“^^^ “ a) - I) < 1 for all w > 0 (4.205) 

If the margin with respect to the degree of stability a derived 
from (4.205) Is combined with those derived using Theorem 4.5 end 
Corollary 4.3, we obtain the following Improved guaranteed margins that 
apply to a given FPDS design employing uncorrelated .observation noise. 

St each output (4.206) 

channel C " Y® » ") 

FM St each output (4.207) 

-1 Ya 

channel C (min (-60®, -2sln“ y- ) , 

. Y 

max (60°, lain" -^ )) 

where Yq defined by 

Y„ ■ win omlnd + T + j«)) (4.208) 

0)^0 " 

4.8.5 Concluding Remarks 

As we have coimnented before (see section 2.4), exponential weighting 
of data Is a technique well known to filter designers for curing the filter 
divergence problems. The general thinking Is that this will prevent the 
old data from saturating the filter. In section 4.8.2 to 4.8.4, Intuition 
of this type Is subject to rigorous examination using the recently developed 
results In robustness analysis for MIMO systems. The results obtained 
here however apply only to those coses where mismatch between the actual 
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and Che design model can be represented as a nulclpllcatlve perturbaclon 
ac Che output matrix £ in Che feedback representation of the filter error 
dynamics (Fig. 4.15). 

It is clear from the discussion of the three previous subsections 
and Che summary remark in section 4.7 for the dual results of RPOS, chat 
improvement of FPOS's ability to tolerate model mismatch of Che type 
depicted in Fig. 4.15 only occur for very special choices of model error 
representation. The exponential weighting technique la therefore not 
a universal cure for every possible type of divergence problem. The 
insights obtained from this section help to identify situations ^ere 
such technique can be effective. 

In many applications of Interest, the above conclusion nay be 
excessively conservative. This follows from the fact that we have only 
used some magnitude information on the model mismatch in our robustness 
analysis. In the case where structural information is available, the 
results In Chapter 4 of [Le 1] may be applicable. 
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CHAFTER V 

APPLICATIOS OP THE REGUIATOR WITH A PRESCRIBED 
DBG8EE OF STABILITY TECHMIOnE TO STATE FEEDBACK 
DESIGS FOR MDLTI-TERMINAL DC/AC POWER SYSTEM 

5.1 Introduction 

The purpose of this chsptsr Is to denonstrste the sppllcstlon of 
the BPDS technique to design of s state feedbrack control law for a 
9-aachlne, 4-termlnal DC/ AC power system. Based on this 
example, the various robustness properties of SPOS designs discussed In 
the previous chapter will be lUustrated. 

Two versions of the 9-machlne^ power system model are available 
around 5 operating points. They differ In their details of machine 
description. In the simple version of the model, every generator Is 
represented by a second-order classical machine. This gives a total 
system order of 18. The open- loop poles of the system are plotted In 
Fig. 5.1. F- ech of the 16 complex poles corresponds to a mode of Inter- 
machlne oscillation. Frequencies of such oscillations range approximately 
from 0.3 to 1.0 Hs* The two remaining poles are located on the negative 
real axis. The one located at the origin corresponds to the mode of clock 
error. The mode located at .375 Is the mode of average frequency. 

^ The power system models employed here are generated by Sherman Chan 
using the Posslm progrm of the General Electric Company 
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Readers Interested In the physical Interpretation of these modes may 
consult Chapter 3 of [Ch 1]. 

In the detailed version of the model, two of the generators are 
modelled with four rotor circuits (two In the direct axis and two In the 
quadrature axis). They are also equipped with IEEE Type 1 exciters and 
third order power system stabilizers. The total order of this systeo 
is 38. The open- loop eigenvalues of this system at operating point 1 
are plotted In Fig. 5.2. The complex poles of the system that are 
associated with machine oscillations are nximbered In descending order 
of their frequency.^ Comparing the oscillatory modes given by the two 
models show that they are In good agreement. It Is also noteworthy 
that 6 ,<alrs of the oscillatory modes have no counterpart In the 18-state 
model, because they are associated with the stabilizing components of 
the power system. 

The basic design objective is to move the open-loop oscillatory 
poles of the system to an appropriate region In the left-half complex 
plane. This region is determined approximately by using engineering 
Judgement on how large the closed-loop bandwidth may be without allowing 
unmodelled high frequency disturbance or dynamics to destabilize the 
system. For physical considerations, the average frequency mode and 
the clock error mode are to be kept Intact. 

In section 5.2, several full state feedback designs are obtained 
by using the RPDS methodology. A different value of the stability 

Only one pole of each complex pair is displayed in Fig. 5.2 . 


1 
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i 



FIG. 5.2 Open-Loop Poles of Che Detailed Model for the 

9-Machine Power System Model at Operating Point 1 
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factor a is employed for each design. The resulting closed-loop systems 
are found to display some desirable properties v^lch are not predicted 
by the theory developed In the previous chapters. 

The robustness properties of the RPDS design are evaluated In two 
ways. In section 5.3 we apply the control laws obtained for the simple 
model at operating point 1 to the other four operating points. The 
resulting movements of the closed-loop poles provide a good Indication of 
the system tolerance to change In operating conditions. The control law 
obtained for the simple design model are also Implemented on the detailed 
38 state model at various operating points . This allows us to study the 
effect of the unmodelled dynamics due to the exciters on RPDS designs - 
an Issue that we feel to be Important since In this case the neglected 
dynamics are no longer separated from the Intermachlne dynamics. 

In section S.4, we turn to the frequency domain robustness analysis 
formulated In Theorems 4.4 and 4.3 with emphasis given to the effect of 
the stability factor a on the quantities omax(l 't'T (Jco)) and 
omaxOL + T^“^(jcu)). 

5.2 RPDS Designs 

The design parameters to be chosen In a RPDS problem are the 
stability factor a t the state weighting matrix ^ and the control 
weighting matrix R. In the absence of Information concerning the relative 
cost of control for the different DC- terminals, we shall simply pick R 
to be an Identity matrix. The state weighting matrix 2. Is chosen to 
satisfy the following objectives: 



-141’ 


(1) Only Che oscillatory modes associatod v*lch the Intecmachlne 
oscillation but not the modes associated with the real poles 
ace penalized. This is done to limit the control amplitude 
and to avoid interaction between the multi- terminal DC 
controller and ezlsclng mechanisms for correcting the 
generation mlsmacch and clock errors. 

(il) The damping ratio of the oscillatory modes should be around 
0.2. This choice is somewhat arbitrary although it serves 
to provide a reasonable LQ regulator (a « 0) design to 
start from. This is not the only method of specifying 
the £ matrix in relation to transient response requirements. 
Indeed, a suitable choice of ^ has to be considered Jointly 
with the stability factor a. 

In cases where a - 0 (i.e. the LQ problem), the previous require- 
ments on ^ can be easily met by application of modal weighting techniques 
such as Solheim’s method [So 1] •Whena is nonzero, an additional trick 
is required to prevent the two real modes from being moved under feed- 
back. To make clear the underlying problem we consider the case where 
a is picked to be 0.5. This choice of stability factor will result in 
a matrix (A + ap with 2 unstable real poles . If the matrix £ Is chosen 
to make these poles unobservable, then it follows from the properties 
of algebraic Rlccatl equations and the discussion in Chapter 3 concerning 
the construction of RPDS root-loci that the resulting RPDS feedback law 
will shift the average frequency mode to -0.625 and the clock error mode 
to -1.0. A scheme to avoid this problem is outlined below. ^ 

^Thls scheme was first suggested to the author by N.A. Lehtomaki 
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Step 1 Compute a positive seml-deflnlte matrix ^ which Is defined as 
S, ^ “6 (£ 3 ^ + *2 ^2^ 

vitere are the right and left eigenvectors respectively of the 

average frequency modeiX^ X 2 right and left eigenvectors respecti- 

vely of the clock error mode and 8 Is set to be any positive real number 
larger than a. 

Step 2 Pick a state weighting matrix £ that makes the two real poles 
(A +• SJ cost-unobservable and the closed-loop oscillatory poles possessing 
a damping ratio about 0 . 2 . 

Step 3 Using A * A + S In place of A In the RPDS design (l.e. solution 
of the RPDS algebraic Rlccatl equation). The resulting feedback control 
law Is given by 

u(t) - - / x(t) (5.2) 

where Is the unique positive seml-deflnlte solution of the equation 

K (A + aI)+A + aI)K-KBB'^K+2 » £ (5.3) 

-ex'— — ' — — -a -o -a “ 

It can be readily verified that the procedure described above 

results In feedback laws that leave the clock error mode and the average 

frequency mode unchanged. First* note that the matrix A has two real 

poles located at -0 and - $- 0.375 with their left and right eigenvectors 

given by ^2 *2 Moreover, by the choice of 6 , 

these two poles will remain In the left-half complex plane upon addition 

of 01 to their real parts (as Is done In the RPDS design with stability 

factor given by a ) • It then follows from the choice of the matrix £ 
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and Che properties of algebraic Rlccacl equations that the matrix 
( A *!■ ~ 1, ^ ) which Is obtained by applying the state feedback law 

(5.2) to x(t) " (A ap x(t) will have two real poles located at 
-3 •t'a and -0.375 -6 -ha. Moreover the eigenvectors of these two poles 
are Identical to those of the clock error mode and the average frequency 

- ij - 

mode respectively. Subtracting ^ from “ ®. 1 . will return 

these real poles back to their original positions (-6 - 0.3.75 to 

-.375 and -6 + a to 0) with their corresponding elvenvectors remaining un- 
changed. It should be emphasized chat Che matrix ^ Is used only for the 
purpose of computing the feedback law. 

The simple model of the 9-machlne power system at operating point 1 
Is chosen to be the nominal model for design. Several values of ci 
ranging from 0.0 to 0.6 are tried In the design. These values of are 
compatible with the damping rate observed In the actual power systems. 

Too large a value of cx will result In faster performance at the cost of 
having Inputs with unacceptably large magnitude. 

The closed-loop poles of Che resulting RFDS design with <x equals 0.0> 
0.2, 0.4 and 0.6 are plotted In Fig. 5.3. Given the Increment of ci equals 
to Aa, all the oscillatory modes move further out to the left by an 
amount roughly equal to -Aoi. This Is by no means a property common to 
all RPDS at all values of positive a. Such behavior of RPDS closed-loop 
poles Is probably a result of the fact that all modes to be controlled 
in the multiterminal DC/AC powe system are of the same nature (i.e. they 
all represent modes of intermachlne oscillations). It Is also note- 
worthy that changes in the Imaginary part of the closed-loop poles with 
respect to a are negligible compared to changes In the real part. 
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Mode 

Open 

Loop 

- 0.0 

- 0.2 

- 0.4 

- 0,6 

1.2 

-0.18 
+j 6.77 

-1.44 
+j 6.77 

-1.61 
+j 6.71 

-1.82 
+j 6.66 

-2.056 
+ 6.60 

3,4 

-0.23 
+J 6.32 

-1.34 
+j 6.32 

BOH 

EBB 

WB3M 

1.92 
+j 6.30 

5,6 

-0.17 
+j 5.32 

-1.06 
+j 5.32 

Bn 

— 

nn 

-1.60 
+j 5.47 


-0.23 
+j 4.18 

-0.97 
+j 4.18 

-1.11 

+4.20 

-1.29 
+j 4.20 

-1.518 i 

+j 4. 21 

9,10 

-0.16 
+j 3.45 

-0.68 
+J 3.45 

-0.83 

+3.39 

-1.01 

+3.32 

-1.32 
+j 3.30 

11,12 

-0.17 
+j 2.80 

-0.60 
+j 2.80 

-0.72 
+j 2.83 

-0.94 
+J 2.87 

-1.25 
+J 2.90 

13,14 

-0.19 
+j 2.60 

-0.52 
+j 2.58 

-0.65 
+j 2.60 

-0.88 
+j 2.60 

-1.18 
+j 2.60 

15,16 

-0.18 
+j 2.20 

-0.46 
+j 2.19 

-0.56 
+j 2.23 

-0.77 
+j 2.25 

-1.09 
+j 2.26 


TABLE 5.1 


Closed-Loop Poles of the Nonninal Design 
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FIG. 5.3 Closed-Loop Poles of Che Nominal Design for 
Various Values of Che ScablllCy Factor oi 
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The may<niiiiii cro88-over frequencies for the various RPDS designs are 
also computed and listed In Table 5.2. The Incremental behavior of this 
quantity with a agrees with our discussion in section 4.6. 

We shall analyze the stability robustness properties of the 
above designs In the next two sections. 

5.3 Behavior of the Closed-Loop Poles In the Face of Perturbations 
In this section, we examine the movement of the closed- loop RPDS 
poles subject to change of operating points and Introduction of un- 
modelled exciter dynamics. In the face of such perturbations, the 
ability of a state feedback, design to hold each closed-loop pole within 
a small neighborhood of Its nominal position Is a good Indication of 
Its robustness properties. The RPDS designs considered In this section 
are those obtained in section 5.2 using the simple model of the 9-machlne 
power system at operating point 1. We shall study the robustness 
properties of such designs in section 5.3.1 and the effect of a on the 
robustness properties In section 5.3.2. 

5.3.1 Robustness Properties with Respect to Change In Operating 
Points and linnpdelled Exciter Dynamics 

In section 5.2 several RPDS control laws were designed for the simple 

model of the 9-machlne power system at operating point 1. We now 

apply these designs to the other four operating points as well as the 

detailed 38 state model of the power system, and study the behavior of 

the resulting closed-loop poles. In Fig. 5.4, the closed-loop poles 

of the nominal design with choice of stability factor a equals 
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TABLE 5.2 


Crossover Frequencies for RPDS Designs 
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0,4^ are plotted together with those obtained from application of the same 
nominal control law to operating points 2 through 5. It can be seen that 
for each of the oscillatory modes, the closed- loop pole locations at 
different operating points agree very well with one another. Indeed, 
with the exception of mode 5,6 and 11,12 locations of the closed-loop 
poles corresponding to each mode are found to lie within circle of diameter 
less than 0.5 In the complex plane. 

In Fig. 5.5 the oscillatory modes of the detailed model at operating 
point 1 are plotted alongside those of the simple model for controller 
design with a ■ 0.4. It Is clear from the figure that there Is a lack of 
agreement between the closed-loop poles of the two models. With the 
exception of modes 5, 6 and 9, 10 the poles of the detailed model are In 
general less well damped than those of the simple model. Six oscillatory 
modes (modes 3,4, 11,12 and 15,16 o . the detailed model In fact 
possess damping ratio of value less than 0.2. Moreover, among these 
six modes, four of which (modes 11, 12 and 15, 16 have real part of 
their closed- loop poles less than the stability factor 0.4. 

Figure 5.6 displays the closed-loop poles obtained by application 
of the nominal control law of a « 0.4 to the detailed model at all five 
operating points. Except for the case of modes 5,6, 7,8 and 9,10, the 
closed-loop poles at different operating points for each mode are very 
close to one another. This Indicates the dominance of the unmodelled 
exciter dynamics over the change of operating points. 

^ Throughout this section and section 5.4.2 only those results of the RPDS 
design with a - 0.4 are displayed. The pattern of the closed-loop pole 
behavior observed for such design Is typical of those obtained from the 
RFDS designs with other values of a. 
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FIG. 5.4 Closed-Loop Poles Result from Application of the 
Nominal Control Law (Based on the Simple Model at 
Operating Point 1; a - 0.4) to the Simple Power 
System Model at Various Operating Points 
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FIG. 5.5 Closed-Loop Poles Related to Machine Oscillation that 
Results from Application of the Nosiinal Control Law 
(Based on the Simple Power Systsa Model at Operating 
Point 1; a«0.4) to the Simple and the Detailed Power 

System Models at Operating Point 1 
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FIG. 5.6 Closed-Loop Poles ReleCed to Hechine Oscllletlons that 
Result from Appllcetlon of the Nominal Control Law 
(Based on the Simple Model at Operating Point 1 with 
a <■ 0.4) to the Detailed Power System Model at Various 

Operating Points 
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It Is clear from the above discussion that the RFDS designs 
obtained In section 5.2 have reasonabley good tolerance of changing 
operation conditions and unmodelled exciter dynamics. This observation 
Is consistent with the excellent guaranteed stability margins derived 
for RPDS In Chapter IV. 

5.3.2 Effect of the Stability Factor a on Robustness Properties 

We next examine the effect of the stability factor a on the 
behavior of the closed- loop poles subject to perturbation of the system 
dynamics. This Is carried out by applying the RFDS deslgxis with different 
values of a to operating points 1 through 5. Recall from the last section 
that the magnitude of the real part of all the closed-loop oscillatory 
inodes display an Incremental behavior with a for values of a between 0.0 
and 0.6. Moreover, such Increment Is uniform In the sense that If a 
Is Increased by an amount equals to Aa , then every closed-loop complex 
pair will be shifted horizontally to the left by a distance roughly 
equal to Aa. 

Figure 5.7 displays the closed- loop poles that result from applying 
the RFDS designs of Section 5.2 (which are based on the simple 9-machlne 
model at operating point 1) to operating point 1 for values of a 
equals to 0.0, 0.2, 0.4 and 0.6. It Is observed that as we Increase 
the value of a , all the osclllary modes move further out to the left 
In the complex plane. However, unlike the case of the nominal design , 
the amount of pole shifting observed here Is no longer proportional to 
the Increment of oi. It Is also observed that the shift with respect 
to Increment of oi tends to be larger for higher values of a. 
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FIG. 5.7 Closed-Loop Poles Result from Application of the 
Nominal Control Laws (Based on the Simple Model 
at Operating Point 1; a » 0.0, 0.2, 0.4, 0.6) to 
the Simple Power System Model at Operating Point 1 
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FIG. 5.8 Closed-Loop Poles Related Co Machine Oscillations Result 
from Application of the Nominal Control Laws (Based on 
Che Sliqile Model at Operating Point 1; a « 0.0, 0.2, 0.4, 
0.6) to Che Detailed Power System Model at Operating Point 1 
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The effect of a on the ability of RFDS to tolerate unmodelled 
exciter dynamics Is examined next. The closed-loop poles corresponding 
to machine oscillations of the detailed model at operating point 1 are 
plotted In Figure 5.8 for state feedback designs with value of a equals to 
0.0, 0.2, 0.4 and 0.6. 

It is clear from the figure that the resulting pole pattern again 
displays an Incremental behavior with <x. In that the closed- loop poles 
associated with BFDS design using larger value of a are hold further 
back from the ju-axls in the face of unmodelled exciter d 3 mamlcs. 

The above observations suggest Improvement of RPDS ability to 
maintain stability In the face of changing operating points and un- 
modelled excltor dynamics with Increasing value of a . This Is a 
consequence of the fact that large value of a will result In the 
closed-loop poles being positioned further away from the J&i-axls. It 
does not however suggest that the ability of RPDS to hold the closed-loop 
poles near their nominal position under perturbations Improves with 
Increasing value of a . 

5.4 Frequency Domain Robustness Analysis 

Based on the unity negatlye feedback representation of RFDS in 
Figure 4.7, frequency domain robustness characterization for the state 
feedback controllers derived In section 5.2 can be readily obtained. 

All the RFDS designs studied In this section are again based on the simple 
model of the 9-machlne power system at operating point 1. The behavior 
of OmlnCl. + and amin(i + T^^(jw)) with changing ct are examined 

In section 5.4.1. Robustness analysis of the type described In Theorems 4.4 
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and 4.5 are given in section 5.4.2. 


5.4.1 Effect of the Stability Factor a on the Return Difference and 
Inverse Return Difference Matrices 


Minimum singular values of the return difference matrices are plotted 
in Figure 5.9 for the RPDS designs obtained in Section 5.2. All the 
singular value plots are found to display a peak in the vicinity of 
0.5 Hz, and thereafter roll off in a first order fashion. Moreover* it 
is clear from the plots that the quantity amin(^ ^ (Jco)) increases 
monotonically with a for all values of (a . This latter observation is 
consistent with our conclusions in Section 4.4.4 regarding the 
behavior of RPDS return difference matrices with increasing 

The plots for the complementary quantity amin(I_+ ^ of 

the same RPDS designs are displayed in Figure 5.10. In view of the 
peak near 0.5 Hz observed in the plots for a min (I^ it is not 

at all surprising that the plots for Omin(I + T ^(jw)) should display 
a valley at about the same frequency. Beyond this frequency* the plot 
for amin(^ + ^ began to rise proportionately with at in a 

first order fashion. It can also be seen from the plots in Figure 5.10 
that amin(_I + is a decreasing function of a for each value of w. 

Consequently, the guaranteed margins (as given in (4.12) and (4.13)} 
for such design will deteriorate with increasing value of a. 

1 

Due to the presence of the matrix ^ In the algebraic Rlccati equation 
(5.3), one may question the applicability of the conclusions in Chapter IV 
to the present design example. A closer look at the problem however 
dismisses such a suspicion. Using the construction of the matrix ^ 
and the property of^the algebraic Rlccati equation, it can be readily 
demonstrated that ^ S, ■■ 0 where ^ is the unique positive definite 
solution of (5.3). This in turn reduces the Rlccati equation (5.3) 
to one identical to (4.39). 
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FIG. 5.9 Plocs of amia(X + T (joi)) for Nominal RFDS Designs 

Based on the Simple^ Power System Model at Operating 
Point 1 with a equals to 0.0, 0.2, 0.4 and 0.6 
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FIG. 5.10 Plots of amln(^ + ]^^(j(o)) for Nominal RPDS Designs 
Based on the Simple Power System Model at Operating 
Point 1 with a equals to O.Ot 0.2, 0.4 and 0.6 
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5.4.2 Robustness Tests Based on Singular Values 

In order to apply the robustness tests described In Theorems 4.4 
and 4.5 to the present example, the multiplicative error matrix L(s) with 
respect to the nominal system model are computed for different perturbed 
system dynamics. The nominal system model consists of the simple model 
for the open-loop power system at operating point 1 which Is regulated by 
a state feedback controller with prescribed degree of stability equal 
to 0.4. This controller Is taken directly from the respective RFDS 
design In Section 5.2. 

The various singular values specified by conditions (4.22)and(4.27) 
are plotted In Figures 5.11 and 5.12 respectively for perturbations 
corresponding to change of operating points. The singular values plots 
for CTmax(^"^(jaj) - p and amax(L(jcj) - ar found to display drastically 
different behavior for different operating points. In particular, the 
values of omax(^ ^(Jw) - H and (7max(Mj<a) - ^ result from switching 
from operating point 1 to operating point 2 are noted for their exceptionally 
large magnitude at low frequencies. This pattern of behavior Is In 
contrast with the Insignificant change of the closed-loop pole positions 
observed In the last section for precisely the same class of perturbations 
(l.e. change of operating points). It becomes clear from this observation 
that robustness characcerlzatlon based on matrix norms can be very 
conservative. Conservatism of the norm-based robustness test Is further 
manifested by the fact that conditions (4.22) of Theorem 4.4 and (4.27) 
of Theorem 4.5 are violated for some nondestablllzlng perturbations 
(such as the change of operating condition from operating point 1 to 
operating point 2) . 
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FIG. S.ll Singular Value Plots for the Robustness Test given In 

Theorea 4.4. The Nominal System is Obtained by Applying 
the RPDS Design in Section S.2 (That are Based on the 
Simple Power System Model at Operating Point 1; a ■ 0.4) 
to the Simple Model of the Open-Loop System at Operating 
Point 1. The Perturbations are due to Change of Operating 
Points 
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FIG. 5.12 Singular Value Plots for the Robustness Tests Given in 
Theorem 4.5. Nominal Design is Obtained by Applying the 
RFDS Design in Section 5.2 (That is Based on the Simple 
Power System Model at Operating Point 1; a ■ 0.4) to the 
Simple Model of the Open-Loop System at Operating Point 1. 
The Perturbations are due to Change of Operating Points 
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FIG. S.13 Singular Value Plots for the Robustness Test Given In 

Theorem 4.4. Tne Nominal System Is Obtained by Applying 
the RPDS Designs In Section 5.2 (That are Based on the 
Slosple Power System Model at Operating Point 1; a ■ 0.4) 
to the Simple Model of the Open-Loop System at Operating 
Point 1. The Perturbations are due to both Change of 
Operating Points and Unmodelled Exciter Dynamics 


-163- 


/ 

/ 

/ 

/ 



FIG. 5.14 Singular Value Plots for the Robustness Tests Given 

in Theorem 4.5. Nominal Design Is Obtained by Applying 
the RPDS Design In Section 5.2 (That are Based on the 
Simple Power System Model at Operating Point 1: a ■ 0.4) 
to the Simple Model of the Open-Loop System at Operating 
Point 1. The Perturbations are due to Both Changes of 
Operating Points and Unmodelled Exciter Dynamics 
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The «bov« obflcrvaclom also apply to parturbatlons that rasulc 
from Introduccloo of unmodellad exciter dynamics in addition to change 
of operating points (Fig. 5.13 and 5.14). The respective value of 
Offlax(L(j(o) - 2) *nd amax(L ^(Jco) - p at each value of cu are found to 
be larger than the corresponding quantities result from changing of 
operating conditions alone (compare Fig. 5.11 and Fig. 5.12 v^th 
Fig. 5.13 and Fig. 5.14 respectively). This again indicates the 
dominance of the unmodelled dynamics over the change of operating point. 

5.5 Coi:cluding Remarks 

We have studied In this chapter the application of RPDS technique 
to design of state feedback control laws for a multl-censlnal DC/AC 
power system. Our major objective is to demonstrate with the aid of a 
nontrivial multivariable design example the various properties of RPDS 
discussed in the previous chapter. 

With regard to the positioning of closed- loop poles t the RPDS 
design obtained in section 5.2 are found to possess an interesting 
property which is not predicted by the results developed in this thesis. 

For values of a between 0.0 and 0.6* a given increment of a will shift the 
doted - loop oscillatory modes horizontally to the left by approximately 
the same amount. This unexpected property of the RPDS design is probably 
a result of the fact that all the open- loop modes to be stabilized are 
of the same nature (i.e. they all correspond to modes of intexmachine 
oscillations). It is by no means a result that applies in general. 

Robustness properties of the RPDS designs In section 5.2 are 
evaluated in two complementary ways. In section S*3t movement of the 
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closed-loop polM thac result from change of operating points and Intro- 
duction of unnodelled e:rclter dynamics Is axaalned. In all cases under 
consideration, the RPDS state feedback controllers are found to maintain 
a reasonable degree of damping for all the closed- loop oscillatory modes . 
Moreover, the closed-loop poles associated with RFDS designs that employ 
a larger value of a are always held further sway from the Jot-axis In the 
face of perturbation. These observations are consistent with the excellent 
stability margins derived for SPDS systems In Chapter IV. In section 5.4, 
the frequency domain robustness tests prescribed by Theorem 4.4 and 4.5 are 
applied to the RPDS designs under consideration. Perturbation of system 
dynamics that produce closed-loop patterns similar to one another are 
found to display drastically different behavior In the plots for their 
respective cmax(L ^(Ju) - I) and amax(L(Ja)} - p . This Is an Indication 
of the conservatism associated with the norm-based robustness tests. 
Conservatism of such, tests is further reflected by the fact that condition 
(4.22)of Theorem 4.4 and (4.27) of Theorem 4.5 ere both violated In case 
of nondestabilizing perturbation. We also examine the behavior of the 
two MIMO frequency domain robtistness measures OmLa(.l + T^(J<>i)) and 
amln(X ^ Cju)) with respect to change In a. They are shown to be 
monotonlcally Increasing and decreasing functions of a respectively for 
all values of oi . The result obtained for the quantity omln(^ -f ^(Ju)) 
agrees with our conclusion In section 4.4.3 while that obtained for 
Cmln(^ ^^(Jo))) confirms a conjecture we made in section 4.5.3 concerning 
the property of KPDS Inverse return difference matrix. 
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CHAPTER VI 


DESKai OF RZGDLATOR WITH A PRESCRIBED DEGREE 
OF STABILITY BASED LOG COMPENSATORS 


6.1 Introduction 

A basic practical limitation associated with the RPDS design Is 
the assumption of full state feedback. In many practical applications 
full state feedback can never be exactly realized and often It Is either 
impossible or too expensive to provide enough sensors for achieving even 
an approximate realization. 

The way this problem is handled In modem control theory Is 
through the use of LQ6 methodology [Ath 1] , In which a Ralman-Bucy filter 
Is used to provide the necessary state estimates using noisy output measure- 
ments. The class of LQG controllers considered here are the RPDS based LQG 
controllers. The state feedback gains for such controllers are obtained 
using the RPDS design methodology discussed In Chapter II. 

The standard configuration for a RPDS based LQG control system 
Is depicted In Fig. 6.1 with various points of the loop marxed. To 
determine the robustness properties of the design, we shall Insert miltl- 
pllcative perturbation of the type considered in Section 4.2 at points 
(2) and (3) and find out the tolerable magnitude of the model error thac 
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wlll not destabilize the system. The other two points (1) and (4) are 
Internal to the LQG compensator and are therefore less significant for 
robustness analysis. However, they have some desirable loop properties 
that can be related to the other two points. 

Following the notation employed In [Le Ij , we shall denote the 
loop transfer matrix at point (K) by T„.(s). Each ^Cs) Is calculated 
by breaking the loop at (K) and treating this point as both the Input 
and output. For the four points marked In Fig. 6.1, we have 


Tj^(s) - ^(s) B (6.1) 

12 (a) « (£"^(s) + BG^ + HC)"^HC £(s) B (6.2) 

T,(s) - C <<>(s) B G (<^"^(s) + B G + H C)“^ H (6.3) 

T^C4) « C Ks) H (6.4) 

where ^(s) ^ (s^ - A)'^ (6.5) 


^ is the RFDS state feedback gain and H Is the KBF gain. 

Note that points (1) and (4) have loop transfer matrix Identical to those 

of RPDS and KBF respectively. Thus at point (1) Che RFDS robustness 

properties apply. Similarly, the KBF robustness properties are valid at 

point (4) . No guaranteed stability margin Is however available at both 

points (2) and (3), which are the actual Interface between the control 

system and the real world. It was demonstrated by Doyle [Do 4] with a 

3lmple example that a reasonable looking LQG design may have arbitrarily 

^It Is important to point out chat only Che stace feedback gain of Che 
RPDS based LQG compensaCors buC not Che KBF gain Is designed wlch a 
prescribed degree of scablllcy. Robuscness propercles of KBF are 
Identical to those of FFDS discussed In section 4.8.2. 
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soall stability margins. 

For conventional LQ6 design (one for which the state feedback 
gain is obtained by using the LQ instead of the RFDS method) .there are 
two dtial robustness recovery procedures that allow us to approximate 
the transfer matrices Tj^(s). and ^(s) with T^Cs) and respectively 

in a systematic fashion. These procedures use the asymptotic pole pro- 
perties of LQ regulators and KBF respectively (see [Do 3] for a review 
of such properties ) , and can be applied only if the plant is minimum 
phase. 

The robustness recovery procedure due to Doyle and Stein [Do 1] 

makes T_ 2 ^®^ approximate T^(s) by using a process noise with spectral 

T - 

intensity of the form p ^ B + c in the KBF Rlccati equation and letting 
p go to infinity. As a result, the LQ robustness margins at point (1) 
can be recovered asymptotically at the input (2) . In a dual fashion, 
the robustness recovery procedure due to Kwakemaak [Kw 3] is used to 
recover the stability margin for (4) at the output (3) . This is 

T 

accomplished by using a state weighting matrix of the form ^.'i’P^ £ 
and letting p go to infinity. 

The objective of this chapter is to study the various Issues 
that arise in the application of robustness recovery procedures to 
design of RFDS based LQG compensators. 

6.2 Design of RFDS Based LQG Compensators Using Robustness Recovery 
Methods 

In this section, we illustrate with the aid of numerical examples 
some considerations that are of importance to design of RFDS based LQG 
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compensators • 

6.2.1 Relationship Between the Stability Factor a and the Minimum 
Phase Condition 

If the stability margins associated with the loop transfer matrix 

T^(s) of a given RFDS based LQG compensator (with prescribed degree of 

stability a) Is found to be unsatisfactory, an obvious way to Improve the 

design Is to employ the robustness recovery procedure of Kualcemaak 

([IStf 1] and [Kw 3]). Recall from the last section that Ruakemaak's method 

requires the adjustment of the state feedback gain G_^ for making T^(s) 

appropriate the KBF transfer matrix T^(s). This Is accomplished through 

T 

the Introduction of a state weighting matrix of the form 2. P£ £ 

Into the LQ Rlccatl equation and letting p go to Infinity. The finite 
closed- loop regulator poles will then asymptotically approach the zeroes 
of £(s^ “ A) ^ Theorem 4.13 of [Ru 1]). However, the LQ regulator 

thus obtained Is no longer guaranteed to possess the same prescribed 
degree of stability as the original RPDS design. This Is the case when 
some zeroes of C(s£ - B have real parts larger than - a. 

If the RPDS Rlccatl equation (4.59) Is used as the design 
equation. It follows from Theorem 3.2 that the resulting regulator poles 
will always lie to the left of <7 - -a for every value of P . But the 
asymptotic location for these finite poles may not coincide with those 
of C{al - B unless C(s£ - A - al)”^ B Is also minimum phase. If 

- a - B falls to be minimum phase, the asymptotic state 

feedback gain will satisfy 
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V G -*■ W C C6.6) 

/p -a - 

as p approaches Infinity where £ C and W Is some orthononaal matrix 
(see Chapter 3 of iKw 1] ) . It can be readily demonstrated with the use 
of elementary matrix manipulations that full state KBF loop cannot be 
recovered asymptotically as a result of £ ^ £. 

Defining 

?(s) - (si - A - H C)"^ (6.7) 

we can rewrite T^(s) as 

13(3) ■ £ i(s) i ^ (£"^(s) + B ^ + H £)“^ H 

- ££(s) B ^[£(s) - |;(s) B (I + ^ f(s) B)’^ £(s)] H 

-r -1 

- C Ms) B (I + ^J(s) B) ^ Ms) H (6,9) 

To go from (6.8) to (6.9), we employ the matrix inversion lemma 
(see for example Appendix A of [Sc 1]). As p approaches infinity the 
following approximation follows lomediately from (6.6) and (6.19) 

C (J)(8) B (I + G Ms) B)"^ G Ms) H 

— — — -Or- — -a — “ 

-*• £M») B (£Ms) B)"^ £Ms) H (6.10) 

It can be shown that the matrix function on the right side of (6.10) Is 
equal to the KBF loop transfer matrix £ ^(s) H when £■ £ • If £ )* £ , 
there exists no useful simplification of this matrix function and the 
robustness properties of the KBF loop is not recovered as a result. 

The following example Illustrates the effect of a on the design 
of LQG compensators using Kwakemaak's robustness recovery procedure. 






Example 6.1 Consider the LQG problem 


min J 

u(.t) e f (y(t)) 



Ix^Ct) £i(t) +u^(t> R u(t)Jdt 

C6.ll) 


subject to the dynamic constraint 


£(t) 


r 0 ll x(t) + r O’! u(t) + r 351 C (t) 

L-3 J LlJ L-61. 


and the observation constraint 


( 6 . 12 ) 


y(t) - 12 1] x(t) +e(t) 

(6.13) 

where _ 

g • 80 [ /35 1/35 1 ] 

L 1. 

(6.14) 


and C (t) * 3 (t) are zero mean white noise processes of spectral Intensity 
equal to 1.0. 


The plant In this example la a stable minimum phase system with 
transfer function given by 


zSsl 

u(s) 


8 2 

(s+1) (S+3) 


(6.15) 


Solution of the above LQG problem for a » 0 results in a state feedback 
gain given by 

G - 150 10] (6.16) 


and a Kalman Bucy Filter gain given by 
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Inspectlng the Nyqulst diagram of the LQG loop transfer function reveals 
poor stability margin ((3M ■ (r6.75 d,b «) and PM ■ C-15®, 15^1) for the 
resulting LQG design. This Is due to the presence of a Instable pole In 
the LQG loop transfer function. In order to Improve the robustness 
properties of the feedback loop, we employ the robustness recovery pro- 
cedure of Kwakemaak. A new choice of the state weighting matrix given 
by 

Sp -a+ppj [2 1] (6.18) 

Is employed In the LQ algebraic Rlccatl equation 


K 

0 1 

+ 

~0 

-4* 

K - K 

• 4 

0 

[0 

1] K + - 0 


-4 -3 

• « 


1 

-3 


1 

• m 


(6.19) 


The state feedback gains for various values of p are computed and the 
respective LQG loop transfer functions t^(s) are plotted In Fig. 6.2. 

It Is observed that the stability margins of the resulting designs Improve 
steadily with Increasing p (see Table 6.1). For a > 0 we replace 
equation 6.19 with the following RFDS algebraic Rlccatl equation 

[0 IJ + ^ • 0 (6.20) 

the resulting LQG design Is always guaranteed to possess a degree of stabi- 
lity a. However, asymptotic recovery of the KEF loop with this modified 
procedure Is only possible for values of a less than 2.0. This Is vividly 
demonstrated by the Nyqulst diagrams depicted In Fig. 6.3 and Fig, 6.4. 

When a Is chosen to be 1.8 (Fig. 6.3), the stability margins of the 


a 1 


a -4 


0 


+ 




-4 -3 + 0 

M m 


1 -3+X 


1 
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FIG. 6.2 Nyqulst Diagrams for Design Iteration of RPDS 

Based LQG Compensators (a ■ 0} Using Kwakemaak's 
Recovery Procedure 
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p 

PM^ 


d" 

0 

15.3® 

-2.0db 

0.29 

100 

21.8® 

-2.7db 

0.39 

500 

35° 

-5.51db 

0.58 

1000 

45.4® 

-9.6db 

0.9 

10000 

72° 

4 

-34db 

0.917 


TABLE 6.1 


Suomary of Stability Margin for Design Iteration of 
a SFDS Based LQG Compensator (a« 4.5) In Example 6.1 


^ This Is In violation of the notation Introduced In Chapter IV. 

The angle given here corresponds to the value of 6 In the 
definition of PH. 

2 

Only the downward gain margin Is given here. The upward gain 
margin Is Infinity for all cases. 

3 

d Is the nearest distance between the MyqulsC diagram and the 
critical point (rl»0) . 

The actual gain reduction margin Is lower than the value indicated 



-176- 



FIG. 6. 3 Myqulsc Diagrams for Design Iteration of RPDS 

Based LQG Compensator (a *1.8} Using Kwakemaak's 
Recovery , Procedure 
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FIG. 6.4 Nyqulst Diagrams for Design Interaction of a 
RFDS Based LQG Compensator (a ■ 2.2} Using 
Kwakemaak's Recovery Procedure 



/I <10000 


FIG. 6.5 Myqulsc Diagrans for Design Tteraclon of i 
SroS Based LQG Compensator (a “ 4.5) Using 
Doyle/Stein's Robustness Recovery Hethod 
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design Improve steadily with Increasing value of p and reasonably good 
stability margins are obtained for valrs of p equals to 1,000, For 
case where a equals 2,2 (Fig. 6,4} the stability margins of the resulting 
design remains unsstlsf story even when large value of p Is employed. 

The above problem Is not shared by the robustness recovery pro~ 
cedure due to Doyle and Stein, Nyqulst diagrams of the transfer functions 
that result from the design iterations of a SPDS based LQG compensator 
using Doyle/Steln's method are depicted In Fig. 6.5. The choice of ot 
chosen for this cose Is 4.5. It Is clear from the figure that the stability 
margins of the design Improve as larger value of p Is used demonstrates 
satisfactory improvement. 

6.2.2 Rate of Robustness Recovery with Respect to the Noise Scaling 

Parameter 

The robustness recovery method of Doyle and Stein [Do 1] 
requires the use of a process noise with spectral Intensity of the 
form 

_ T 

1 “ £ + P i B (6.21) 

and letting p go to Infinity for asymptotic recovery of the fttU state 
RFDS stability margins. This procedure provides characterization of the 
LQG feedback loop as the noise scaling perasieter p approaches Infinity. 

It does not, however, give us any clue as to the behavior of the LQG loop 
when p varies. We shall ezomlne in this chapter the effect of a on the 
rate of recovering the RFDS stability margins with respect to changes In p. 
In the SISO case, this can be simply accomplished by Inspection of the 
respective Nyqulst diagrams. In the MHfO case, the singular value 



plots of the respective loop transfer matrices have to be used (For an 
example of this, see [So 1]). For ease of exposition^ only SISO systems 
will be considered In the following exai^le. 

ft. 2 Consider the RFDS based LQG design problem defined In Example 

6.1. Design Iterations with noise scaling factor P ■ 0, 100, 500, 1000 and 
10000 are performed for several choices of the stability factor . The 
Nyqulst diagrams of the LQG design with 0l chosen tobeO, 3, 4.5, 6 and 9 
are plotted for each of the five Iterations (Fig. 6.6) and the resulting 
stability margins tabulated (Table 6.2 and Fig. 6.7). Several Interesting 
observations are In order. 

First It Is noticed that the Nyqulst diagrams of the LQG designs 
with value of a equals to 0, 3 and 4.5 move further away from the critical 
point (-1,0) than those designs with a equals to 6 and 9 as the value of p 
Increases. Consequently, the stability margins of the designs with a 
equals to 0,3 and 4.5 are superior to those with values of <x equals to 
6 and 9 for a fixed value of p (see Table 6.2 and Fig. 6.7) . Secondly, 
the rate of robustness recovery with respect to p does not strictly 
decrease with a for the values of a under consideration. It Is clear 
from the Nyqulst diagrams In Fig. 6.6 and the data presented In Table 6.2 
that the stability margins for LQG designs with a eqtials to 4.5, 6.0 and 9.0 
decreases with Increasing value of a for a fixed value of p. This is 
however not the case for designs with a equals to 0 and 3. 

In view of the above observation. It does not seem likely that one 
can make any conclusive statement concerning the effect of a on the rate of 
robustness recovery. However, these observations do Illuminate a potential 



( 'a»6.0 

'^a»0.0 

P 

Fig 6.6b 

Nyqulst Diagrams for Design Iterations of RFDS Based 
LQG Compensators (a « 0.0, 3.0, 4.5, 6.0 and 0.0} 
Using Doyle/Steln Method 
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Slnce the Nyqulst Diagram of the RFDS Transfer Functions 
Corresponding to a * 0.0, 0.3 and 0.45 are very closed 
to one another, only that of a ■ 0.0 Is Displayed 
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Table 6.2 to be continued 
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TABLE 6.2 

Summary of Stability Margins for Example 6.2 


^ This is a violation of the notation Introduced in Chapter IV. The angle 
given here corresponds to the value of 8 in the definition of FH 

2 

Only the dotnarard gain margin is given here. The upward gain margin 
is infinity for all cases. 

^ d is the nearest distance betvreen the Myqulst diagram and the 
critical point (-1»0). The actual gain reduction margin is lower 
the value indicated. 

4 

The actual gain reduction margin is lower than that Indicated. 
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problem that may occur In the design of 8PDS based LQG compensator using 
robustness recovery procedures. It is apparent from this example that 
large values of a (p( ■ 6,9) may slow down the recovery of RPDS stability 

\ 

margins and thus necessitate the introduction of high filter gains for 
securing good robustness peroperties. The use of high filter gain in 
feedback loop often leads to degradation of observation noise rejection 
that is undesirable from the performance point of view. The resulting 
situation is one where a satisfactory compromise among the need for 
speed of response (as specified by the choice of a) , noise rejection and 
stability robustness has to be made by appropriate choices of a and p. 

6.2.3 Effect of the Stability Factor a on the Noise Relection Properties 

To conclude the dlscxtsslon in this section, we examine the effect 
of the stability factor a on the noise rejection properties of RPDS 
based LQG compensaturs . Such effect can be characterized in terms of the 
behavior of the state covariance matrix (assuming zero mean for the random 
vector x(t) 

1 - E[5(t) J(t)3 (6.22) 

— * 

with different values of a. 

Recall £.rom the definition of the state estimate error £(t) in 
(4.171) that 

£(t) ■ x(t) - x(t) (6.23) 

It is useful to point out that x(t) and e,(t) are uncorrelated 
random vectors (Chapter III of [An 2] and Chapter IV of [Kw 1]). As a result 
of thlst we can express I as a sum of two covariance matrices 
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I , • Z ; + ^ (.6.24) 

where 

^ • Elx(t) x(.t)^] (6.25) 

and 

Z - E[e(t) e(t)^] (6.26) 


Moreover, the dynamics of the LQ6 control system can be conveniently 
described by using [si(c) x(t) ] as a state vector. The underlying 
differential equation in this case is given by 


e(t) 


A - H C 0 


e(t) 


-1 H 

r 

A 

■ 



A 



L”. 

x(t) 

^ m 


-H C A - B G(j 

m • 


x(t) 

m <m 


1 

io 

f— 

|<E 

ft 

1 


(6.27) 

which can be readily derived by inspection of the block diagram in 
Fig. 6.1. Using the above d 3 rnamlcal equation, and Theorem 1.53 of 
(Kw 1] (which is a standard result in stochastic process theory) leads 
directly to the following restilt. 


Learns 6.1 Given a LQG control system described by (6.27) the steady 

state covariance matrix 


E 


( 


•(t)' 

x(t) 




(6.28) 


satisfies the algebraic matrix equation 
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(6.29) 


It is strslghtfoxward to show that ths solution of equation 6.29 can be 
obtained by solving the two following decoupled matrix equations 


and 


(A-HC) Z + Z(A-HC)^+ (£+ H 9 H^) -0 


<A - B ^ + ^(A - B + I 9 - 0 


(6.30) 


(6.31) 


This is a direct consequtmce of the fact that x(t) and £(t) are uncorrelated. 
It beconap clear from (6.30) and (6.31) that the RPDS state feedback gain 
(and hence the stability factor a) only affects the matrix Z^. Differen- 
tiating both sides of (6.31) with respect to a and rearranging* we obtain 


where 


* 3T<i-iV - £ 




(6.32) 


(6.33) 


The matrix 


3 Z'' 
” tx 

3a 


allows us to access the effect of stability factor 


a-a 


on the noise covariance matrix Z ^ nt a ■ a * 


If the ith diagonal 
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3£. 

clement of la poaitlye, then the variance of the 1th atate will 

Increase with a at a « a , Likewise, If the 1th element of 


3Z'' 

— X 

3a 


Is negative, Chen the variance of the 1th state will decrease 


a-o 


with a, at a*V.Slnce(^ “ 1. £ ) Is a stable matrix, we can thus conclude 

from the well known properties of algebraic Lyapunov equation that 

3r 

-X ' 

“3a 


_ Is positive (negative) semldeflnltc (l.e the variance of all 
a-a 

states Increases (decreases) with a at a ■ a). If the matrix D Is 
negative (positive) semldeflnlte. If D turns out to be indefinite. 


the resulting 

3a 


aau 


Is also Indefinite. In this case, the 

Increase of will Increase the variance for some states and decrease 
that for others. 

The following example Is again based on LQG problem considered 
In Exaiq>le 6.1. It demonstrates the behavior of ^ with Increment 
of a . 


Example 6.3 Consider the RFDS based LQG controller design 

problem in Example 6.1. The noise covariance matrlx£^ Is cosiputed for 
the resulting designs at five different values of a and p Is chosen 
to be 0 (see Table 6.3). It Is clear from the tabulated results that 
the variance of the first state decreases with awhile that of the 
second state Increases trlth a . 

6.3 Concluding. Remarks 

We have Identified in this chapter several potential problems 
that may result from the use of robustness recovery procedure in design. 
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of RPDS based LQG componsator. Too large a choice of the stability factor 
Is found to prevent the recovery of the full state KBF stability 
margins with the Kwakemaak's procedure and slow down the rate of 
robustness recovery with the Doyle/Steln's procedure. However, It Is 
Important to point out that the result on rate of robustness recovery 
Is obtained only for a particular single-input system. More practical 
experience with application of robustness recovery method to MSfi) RPDS 
based LQG compensator design Is needed before the nature of this problem 
can be fully comprehended. 
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j 

i 

CHAPTER VII 

SHMHARY. CONCLUSIONS AND SUGGESTED DIRECTIONS FOR 
FUTURE RESEARCH 

7.1 Suaaary | 

We have examined in this thesis a wide range of problems related to 
RPDS methodology and Its applications. This Includes attesipts to 

(1) explore the use of RPDS methodology for time-varying systems 

(11) adapt RPDS methodology for LQ regulator problem with design 
specifications other than prescribed degree of stability 

(111) develop methods of elgenstructure analysis for RPDS control 
systems 

(iv) clarify the robustness properties of RPDS In the multiple- 
input case 

(v) Identify potential problems that may occur in the design of 
RPDS based LQG compensators. 

We began our Investigation in Chapter II with a formulation of the 
RPDS problem for linear time-varying systems. A generalized notion of 
’degree of stability' that applies to all flnlte-dlmenslonal linear 
systems Is introduced. This definition has the desirable property that 
It reduces to characterization In terms of eigenvalues for the LTI systems. 

It turns out that the exponential weighting technique for solving the 
time- Invariant RPDS problem Is equally applicable to the time-varying 
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case. The related problem of designing Kalman Bucy filters with prescribed 
degree of stability (FPDS) Is fomnilated and solved In a dual fashion. 

For case of LTI systems with stationary noises, the optimal filter solution 
admits Interesting Interpretations In terms of the special form of noise 
Intensity matrices required for speeding up the error dynamics. 

In Chapter III, several eigenstructixe characterizations of the 
time- Invariant RFDS are derived. The sensitivity equations for the 
SFDS poles are obtained In two different ways. The first approach Is 
a direct application of the classical eigenvalue sensitivity result. 

The second approach utilizes the special elgenstructure properties 
of Hamiltonian System associated with RPDS problem. The computational 
requirements for these methods are briefly compared. Asirmptotic 
behavior of the RPDS root-loci is studied next. It Is shown that the 
properties of the RPDS root-loci can be readily derived from the optimal 
root-loci properties of a related LQ regulator problem. Based upon the 
behavior of RPDS poles as the control weighting on the states becomes 
vanishingly small, a novel algorithm for designing regulators with 
prescribed damping ratio (RPDR) Is developed. 

The important design issue of robustness Is considered In Chapter IV. 
Based on the framework of frequency domain robustness analysis due to 
Lehtomakl [Le IJ , various robustness properties of RPDS are characterized 
In terms of the minimum singular value of the RPDS return difference and 
Inverse return difference matrices. In particular, the RPDS designs 
with R chosen to be diagonal are found to possess excellent gain and 
phase margins with respect to the stability and degree of stability 
property. However, tolerance of uncertainties for RPDS will Improve 
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wlch increasing value of atabllity factor only for specific types of 
aodel error representation such as E^(s) ■ (T ^(s) “ T(s)) X 

Chapter V continues the discussion of Chapter IV with a milti- 
termlnal DC/AC power systaa example. SFDS state feedback design for 
several choices of the stability factor a are applied to a 9«inachlne, 

4 terminal DC/AC power network. Robustness properties of the closed- 
loop system thus obtained are evaluated in two ways. First, the closed- 
loop pole pattern results from change of operating points and Introduction 
of unmodeUed dynamics are studied. Second, the minimum singular value 
of the return difference and the inverse return difference matrices are 
computed and compared with the magnitude of the respective type of model 
errors as specified in the robustness tests of Chapter IV. 

Output feedback realization of RFDS using LQG methods Is considered 
In Chapter VI. In view of the lack of gtiaranteed stability margins for 
such compensators, only those RFDS based LQG control systems designed with 
the robustness recovery procedures are considered. Particularly, we examine 
the effect of a on the recovery of stability margins for full state 
feedback loop. It is found that too large a value of ^ may prevent 
the recovery of the KBF loop stability margin using Kwaakemaak's method 
and slow down the rate of recovering the RFDS stability margin using 
Doyle/S teln's method. 

7.2 Conclusions 

The major contributions of this thesis are basically of two 
categories. In the first category are results related to the extension 
of RFDS methodology. The classical RFDS problem formulation and Its 
solution technique (due to Anderson and Moore) Is found to be useful 
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£or solving a number of RFDS related problems not previously considered 
in literature. These Include 

(1) tha extension of exponential weighting technique to 
solution of the time-varying RPDS problem 
(11) the formulation of RFDR problem as a special case of 
RPDS problem 

(ill) the adaptation of optimal root-loci results to RPDS root-loci 

(iv) the adaptation of LQ eigenvalue sensitivity results to 

the respective problem of RPDS 

While most of these results are of theoretical Interest, those obtained 
for RPDS root-loci are also of Importance from the design point of view. 

The contributions in the second category are mainly related to design 
implications of RPDS. While a large value of a can lead to regulators 
with good damping properties, other design considerations will put an 
upper limit on the actual value of ^ to be used. Some of such design 
considerations discussed in this thesis are 

(i) effect of a on the cross-over frequency - too large a value 
of a may extend the cross-over frequency of RPDS well into 
regions dominated by unknown and/or unmodelled dynamics 
(11) effect of a on robustness properties - Increasing the value 
of a will lead to Improvement in tolerance of modelling error 
only in very specific context. In other words, such improve- 
ment is only valid for certain type of model error representations. 
For Instance, we have demonstrated that the guaranteed stability 
margins of a RPDS may in fact deteriorate with Increasing value 


of a. 
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(111) effect of the stability factor a on robustness recovery 
of LQG - too large a choice of oi may lead to failure In 
recovering the loop with Kwakemaak's procedure and 
Impede the rate of recovery of BFDS stability margin with 
Ooyle/Sceln's procedure. 

To summarize, kPDS method Is not merely a procedure to be used 
blindly for design of fast response sytems. It Is clear from our 
discussions In Chapters IV, V and VI that a fair amount of Iteration on 
the design of RFDS, with due regard given to various design considerations 
such as stability robustness and noise rejection Is necessary to obtain 
satisfactory results. 

Like the LQ regulator method of which It Is a special case, the 
RPDS design procedure Is basically a multi- loop procedure. With the 
aid of the various recently developed frequency domain tools (such as the 
singular value-based robustness tests for Mllfi) systems and the robustness 
recovery procedures). It should provide a reasonable starting place to 
design feedback systems with a prescribed degree of stability. 

7.3 Suggestions for Future Research 

Due to the lack of time, we were not able to pursue In depth the 
many Interesting avenues of research opened up by this work. As 
suggestion for future research, we list the following: 

(1) The effect of the stability factor on the tolerance of 
structured model errors 

The robustness tests considered In this thesis use only the magni- 
tude Information of the error. As a result, they can be unduly conservative 
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In that some of the small pertubatlons that will theoretlced .17 destabilize 
the system will not occur physically. These teats can be further refined 
to take into account the difference between model errors that Increase 
the stability margins of the feedback systems and those that decrease It 
(see [Le 1] Chapter 4). A useful avenue of research Is to study the 
effect of a on the robustness properties of RFDS when structural Infor- 
mation of the error &(s) (l.e. numerical relations among elements of ^(s)) 
Is taken Into account. 

(ID Discrete Time RFDS 

The formulation and solution of discrete time RFDS problems have 
been considered by several authors ([Sa 2], [An 2J). In the LTI case, 
the corresponds to picking a state feedback gain ^ such that the 
quadratic performance index 


n*0 

Is minimized subject to 


E [x^ (n) 2, *(n') + 2 u^ M x(n) + u^ (n) R u(n) ] 


(7.1) 


x(n+'l) ■ A x(n) + ^ ]i(n) 

” n « 0,1,2, . . . (7.2) 

where £ > M, R are some weighting matrices having the property that 

positive definite, and a ^ 1 Is the stability factor. Frovlded 


r** -Tn 
M R 


that the system is controllable and cost observable, then the optimal 
state feedback gain is given by [Sa 2] 


% ■ (R + Sx <1^ J5o A + M) 


(7.3) 
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f^ere ^ is the ualque symnetrlc positive definite solution of the 
discrete RPDS Riccati equation 

^ K A + ^ - (B^ K A + M)^ K B) (b’' K A + M) ] 

•XL ~ ~ci ”“a~ -o— - -a— — 

(7.4) 

The eigenstructure and robustness properties for this class of 
regulators have not been treated in the literature thus far. In view 
of the similarity in structure between the continuous time and the 
discrete time RPDS problem, we expect the eigenstructure characterization 
of discrete time RPDS to follow from that of continuous time RPDS given 
in Chapter III. 

The generality of the framework for robustness analysis developed 
in [Le 1] (which we have extensively applied in Chapter IV) also allows 
us to characterize the robustness properties of discrete time RPDS in 
terms of the minimum singular value of its return difference and Inverse 
return difference. matrices. The derivation of the respective robustness 
theorems are based on the discrete multivariable Nyqulst theorem, which 
is essentially identical to Theorem 4.1 with the Nyqulst contour replaced 
by the unit circle centered at the origin. Based on such a framework, 
the effect of a on the robustness properties of discrete time RPDS can 
be studied in a fashion similar to that of Chapter IV. • 

A significant difference between the continuous time and discrete 
time RPDS occurs when one tries to solve the discrete time RPDS problem 
using only time- invariant weighting matrices Ccompare with corollary 2.5). 
It can be readily shown by algebraic manipulation of the discrete Riccati 
equation (7.4) that an appropriate choice of constant matrices £, M and 
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R for such purpose is given by 


a - i + (a^ - 1)A^ A 

(7.6) 

M - M + (O^ - l)t^ K A 

(7.7) 

R - R + (a^ - 1)B^ K B 

— — — . 

(7.8) 


t^ere ^ is the unique positive definite solution of the algebraic Rlccatl 
equation. The Intrestlng point to observe Is the expression for M as 
given by (7.7). M turns out to be a nonzero matrix even if the cross- 
weighting matrix M in the time-varying cost furctional (7.1) is chosen to 
be zero. In the dual problem of FPDS design, this Implies the need to 
use correlated process noise and observation noise that are second order 
stationary. A proper explanation of the above observation is not obvious. 
Understanding of this problem is probably Important to a better appreciation 
of the robustness properties of discrete time RPDS, for it is well known 
in the case of continuous time optimal regulators that the use of cross- 
weighting matrix between u(t) and x(t) can lead to deterioration of the 
stability margins. 

(ill) Alternative formulation of the Regulators with Prescribed 
Damping Ratio (RFDR) problem. 

In Chapter III of this report, the RFDR problem is formulated 
and solved as a special case of the RPDS problem. It is of Interest to 
know if there exists a more direct approach to formulate and solve the 
RPDR problem. 
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